Distribution Functions of tlie Nucleon and Pion in tlie 
Valence Region 



Roy J. Holt* 

Physics Division, Argonne National Laboratory, 
Argonne, Illinois 60439, USA 

Craig D. Roberts^ 

Physics Division, Argonne National Laboratory, 
Argonne, Illinois 60439, USA and 

Department of Physics, Peking University, Beijing 100871, China 

Abstract 

We provide an experimental and theoretical perspective on the behavior of unpolar- 
ized distribution functions for the nucleon and pion on the valence-quark domain; 
namely, Bjorken-x ^ 0.4. This domain is key to much of hadron physics; e.g., a 
hadron is defined by its flavor content and that is a valence-quark property. Fur- 
thermore, its accurate paramctrization is crucial to the provision of reliable input 
for large collider experiments. We focus on experimental extractions of distribution 
functions via electron and muon inelastic scattering, and from Drell-Yan interac- 
tions; and on theoretical treatments that emphasize an explanation of the distri- 
bution functions, providing an overview of major contemporary approaches and 
issues. Valence-quark physics is a compelling subject, which probes at the heart of 
our understanding of the Standard Model. There are numerous outstanding and 
unresolved challenges, which experiment and theory must confront. In connection 
with experiment, we explain that an upgraded Jefferson Lab facility is well-suited to 
provide new data on the nucleon, while a future electron ion collider could provide 
essential new data for the mesons. There is also great potential in using Drcll-Yan 
interactions, at FNAL, J-PARC and GSI, to push into the large- a; domain for both 
mesons and nucleons. We argue furthermore that explanation, in contrast to mod- 
eling and paramctrization, requires a widespread acceptance of the need to adapt 
theory: to the lessons learnt already from the methods of nonperturbative quantum 
field theory; and a fuller exploitation of those methods. 
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I. INTRODUCTION 



From the first deep inelastic scattering experiments (Bloom et ai, 1969; 
Breidenbach et ai, 1969) and the advent of the parton model (Bjorken, 1967; 
Bjorken and E.A.Paschos, 1969; Feynman, 1969) there has been a tremendous effort 
to deduce the parton distribution functions of the most stable hadrons - the proton, 
neutron and pion. The long sustained and thriving interest in these structure functions is 
motivated by the necessity to understand hadron structure at a truly fundamental level. 
While it is anticipated that quantum chromodynamics (QCD) will provide the explanation 
of hadron structure, a quantitative description of hadrons within QCD is not yet at hand. 
An explanation and prediction of the behavior of parton distribution functions in the 
valence-quark region;^ viz., Bjorken-x ^ 0.4, poses an important challenge for QCD and 
related models of hadron structure. This is the domain on which the transition takes place 
from deep inelastic scattering, with incoherent elastic scattering from numerous loosely- 
correlated small-x partons, to scattering from dressed-quarks that become increasingly well 
correlated as x — 1. We will focus our attention herein on the behavior of unpolarized 
structure functions on the valence-quark domain. 

The parton distributions are also essential to an understanding of QCD's role in nuclear 
structure. It has been known (Aubcrt et ai, 1983) for more than two decades that the par- 
ton distribution functions in a nucleus cannot simply be obtained by adding together the 
distributions within the constituent nucleons. This mismatch is the so-called EMC effect 
and, although concerted efforts have led to the identification of some of the ingredients nec- 
essary to an explanation, we still lack a completely satisfactory understanding of the nuclear 
dependence of parton distribution functions. This hampers us enormously. For example, an 
accounting for the EMC effect in light nuclei is key to extracting the neutron structure func- 
tion. That is a necessary precursor to a veracious determination of the differences between 
the parton distribution functions of the light-quarks. 

An accurate determination of the pointwise behavior of distribution functions in the va- 

^ As will become clear, for a given hadron in its infinite momentum frame, a parton's Bjorken- a; value 
specifies the fraction of the hadron's momentum carried by this parton. There is no unambiguous beginning 
to the valence-quark domain. We choose x > 0.4 because thereafter the gluon distribution can be said to 
be much smaller than the valence u-quark distribution in the proton. 
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lence region is also important to very high-energy physics. Particle discovery experiments 
and Standard Model tests with colliders are only possible if the QCD background is com- 
pletely understood. QCD evolution, apparent in the so-called scaling violations by parton 
distribution functions, entails that with increasing center-of-mass energy, s, the support at 
large-x in the distributions evolves to small-x and thereby contributes materially to the 
collider background. 

Deep inelastic scattering (DIS) of electrons from protons and bound neutrons at the Stan- 
ford Linear Accelerator Center (SLAC) led to the discovery of quarks. These experiments 
observed more electrons scattering with high energy at large angles than could be explained 
if protons and neutrons were uniform spheres of matter.^ In the approximately forty in- 
tervening years, electron DIS has played a central role in measuring structure functions, 
and SLAC has been especially effective in mapping the proton's structure functions in the 
valence- quark region. 

More recently, muon scattering experiments - performed by the European Muon Col- 
laboration (EMC), the New Muon Collaboration (NMC) and the Bologna-CERN-Dubna- 
Munich-Saclay Collaboration (BCDMS) - have contributed to our store of information. 
There have also been a substantial number of neutrino scattering experiments. However, 
they have generally used nuclear targets rather than pure hydrogen targets. Drell-Yan ex- 
periments have been effective at measuring the anti-quark distributions in the proton and 
nuclei. Putting all this together, it can be said that as a consequence the proton structure 
function is extremely well known, at least for x ^ 0.7. 

Herein we will discuss the status of charged-lepton and Drell-Yan experiments, as well 
as prospects for new experiments at, e.g., the Thomas Jefferson National Accelerator Fa- 
clity (JLab) and FermiLab (FNAL). On the other hand, although a number of experiments 
at the Conseil Europeen pour la Recherche Nucleaire (CERN) and Deutsches Elektronen- 
Synchrotron (DESY) have also provided measurements of the proton structure function, 
many of these efforts focused on the low-x behavior and hence they will not be discussed. 

An experimental determination of the neutron structure function at high Bjorken-x has 
proved especially troublesome, the main reason being that most of our information about the 

^ In both method and results this scries of experiments, conducted from 1966-1978 and for which Taylor, 
Kendall and Friedman were awarded the 1990 Nobel Prize in Physics, was kindred to that which led 
Rutherford to discovery of the nucleus in 1911. 
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neutron structure function is obtained from DIS experiments on a deuteron target, for which 
the nature of the EMC effect is simply unknown. Consequently, inference of the neutron 
structure function from proton and deuteron measurements is model-dependent. We will 
treat this topic in some depth and, in addition, canvass prospects for the future. 

On the theoretical side the challenge is first to parametrize and ultimately to calculate 
the parton distribution functions. We write ultimately because the distribution functions 
are essentially nonperturbative and therefore cannot be calculated in perturbation theory. 
Thus, absent a truly-accurate, quantitative and predictive nonperturbative tool, the QCD 
calculation of parton distribution functions remains an alluring but distant prospect. Herein 
we will provide both an historical perspective on past attempts at a theoretical interpre- 
tation and an overview of recent progress towards this goal. Notably, today, there are no 
computations of the pointwise behavior of the nucleon's valence-quark distribution function 
that agree with the predictions of the QCD parton model. 

There has, on the other hand, been much success to date with parametrization. A number 
of independent groups have analyzed and performed a global fit to the vast body of extant 
DIS and Drell-Yan nucleon structure function data. In this approach the physical processes 
are factorized into the product of short- and long-distance contributions. The so-called 
short-distance parts are calculable in perturbation theory, whereas the long-distance parts 
are determined by the parton distribution functions. In this way, via the operator product 
expansion in QCD, scattering processes involving hadrons are connected to subprocesses 
involving partons. The utility of this approach is grounded on the fact that so long as 
factorization is valid, the parton distribution functions are universal; namely, all hadron- 
level interactions for which factorization applies are described by the same small body of 
parton distribution functions. This being true, the nonperturbative problem is reduced to 
that described above: instead of calculating all cross-sections from scratch, one has merely to 
calculate the distribution functions. Herein we will introduce the modern parametrizations. 

At its simplest, the nucleon is a three-body problem. As ostensibly a two-body problem, 
it would appear theoretically simpler to calculate properties of the pion and hence the parton 
distributions therein. However, the pion is in fact both a bound state and the Goldstone 
mode associated with dynamical chiral symmetry breaking (DCSB) in QCD. This amplifies 
the importance of understanding its properties. However, it also significantly complicates 
the calculation of the pion's distribution functions and places additional constraints on 
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any framework applied to the task: the pion simply cannot veraciously be described as a 
constituent quark-antiquark pair. Instead, key features of nonperturbative quantum field 
theory must be brought to bear. 

Although it is not presently possible to perform deep inelastic scattering on a free pion, 
the pion structure function has been measured in pionic Drell-Yan experiments at CERN 
and FNAL. Owing to predictions from the QCD parton model, the form of the pion's 
valence-quark distribution function provides a stringent test of our understanding of both 
QCD and QCD-based approaches to hadron structure. The FNAL Drell-Yan experiment 
(Conway et al, 1989) measured the pion structure function up to high Bjorken-x, where 
the QCD predictions were expected to be realized. However, they were not manifest, so 
that a very puzzling discrepancy remains at this time. The status of structure function 
measurements of the pion as well as prospects for future experiments at JLab and a possible 
future electron ion collider (EIC) will be presented, as will a full description of the theoretical 
perspective. 

The purpose of this article is not to present a comprehensive review of all the issues 
involved in extracting structure functions from high energy data, as in (Sterman, 1995), nor 
an exhaustive digest of theoretical developments in the field. Rather the focus will be on 
structure functions in the valence region, which have chiefly been extracted from electron and 
muon inelastic scattering, as well as from Drell-Yan interactions. Moreover, we will focus on 
the theoretical treatments that emphasize an explanation of these structure functions and 
summarize the major contemporary approaches and issues. To end this Introduction, a brief 
remark on notation: unless denoted otherwise, when addressing a distribution function, it 
will be that associated with the proton; e.g., ^^^,(x) will mean the proton's valence -u-quark 
distribution. 

II. NUCLEON STRUCTURE FUNCTIONS FROM ELECTRON AND MUON 
SCATTERING 

A. Kinematics 

The Feynman diagram for lepton inelastic scattering from a hadron is given in Fig. II. 1. 
Practically, the leptons in a laboratory scattering process can be electrons, muons or neu- 
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FIG. II. 1 Feynman diagram representing deep inelastic lepton scattering from a hadron. 

trinos, as illustrated in the figure, and the current can be a photon, a charged weak current 
W^, or a neutral weak current Z° of four momentum transfer, q. The incident lepton has 
incoming and outgoing four momenta of k and k', respectively, so that the momentum- and 
energy-transfers are given by 

q =k-k', (III) 
u = E-E' , (II.2) 

where the incident momentum and mass of the hadron are given by p and M, while the 
outgoing debris after the inelastic event has a momentum p'. 

The invariants in the scattering process are: the square of the four momentum transfer; 
the square of the invariant mass; p ■ q, which is related to the energy transfer in the target 
rest-frame; and the square of the total energy in the center-of-mass frame. The square of the 
four-momentum transfer depends on the energy of the beam and scattered lepton as well as 
the lepton scattering angle: neglecting the lepton's mass, 

Q2 := -g2 = ^EE'sin^^ . (II.3) 
8 



The square of the invariant mass, W, is given by 



= {p + qf = M^ -Q^ + 2p-q = M^ + -Q\ (11.4) 

X 

where M is the mass of the hadron and 

X = (II.5) 
2p ■ q 

is the Bjorken variable. The square of the total energy in the center of mass is given by the 
Mandelstam variable, s: 

s = {k + pf (II.6) 

A final, useful invariant is 

which, in the target rest-frame, is a measure of the fractional energy loss by the incident 
lepton. 

The Bjorken limit is theoretically defined as 

Q^— j-cxD, 2p ■ q oo ,x = fixed . (II-8) 

In this limit, x can be shown to equal the fraction of the hadron's momentum carried by the 
struck quark. Empirically, in order to avoid complications associated with the production 
of hadron resonances, the Bjorken scaling regime is explored via 

W^^oo, g^^oo, x = ^^—^^;p^-^ = fixed. (II.9) 

In principle, the deep inelastic approximation is only good when 

However, it is empirically well established that scaling in deep inelastic scattering (DIS) sets 
in at tractable values of W and Q^. As a practical matter, evaluations of data typically 
impose the requirement that > A GeV^ and W > 3.5 GeV. As discussed and explored, 
e.g., in (Accardi et ai, 2009; Malace et ai, 2009), data at high Bjorken x but not meeting 
these kinematic requirements would be subject to large target-mass corrections, which are 
kinematic and owe to binding of partons in the hadron. Before discussing these topics in 
more detail, the cross section for inelastic electron scattering will be presented. 
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B. The cross section for charged lepton scattering 



The cross section for inelastic charged-lepton scattering from a hadron via a one-photon 
exchange process, as indicated in Fig. II. 1, is well known and given by 

(.1.11) 

where L^^^ is the leptonic tensor and W'^'^ is the hadronic tensor. The leptonic tensor has 
the following straightforward dependence on the kinematic variables: 



k^k'" + k^k'^ + g^" (m^ - {k ■ k'))] , (11.12) 



where m is the lepton's mass. Conserving current and parity, the hadronic tensor for a 
spin-| charged particle depends on two structure functions, Wi{h',Q'^) and W2{i^,Q'^); viz., 

^ + ^) + ^ (v. - [v. - ^,.) . (11.13) 

Contracting the leptonic tensor with the hadronic tensor yields the well-known cross section 
of the form: 



(11.14) 



dQdE" ' ' ' 

The quantities Wi and W2 in Eq. (11.14) are often written in terms of dimensionless 
structure functions, Fi and F2, for the nucleon: 

MW^iu,Q^) =Fi(x,g2), (11.15) 
uW2iiy,Q^) =F2(x,Q2). (11.16) 

In turn, these structure functions are frequently written in terms of the transverse {Ft) 
and longitudinal {Fl) structure functions, which correspond to absorption of transverse and 
longitudinal virtual photons, respectively: 

FT(x,g2) = 2xFi(x,g2), (11.17) 
Fl{x,Q^) = F2{x,Q^)-2xFi{x,Q^). (11.18) 

They can, in principle, be measured separately by performing a Longitudinal/Transverse 
(L/T) separation, as will be discussed in Sec. II. H. The most intuitive physical picture of 
the interaction in the deep inelastic scattering regime is provided by the parton model. 
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C. The parton model 



It was discovered at SLAC in the late 1960s that the structure functions appeared to 
scale; i.e., their evolution with x is nearly independent of over a very large kine- 
matic range. These important points are well documented, not only in the original pa- 
pers (Bjorken, 1967; Bjorken and E.A.Paschos, 1969; Bloom et al, 1969; Breidenbach et ai, 
1969; Fcynman, 1969) but also in the Nobel lectures summarized in (Friedman, 1991; 
Kendall, 1991; Taylor, 1991), and here will only be reviewed briefly. 

The observed scaling feature means that the cross-section is not separately a function 
of the two kinematic variables, energy transfer, z/, and . Instead, the cross-section's 
behavior can be expressed in a dependence on only one scaling variable; namely, Bjorken x. 
Indeed, in the Bj or ken-limit, Eq. (II. 9), the proton structure functions Wi and W2 assume 
a form consistent with elastic scattering from a point fermion; viz., 

2MW,M^) ^ 2F,(.) ~ (1 - J^) . (11.19) 

^W2(^, Q-') ^ F2{x) ~s(^l- ^Jl-^ , (11.20) 

where x is given by Eq. (II. 5). This near-dependence on only a single dimensionless variable 
is referred to as Bjoken scaling. Assuming this to be the case, then the expression for the 
cross section, Eq. (11.11), becomes 



da .azE 



2 77/2 



-{E,E\e)=A 



e 



F^ix) 2 e 2F,{x) . 2 9 

cos — sm - 

V 2 M 2 



(11.21) 



dndE" ' ' ' 

Capitalizing on these results and assumptions, further analysis leads to the relation 
(Callan and Gross, 1969): 

F2{x) = 2xFi{x). (11.22) 

This "Callan-Gross relation" is observed to be approximately valid, a feature that is often 
cited as evidence that the pointlike constituents of the proton; i.e., the quarks, are spin-| 
degrees-of-freedom. 

An intuitive understanding of the observed scaling behavior is provided by the quark- 
parton-model (QPM) or QCD parton model. In this model the electron collides with a 
collection of partons within the hadron. Owing to Lorentz contraction of the hadron at ex- 
tremely high momentum, these partons are practically frozen in time during the extremely 
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brief collision process. Moreover, the probability of finding two partons near enough to inter- 
act with each other drops as l/Q^. Thus, initial and final state interactions are suppressed 
in the DIS regime. Within this model, the photon interacts incoherently with individual 
partons and the cross section depends in a simple way on the probability of finding a quark 
of flavor "i" with fraction x of the proton's momentum. These properties are realized in the 
following simple expression for the electromagnetic structure function 

F2{x) = 2xF,{x) = Y.e^x + , (11.23) 

i 

where the Cj are the charges of the individual quarks and qi{x) are the probability densities 
for finding quark-i with fraction-x of the hadron's momentum. Assuming that only the 
relatively light quarks {u, d, s, c) and their respective antiquarks contribute to the proton 
structure function, then the proton structure function can be written as 

14 1 r - 1 

-F^{x) = - \u{x) + u{x) + c{x) + c{x)] + - d{x) + d{x) + s{x) + s{x) . (11.24) 

X 9 9 L J 

In Fig. 11.2 we depict the contributions of these quark flavors to the structure function 
in the valence region. The CTEQ6L leading order (LO) evaluation (Martin et al, 2004a; 
Stump et ai, 2003), a modern global evaluation of existing data, was used for illustration. 
The distribution of up and down quark flavors is fixed primarily by electron and muon 
scattering from the proton and neutron, while the distribution of the u and d quark sea is 
largely determined by Drell-Yan processes. The strange and charm quark sea is determined 
by neutrino interactions. The w-quark distribution is the dominant contribution to the 
proton structure function, while the d- and anti-u-quarks are significant. The charm and 
strange quark distributions are small and are not easily distinguishable in the figure. We 
note that semi-inclusive electron scattering from the proton and deuteron is emerging as 
an important tool for determining the strange quark distribution, thereby providing an 
important complement to the method of neutrino scattering from complex nuclei as a probe 
of non- valence proton structure functions. Of course, the u and d quark distributions in the 
valence region are particularly important for testing descriptions of nucleon structure, and 
it is therefore vital to obtain accurate and precise data, and a reliable flavor separation in 
this region. 

A primary goal of modern theory is the computation of valence-quark distribution func- 
tions, using QCD- motivated or QCD-based models and nonperturbative methods in QCD. 
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FIG. II. 2 (Color online) Quark flavor dependence of the proton structure function according to the 
CTEQ6L (leading order) evaluation at a scale = lOGeV^ (Martin et ai, 2004a; Stump et ai, 
2003). We plot the quantities ef xqi{x) in order to emphasize the contribution to the structure 
function from from the various quark flavors. 

Naturally, valence quarks have no strict empirical meaning because experiment cannot read- 
ily distinguish a valence- from a sea-quark. However, valence-quark distributions for the u- 
and (i-quarks in the nucleon can be defined as the following flavor-nonsinglet combinations: 

Uv{x) := u{x) — u{x) , (11.25) 
d^,{x) := d{x) - d{x) . (11.26) 



Conservation of charge leads to important normalization conditions: 

j dxUy{x) = 2, [ dxdy{x) = 1. (11.27) 

Naturally, 







dxSy{x) = Q, (11.28) 

JO 



with the same result for all heavier quarks. However, it does not follow that Sv{x) must be 
identically zero. 
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Consistent with nomenclature, the combinations in Eqs. (11.25) and (11.26) have nonzero 
values of various flavor quantum numbers, such as isospin and baryon number. Singlet 
distributions also play a role and are constituted from a sum of quark distribution functions: 
S(x) = T,i[qi{x) + qi{x)]. Under QCD evolution, discussed below, singlet distributions mix 
amongst themselves and with the gluon distribution, but the nonsinglet distributions do 
not. 

D. QCD and scaling violations 

The proton structure function in the valence region is shown in Fig. II. 3 as a function 
of Bjorken x and for three values of for the CTEQ6L LO evaluation. Although the 

range illustrated is relatively large, from 2 GeV^ to 100 GeV^, the structure function 
does not exhibit a strong dependence on Q^. This result validates the parton model as an 
approximation to a QCD description of the structure function. Of course, the structure 
functions have a dependence on Q^, indicating the necessity of QCD rather than a parton 
model to describe the data. This relatively small dependence is known as a scaling 
violation. 

The validity of the QPM led to the idea of factorization in deep inelastic scattering. In 
factorization, the cross section is separated into two distinct contributions: a short-distance 
part, described by perturbative QCD (pQCD); and a long distance part. The Q^- dependence 
driving scaling violation originates in the pQCD part and owes to the same phenomena that 
generate running of the strong coupling constant, a^, and hence asymptotic freedom. This 
factorization for the moments of the structure functions can be shown rigorously using the 
operator product expansion and the renormalization group equations. Alternatively, a dis- 
tribution function where the moments can be written as products can be expressed in terms 
of a convolution of two distributions. Naturally, the long-distance part is only calculable 
using nonperturbative methods and presents a challenge for modern hadron theory. 

At small distances or, equivalently, large momentum transfers in QCD, the effec- 
tive coupling constant becomes small. This feature is known as asymptotic freedom 
(Gross and Wilczek, 1974; Politzer, 1974). In terms of the renormalization scale, Aqcd, 



14 



0.5 - 




0.0, 







0.2 



0.4 



0.6 



0.8 



X 

FIG. II. 3 (Color online) Structure function of the proton at three values of Q^, computed from 
the CTEQ6L leading-order evaluation. 

the strong coupling constant has the perturbative form: 



where: t = In {Q"^ / Aq^j^); /3o = H — 2nf/3, f3i = 102 — 38nj/3 and Uf is the number of quark 
flavors involved in the process. From Eq. (11.29) it is apparent that — as — oo, 
so that a parton-model limit is contained within QCD; e.g., as will become clear below, 
Eq. (11.23) is obtained from Eq. (11.37) in this limit. 

Plainly, the deflnition of short- and long-distance phenomena is ambiguous, and associ- 
ated with a separation of the cross-section into these two contributions is a factorization 
scale; namely, a mass-scale chosen by the practitioner as the boundary between hard and 
soft. If that scale is large enough, then the hard part of the scattering cross-section can be 
calculated in pQCD. This, however, introduces another scale, associated with renormaliza- 
tion. For convenience and economy, the factorization and renormalization scales are usually 
chosen to be the same. It follows, in addition, that factorization is scheme-dependent be- 
cause the choice of renormalization scheme implicitly specifles a division of the flnite pieces 
of the cross-section into those that are retained in the hard contribution and those under- 
stood to be contained in the soft piece. The part identifled as owing to long-distance effects 




(11.29) 
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is basically the parton density distribution function. 

In fitting data to construct PDFs, it is usual to use the MS (modified minimal subtrac- 
tion) renormalization scheme whereas, in regard of factorization, there are two commonly 
used schemes. The DIS scheme (Altarelli et ai, 1978) was designed to ensure there are no 
higher-order corrections to the expression for the F2 structure function in terms of the quark 
PDFs; i.e., all finite contributions are absorbed into the PDF. On the other hand, in the 
more widely used MS scheme (Bardeen et ai, 1978), in addition to the divergent piece, only 
the usual (ln47r — 7^;) combination is absorbed into q{x) and hence the expression for F2 
exhibits explicit 0{as) corrections. [See (Brock et ai, 1995) for more on these points.] 

The CTEQ6L (leading-order) structure functions shown in Fig. II. 3 were evaluated using 
the MS'-scheme (modified minimal subtraction). At leading-order there is no difference 
between factorization schemes. As we will subsequently see (Fig. II. 4), radiated gluons give 
rise to the scaling violation or Q^-dependence of the structure functions. They, and the 
quark distribution functions defined therefrom, then depend on both x and Q^. The 
dependence is now routinely described within the framework of next-to-leading-order (NLO) 
QCD evolution. 

An important and extremely useful feature of factorization is that a measurement of a 
structure function at relatively-low permits the calculation, through the use of pQCD, 
of the structure function at high Q^. In leading order, a set of integro-differential equations, 
now known as the DGLAP evolution equations (Altarelli and Parisi, 1977; Dokshitzer, 1977; 
Gribov and Lipatov, 1972; Lipatov, 1975), are used for this purpose. Intuitively, one may 
think that as increases, a parton can sometimes be resolved into two partons; e.g., a 
quark can split into a quark and a gluon, or a gluon into two gluons - see Fig. II. 4 for a 
graphical representation of the leading order DGLAP equations. The two resolved partons 
then share in the fraction of the nucleon momentum carried by the initial unresolved parton 
at lower Q^. The parton distribution thus becomes a function of Q^. Such evolution of the 
parton distribution function (PDF) can be seen in Fig. II. 3. At the highest values of the 
structure function is shifted toward lower values of x. This work has been generalized to 
other reaction processes (Furmanski and Pctronzio, 1982). 

We note that in principle the domain of "relatively-low Q^" means, nevertheless, 
^ ^QCD- Notwithstanding this, in some empirical determinations of PDFs through 
data fitting [see Sec. V.C] and in the application of models to their calculation [see Sec. VI], 
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FIG. II. 4 Graphical representation of the leading order DGLAP equations. 

evolution is sometimes applied from low ~ AAqqj^. In such cases an interpretation of 
the \ow-Q^ distribution as a parton model density is questionable owing to the probable 
importance of nonperturbative corrections. However, contributions from such corrections 
are suppressed by DGLAP evolution to larger-Q^ and hence one anticipates that in these 
applications there is always a whereafter the desired interpretation becomes valid. 

As stated above, the evolution equations enable the accurate calculation of the PDFs at 
a general Q^, provided they are known at another scale, Q^, so long as pQCD is a valid 
tool at both scales. At next-to-leading order the QCD evolution equations have the form 
(Herrod and Wada, 1980) 

dt 

where the convolution is defined as 



PNS±®qNS±{x',Q^), (11.30) 



C®f:=f ^Ciy)fi-), (11.31) 

Jx V V 



and 

dt 



^'^M=p® ^'^M. 11.32 

In these formulae, qNS- = u—u, d—d and qNS+ = {u+u) — {d+d) or {u+u) + {d+ d) — 2{s + s) 
are the non-singlet PDFs, while S = X)i(9i + qi) is the singlet combination. The splitting 
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functions are given by 



Pns± = ^Pl{x) + {^)'p'ns±{^) , (n.33) 

p = ^P\^) ^ {^)' P\^) ^ (n-34) 

where 

= I ^'^'^ 2n/P4 I ^^^^^^ 

The leading order DGLAP equations are obtained simply by eliminating the a1 terms in the 
equations above. 

It will be apparent from Eqs. (11.30) and (11.31) that x = 1 is a fixed point; namely, that 
the X = 1 value of each distribution function is invariant under evolution - it is the same at 
every value of the resolving scale Q^. This is because the right-hand-side of each evolution 
equation vanishes at a; = 1. Naturally, when Bjorken-x is unity, then q^ + 2p-q = and hence 
one is strictly dealing with the situation in which the invariant mass of the hadronic final 
state is precisely that of the target: W = M in Fig. 11.1; viz., elastic scattering. The structure 
functions inferred experimentally in the neighborhood of x = 1 are therefore determined by 
the target's elastic form factors, which all vanish in the limit oo. This fact per se is not 

interesting. However, the rate at which distribution functions vanish with x is, because it 
can lead to nonzero renormalization-scale-invariant distribution-function-ratios at x = 1. For 
example, the value of limx^idy{x)/uv{x) in the nucleon is an unambiguous, scale invariant 
feature of QCD and hence a discriminator between models. [For a concrete illustration of this 
point, see also Fig. VI. 42, which depicts {x)/u'^{x).] This is an important consequence, 
and one much discussed and explored in experiment and theory. (See Sees. 11. G and VI.) 

We reiterate, however, that the neighborhood of x = 1 poses difficulties for experiment, 
and for theory, too, owing, e.g., to the importance on this domain of target-mass corrections 
and higher- twist contributions (explained below). Whilst the equivalence between x = 1 
and elastic scattering cannot be avoided, one might choose to approach the limit obliquely. 
For example, noting from Eq. (11.9) that 

Q^»wj-M^ {W-M){W + M) 

one would ideally choose a path such that both and W — M are simultaneously kept 
as large as kinematically possible. This is experimentally very difficult, as discussed in 
connection with Fig. II. 7. Future experimental possibilities are outlined in Sec. IV. 
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FIG. 11.5 (Color online) Comparison of the proton's u{x) distributions in LO, NLO and NNLO at 
g2 = 10 GeV^ 

A great triumph of QCD is the good agreement between the evolved structure function for 
the proton and experiment over many orders of magnitude in Q^. Of course, the structure 
function should be calculated at the same order as the PDFs. The equation for the structure 
function at NLO in the DIS scheme is given by 



where g{x,Q^) is the gluon distribution and the Wilson coefficients, {Ci^g}, are given in 
(Furmanski and Petronzio, 1982; Gliick et al, 1995a). 

During the past several years rather complete NNLO fits (Martin et ai, 2007) in the MS 
scheme have been performed for the data. These fits have made use of NNLO splitting func- 
tions (Moch et ai, 2004; Vogt et al, 2004). The LO, NLO and NNLO u-quark distribution 
functions are plotted at a scale of = 10 GeV^ in Fig. II. 5 for comparison. Naturally, the 
structure functions themselves are measured quantities and do not depend on the order of 
the fit. In the valence region, the largest difference occurs between LO and NLO as expected. 
There is some evidence (Yang and Bodek, 2000) from early NNLO analyses that evolution 
at NNLO largely offsets the effect of high twist in an NLO analysis on the large-x domain. 




(11.37) 
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However, a more recent analysis (Blumlein et al, 2007) performed up to N3L0 contradicts 
this finding. 

E. High twist effects and target mass corrections 

As we have seen, the interactions between partons at short distances that give rise to scal- 
ing violations are reasonably well described in a pQCD approach through the Q^-evolution 
equations. Probes of a hadron at intermediate values of might expose correlations among 
the partons. Thus far, we have only considered the case in DIS where structure functions are 
governed by incoherent scattering from the partons. As we move to lower where nucleon 
resonances could dominate or to very high x where the elastic scattering limit and complete 
coherence dominates, then correlations between the partons must be taken into account. 
Ultimately a theoretical treatment of these correlations demands an understanding of the 
way that quarks and gluons bind to form hadrons. Although perturbative QCD cannot 
explain binding effects, it does permit one to construct a model of power corrections to the 
perturbative result. It was found (Ellis et a/., 1982, 1983) that the first power correction, 
l/Q"^, to the structure function is governed by the intrinsic transverse components of the 
parton momentum. 

The operator product expansion is normally used to order contributions to a DIS structure 
function according to their twist, r, which is defined as the difference between the naive 
mass-dimension of an operator and its "spin". Here "spin" means the number of vector 
indices contracted with the configuration-space four- vector, z. For example (Sterman, 1995), 
consider the expansion of a dimension-six current-current correlator, typical of DIS, 

oo oo 

r{z)UO) a,n(^^/x') z^^... 0:^...,J;^ = 0,/i2) , (n.38) 

r=2 n=0 

where yU. is the factorization scale. One is interested in the behavior of the coefficient 
function Cr,n(^^, A*^) = {z'^)~'^'^^'^ C{z'^fi'^) because the more highly singular is this function 
on the light-cone, the more important is the associated operator at large Q^. Now, suppose 
the operator O^^ unif^'^) dimension do„, then, in order to balance dimensions on both 
sides of Eq. (11.38), one must have cic = 3 — {do„ — n)/2, since each factor of z on the 
right-hand-side contributes mass dimension negative-one. Plainly, the largest value of dc is 
obtained with the smallest value of twist r = do„ — n, and the associated operators are those 
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FIG. II. 6 (Color online) Empirical nonsinglet higher-twist coefficient, denoted Cht, as a function 

of Bjorken-x, for the proton (left panel) and deuteron (right panel), as obtained when the twist-2 
contributions are treated at NLO, NNLO and N3L0. [Figure adapted from (Blumlein et al, 2007).] 

which dominate the large-Q^ behavior of the correlator. As Eq. (11.38) indicates, numerous 
operators with different dimensions have the same twist and all such operators are associated 
with the same degree of light-cone singularity. 

As evident in the example above, the leading- twist contributions in DIS are twist-2. In 
unpolarized DIS, the higher twist components; i.e., twist = 4, 6, . . . , are suppressed by l/Q^, 
l/Q^ , . . . , respectively. Thus, in general a structure function should be expressed in the 
form (Alekhin et al, 2004) 



F2{x,Q'') = F,^^{x,Q'') + 



0{1/Q'), 



(11.39) 



where Fg^"^ refers to the leading- twist part. As a practical matter, the magnitudes of higher- 
twist terms are generally unknown and the estimates are therefore somewhat controversial. 
Before one can claim discovery of higher twist components, it is essential that these be 
distinguished from target mass corrections and Q^-evolution effects. 

Several analyses have shown that higher-twist coefficients can become quite large at 
high-x and relatively low Q^. An analysis (Virchaux and Milsztajn, 1992) of BCDMS and 
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SLAC data indicated that the twist-4 coefficient is sizable for the proton and deuteron. 
Here, however, the inclusion of target mass corrections reduced the magnitude of the twist-4 
coefficient. This finding was supported in an analysis that included NNLO (Yang and Bodek, 
2000). It was found that the NNLO contribution can partially compensate for a large higher- 
twist effect. More recently, an analysis (Blumlein et ai, 2007) of world data indicated that 
the twist-4 coefficient is sizable for NLO, NNLO and N3L0 for the proton and deuteron. 
While increasing the order of the QCD analysis of the data reduces the twist-4 coefficient, the 
coefficient is significant at high x. This is displayed in Fig. II. 6, which plots an empirically 
determined correction defined via 

F,{x, Q') = OTMc[Ft'''-\^, Q')] (l + ^i^I^^) , (11.40) 

where Otmc describes inclusion of target mass corrections of the twist-2 contributions to 
the structure function. At present, there exists no statistically significant determination of 
the magnitude of the twist-6 coefficient. 

In connection with the size of higher-twist contributions, it is relevant to consider their 
impact on Bloom-Gilman duality (Bloom and Oilman, 1971). This is an expectation that the 
nucleon structure function, when measured in the region dominated by low-lying nucleon 
resonances, should follow a global scaling curve that is determined by high energy data. 
NB. Whilst it has been observed that the standard application of DGLAP evolution to deep 
inelastic structure functions can appear to be inconsistent with Bloom-Gilman duality at 
fixed W , this apparent conflict is resolved if one takes into account the fact that the struck 
quark is far-off shell in the x — )■ 1 domain (Brodsky, 2005; Lepage and Brodsky, 1980). 

Underlying Bloom-Gilman duality is the assumption of quark-hadron duality; namely, a 
belief that one obtains equivalent descriptions of physical phenomena irrespective of whether 
one uses partonic or hadronic degrees-of- freedom (Melnitchouk et al, 2005; Shifman, 2000). 
While this assumption is probably correct in principle, in comparisons between computations 
that use finite-order truncations, it can and is violated in practice. The presence of large 
higher-twist contributions can therefore entail an apparently strong violation of Bloom- 
Gilman duality. 

Experiments, particularly at the highest values of x, might not be performed at suffi- 
ciently large momentum transfer to avoid high-twist effects and target-mass corrections. A 
pQCD treatment is truly valid only in the limit that the mass of the hadron is negligible 
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in comparison with the dominant scale, . The target-mass corrections are kinematic and 
owe to binding of partons in the hadron. Nachtmann first pointed out (Nachtmann, 1973) 
that at finite and non-vanishing target mass, the scahng variable should be the fraction, 
^, of the nucleon's light front momentum carried by the parton; viz.. 



Naturally, this Nachtmann variable reduces to Bjorken-x for M^/Q^ — )■ 0. 

It was subsequently shown (Georgi and Politzer, 1976) that at leading-order, if the quan- 
tity x^M^/Q^ is reasonably small, then it is straightforward to apply corrections for finite 
mass to the structure function. If the leading-twist structure function is denoted by F^^^ 
and the target-mass corrected structure function by F^^^ , then one has the expressions: 



Q^) = + dx'Frix) + 2^ £ d."Ft\x") , (11.42) 

F™(x,g^) = ^^^(0 + 6^ fdx'F^^{x') + l2^ [\x"Fr{x").{llA3) 

ti KcJ hi Jf KJ hi J x' 



A concern with the Nachtmann or the Georgi- Politzer approach is the so-called "threshold 
problem" It results from the fact that the maximum kinematic value of ^ is less than unity, 
which means that the corrected leading twist structure function does not vanish at x = 1. 

Target-mass corrections have been extended to NLO QCD (DeRujula et al, 1977), 
and generalized to include charge and neutral weak current deep inelastic scattering 
(Kretzer and Reno, 2004). De Rujula et al argued that the threshold problem could be 
solved by considering higher-twist effects. 

A recent interesting approach (Steffens and Melnitchouk, 2006) to the "threshold prob- 
lem" followed the original Georgi-Politzer approach but with one critical difference. Namely, 
the upper limit of the integrals in Eq. (11.42) becomes the upper limits of x' and x" that 
are permitted by kinematics at x = 1 rather than merely unity. In terms of existing data, 
this solution does not result in a large effect for > 2 GeV^, but can have a significant 
effect below 1 GeV^ where the underlying assumption of factorization becomes questionable 
anyway. 

In addition to target-mass corrections, there are also effects from nonzero current-quark 
mass (Barbieri et al, 1976). Such corrections will not be discussed herein because of our 
focus on high-x, whereat heavy quarks have little effect in nucleons and Goldstone bosons. 
These effects, in addition to those arising from target-mass corrections, are discussed in a 




(11.41) 
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FIG. II. 7 vs. X for three values of W . Measurements at very high x demand very high to 
remain in the DIS region. 

recent review (Schienbein et al, 2008). Furthermore, fresh work (Accardi and Qiu, 2008) 
has pointed out that changing the upper hmit of the integrals in the Georgi-Pohtzer equations 
leads to an abrupt cutoff at x = 1, and proposed a "quark jet" mechanism to give a more 
reasonable approach to the very large x behavior. 

F. The proton structure function 

Measurements of structure functions at very high values of x are extremely challenging. 
The main problems can be seen at a glance from Fig. II. 7. Here the is plotted from 
Eq. (II. 4) as a function of x. Typical evaluations of the parton distribution functions require 
that ly > 3.5 GeV. Clearly, the necessary to meet this demand is extremely high and 
impractical, at present, for x > 0.7. 

The first DIS measurements, performed at SLAG, were limited to a beam energy of 
24.5 GeV. Although this beam energy limit is a serious problem for very high x measure- 
ments, the SLAG experiments had good control of systematic errors. Relatively high lumi- 
nosities were possible, the incident electron energy was very well known and the spectrometer 
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FIG. II. 8 Proton structure function vs. Bjorken-x with W > 2.0 GeV. For ease of representation: 
the upper band of data were multipUed by a factor of 10 and correspond to the domain 25 GeV^ < 
q2 < 40GeV^; the next band are scaled by a factor of 3 and cover the domain: 15 GeV^ < < 
25GeV^; the third-lowest band are unsealed and cover: 7GeV^ < < 15GeV^; and the lowest 
band are scaled by a factor of 1/5 and cover: 2GeV^ < < 7GeV^. 

properties were relatively well understood. A consistent re-analysis of all the early SLAG 
experiments was performed in (Whitlow et ai, 1992), and the result is displayed in Fig. II. 8. 

After the importance of the early SLAG experiments was realized, the quest for addi- 
tional data, particularly at high Q^, led to muon scattering experiments at GERN and 
FNAL. At GERN a series of experiments was performed at muon energies between 90 and 
280 GeV by the NMG collaboration (Amaudruz et al, 1992), between 120 and 280 GeV by 
the BGDMS collaboration (Benvenuti et al, 1990b) and at 280 GeV by the EMG collabo- 
ration (Aubert et ai, 1987). In this case, the combination of very high beam energies and 
relatively low luminosities makes it difficult to obtain structure function data at very high 
X. The practical limit is about x = 0.75 as shown in Fig. II. 8. The FNAL E665 experi- 
ment provides data at even lower x than the GERN experiments since the beam energy was 
470 GeV. 

Data for various ranges from the BGDMS (Benvenuti et ai, 1990b), EMG 
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FIG. II. 9 (Color online) The proton structure functions vs. Bjorken-x with W > 2 GeV. The 
dashed black curve represents the CTEQ6L evaluation with = 10 GeV'^, while the solid red 
curve indicates the result obtained when target mass corrections are incorporated (Malace et al, 
2009). 

(Aubert et al, 1983, 1987), NMC (Amaudruz et al, 1992; Arneodo et al, 1997a) and SLAC 
(Whitlow et al., 1990) experiments are plotted in Fig. II. 8. The effect of the kinematic con- 
straints on the reach in x is evident throughout the relatively wide range of values for Q^. 
The SLAC data persist out to the highest values of x, but only for relatively low values of 



If the condition on W is relaxed to > 2 GeV, then data from SLAC exist up to x = 0.85 
as indicated in Fig. II. 9. Clearly, the CTEQ6L evaluation at the very highest values of x 
deviate from the data. However, if target mass corrections are taken into account, as shown 
by the dashed curve in the figure, then good agreement is restored. This illustrates that 
target mass corrections are indispensable in order to fully understand and extract parton 
distributions at very high x. 
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G. The neutron structure function 



Interest in the neutron structure function at very high x has burgeoned during the past 
decade. The neutron structure function, or at least the ratio of the neutron structure 
function to that of the proton, is beheved to be one of the best methods to determine 
the d{x)/u{x) ratio. A knowledge of the this ratio in the valence region would provide an 
important constraint on models for the nucleon (Brodsky et al, 1995; Isgur, 1999). For 
example, if one were to assume that a simple SU{6) flavor symmetric model for the nucleon 
were valid, then the d/u ratio would just be 1/2 for the valence quark distributions. In this 
model, the nucleon and A would be degenerate, so this is clearly not a correct description of 
the nucleon. Other interesting limits occur, for example, where the d-quark in the proton is 
"sequestered" in a pointlike scalar diquark. In this model, d/u\^^-^ = and F^/Ffl^.^^ = 
1/4. Furthermore, in a world where pQCD could be applied naively and scattering at large-x 
involves only quarks with the same helicity as the target hadron, then d/u\^^-^ = 1/5 and 
-^2"/-^2 L->i ~ 3/7 = 0.43. [An extensive discussion of these issues is provided in Sec. VI, 
where it is argued, e.g., that FJ'/Ffl^^-^ = 0.36 and only with reliable data at a; > 0.8 will 
one be empirically able to determine the x — t- 1 behavior.] 

The parton model is a good starting point for introducing the neutron structure function. 
Under the assumption of isospin symmetry and recalling Eq. (11.27), the neutron structure 
function is simply obtained from Eq. (11.24) by the operation u ■v^ d; viz.. 



d{x) + d{x) + c{x) + c(x)j + \ \u{x) + u{x) + s{x) + s(x)j (11.44) 

J y L J 



Neglecting the c-quarks, the expression for the neutron-to-proton structure function ratio is 

F^i^x) u{x) -\- u{x) + A d{x) + d{x) + s{x) + s{x) 



(11.45) 



F^{x) 4 [u{x) + u{x)] + d{x) + d{x) + s{x) + s{x) ' 

From this equation it is readily apparent that the ratio is: 1/4 when u quarks dominate; 4 
when d quarks dominate; and, if the sea dominates, the ratio approaches unity. Within the 
parton model, Eq. (11.45) leads to the Nachtmann limit (Nachtmann, 1973), 

As interest has grown, a number of new methods for measuring the / ratio have been 
proposed. One approach (Fenker et al, 2003) is to tag very low momentum protons emerging 
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from the deuteron when deep inelastic scattering from the neutron is performed. This 
technique minimizes off-shell effects since at very low momentum the neutron is practically 
a spectator in the deuteron. Another method (Afnan et al, 2000, 2003; Bissey et al, 2001; 
Pctratos et al, 2006) is to perform deep inelastic scattering from "^He and "^H at very high x. 
In forming the ratio of the scattering rates from these two nuclei, the nuclear effects cancel 
to a high degree in extracting the F^/Ff ratio. 

Two proposed methods avoid nuclear effects altogether. Both a ratio of charge- current 
neutrino to antineutrino scattering from the proton and parity- violating deep inelastic elec- 
tron scattering at high x from the proton are sensitive to the d/u ratio. Nevertheless, both 
of these methods are fraught with technical difficulties, particularly at high x. 

The primary problem with measuring the neutron structure function is that there exists 
no practical free neutron target. Typically, a deuteron target is employed in experiments. 
The neutron structure function is then extracted from a measurement of the proton and 
deuteron structure functions by employing a model for the deuteron wave function. As 
an example, we consider the covariant approach of (Thomas and Melnitchouk, 1998). In 
this case, the proton and neutron structure functions are convoluted with a nucleon density 
function in the deuteron. The neutron is also off shell, so that a model for the off-shell 
behavior is necessary. Then the expression for the deuteron structure function is given by 



where PN/niv') is the probability of finding a nucleon of momentum y' in the deuteron 
and 5°^^ F2 is the off shell correction. The extraction process is iterative. First, the off- 
shell effect is subtracted from the measured deuteron structure function. Then the proton 
structure function, convoluted as indicated by the first term in Eq. (11.47), is subtracted from 
the remainder, leaving the convoluted neutron structure function. The neutron structure 
function is then deconvoluted to give the neutron structure function. The last three steps 
are then repeated until convergence is achieved. 

The results (Thomas and Melnitchouk, 1998) of this procedure for the ratio of the neutron 
to proton structure functions are shown in Fig. 11.10. In this work a covariant deuteron 
wave function and a consistent off-shell correction were used (Melnitchouk et al, 1994). 
Another application (Burov et al, 2004) of the covariant approach indicates that data for 
the deuteron structure function at very high x are essential for constraining the high x 




(11.47) 
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FIG. 11.10 The ratio of the neutron to proton structure functions as a function of x. The results 
of three extraction methods are shown. 

neutron structure function. Although the convolution approach has been used by a number 
of authors as a step toward explaining the EMC effect, the convolution method has no firm 
theoretical basis. 

Two other extractions, also shown in the figure, yield very different results at high x. The 
extraction, labeled "Whitlow et al (Paris)" (Whitlow et ai, 1992) made use of light- front 
kinematics with the null vector along the incident beam direction and an early deuteron 
wave function. The extraction labeled 'Whitlow et al (EMC)" (Whitlow et ai, 1992) made 
a density-dependent extrapolation of the EMC effect in heavy nuclei to that in the deuteron. 
There is significant controversy (Melnitchouk et al, 2000; Yang and Bodek, 1999, 2000) sur- 
rounding the density- dependent extrapolation method. 

A recent approach (Arrington et ai, 2009) made use of the light-front, with null-plane 
dynamics and modern deuteron wavefunctions. This framework, which can be applied even 
when the the constraints of DIS kinematics are not met, greatly simplifies the convolution 
equation in terms of the nucleon structure functions. The analysis established that it is 
important to account properly for the Q^-dependence of the data and incorrect to use a 
simple convolution formula to analyze data, if that data does not satisfy the DIS kinematic 
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constraints. When the CD-Bonn potential is used to determine the deuteron wave function, 
the method of (Arrington et ai, 2009) produces a ratio similar to that labeled "Whitlow et 
al (Paris)" in Fig. 11.10, which was obtained using a light-front approach with the null vector 
aligned along the direction of the incident electron beam (Frankfurt and Strikman, 1981). 
Whilst this particular choice of null-vector orientation is consistent with DIS kinematics, 
the very high x SLAC data, included in the analysis, are not. They are usually excluded 
when fitting PDFs. 



H. The longitudinal structure function 



Before the F2 structure function can be extracted from data, one must either know 
or measure the R structure function. By combining Eqs. (11.14), (11.15) and (11.16), the 
expression for the cross section in terms of -Fi(x,(5^) and F2{x,Q'^) is given by 



da 



2z?'2 



dndE 



4- 



F2{x,Q^) 2^ , 2Fi(x,g^) . 2O 

cos - H — sm - 

u 2 M 2 



(11.48) 



This cross section is often expressed in terms of the F2 and R structure functions in the 
following way: 



da 



dndE 



-{E,E',e) = A 



alE''F2ix,Q' 



2d 
cos - 

2 



1 



1 



e l + R{x,Q^) 



where 
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i + 2(i + i^Vg2)^^j^2^ 
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V J 



1 . 



(11.49) 



(11.50) 



(11.51) 



ar 2xFi{x,Q^) 

Here, ai and ut refer to the longitudinal and transverse cross sections, respectively. 
Since the quarks are pointlike spin-1/2 objects, the cross-section for absorption of a longi- 
tudinal photon is small in comparison with that for a transverse photon. Furthermore, a 
vector-vector interaction preserves chirality, so the longitudinal cross-section will depend on 
violations of this chirality. Violations of order m'^/Q^ would be expected, where rriq is the 
struck quark's current-mass. Intrinsic transverse quark momentum, kx, as well as higher 
orders in from the gluon distribution also can give rise to an increase in the longitudinal 
cross section. For example, the Bjorken-Feynman model gives the result (Feynman, 1972) 
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(11.52) 



The QCD contribution to F^, and consequently to i?, of order is given by (AltareUi, 
1982; Reya, 1981) 



where G is the gluon distribution function, and only four flavors of quarks are considered in 
the second term. An excellent extraction (Whitlow et ai, 1990) of F2 and R was performed 
from a global analysis of the SLAC data. 

It has long been believed that a tight constraint on the gluon structure function can 
be built from good measurements of Fl. The best extractions for R in the valence region 
are from the lower energy SLAC and JLab data (Tvaskis et ai, 2007; Whitlow et ai, 1990; 
Yang and Bodek, 2000). The best extractions of Fl are from HERA data, but at lower x 
than the subject of this review. Nevertheless, the deduced gluon distribution is very sensitive 
to the order of used in the analysis (Martin et ai, 2006). 

I. The gluon structure function 

Interest in the gluon structure function has grown markedly during the past decade. Part 
of the interest resides in the fact that gluons comprise more than 98% of the rest-mass of 
the nucleon as well as more than half of it's light-front momentum. This interesting aspect 
of the role of glue in the nucleon can be seen from Fig. 11.11 . Here the mass of a light-quark 
is plotted as a function of the quark's momentum. The "data" in the plot are results from 
lattice simulations, while the curves are from Dyson-Schwinger equation (DSE) calculations. 
The curve labeled "chiral limit" is the DSE result obtained with the current-quark mass set 
to zero. Clearly, even if the current-quark mass vanishes, the quark mass rises rapidly to 
near 0.3 GeV at low quark momentum. This is dynamical chiral symmetry breaking (DCSB) 
and a clear demonstration in QCD of the effect on the quark mass produced by the presence 
of gluons with strong self-interactions. 

Sensitivity to the gluon structure function via lepton beams can be produced by a careful 
measurement of scaling violations; i.e., measurements of the longitudinal structure function 
or the partial derivative of F2{x,Q'^) with respect to InQ^. This is particularly effective 
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FIG. 11.11 (Color online) Dressed-quark mass function, M{p): solid curves - Dyson-Schwinger 
equation (DSE) results, obtained as explained in (Bhagwat et al, 2003; Bhagwat and Tandy, 
2006), "data" - numerical simulations of unquenched lattice-QCD (Bowman et al, 2005). (NB. 
m = 70MeV is the uppermost curve and current-quark mass decreases from top to bottom.) One 
observes the current-quark of perturbative QCD evolving into a constituent-quark as its momen- 
tum becomes smaller. The constituent-quark mass arises from a cloud of low-momentum gluons 
attaching themselves to the current-quark. This is dynamical chiral symmetry breaking (DCSB): 
an essentially nonperturbative effect that generates a quark mass from nothing; namely, it occurs 
even in the chiral limit. [Adapted from (Bhagwat et al, 2007).] 

at low values of x where this derivative is directly proportional to the gluon distribution 
function at first order. 

III. DISTRIBUTION FUNCTIONS FROM DRELL-YAN INTERACTIONS 

Informative reviews on the Drell-Yan process exist (Kenyon, 1982; Reimer, 2007). In 
this section we emphasize that the Drell-Yan interaction provides: (i) the cleanest access 
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FIG. III. 12 Feynman diagram representing the Drell-Yan interaction 



to the antiquark distributions in the valence region in hadrons; (ii) a means to determine 
the quark distributions in the proton at very high x; and (iii) the most information on the 
parton distributions in mesons. Naturally, in order to isolate the valence quark distributions 
in the proton, it is imperative to measure the antiquark distributions. 

A. The Drell-Yan interaction 

The Drell-Yan process was devised (Drell and Yan, 1970a, 1971) to explain hadron-hadron 
collisions where an anti-lepton lepton pair is produced. For example, two hadrons, A and 
B, collide and produce the lepton pair 



At relatively large values of momentum transfer, say greater than a few GeV^, the underlying 
process is believed to be dominated by antiquark-quark annihilation: 



The Feynman diagram for this underlying Drell-Yan interaction is given in Fig. III. 12. 
Here, two hadrons collide and a time-like photon of mass Mi+i- is formed from the annihi- 
lation of a quark in one of the hadrons with an antiquark in the other. The photon then 
decays by emitting an antilepton- lepton pair. In a typical experiment, the lepton pair is 
often an antimuon-muon pair, as a matter of technical convenience. 

In this process the square of the four momentum transfer is just given by the square of 
the dilepton mass, M^+;_. At leading-order, the Drell-Yan cross-section can be determined 




(III.l) 




(III.2) 
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FIG. III. 13 Feynman graphs of NLO corrections to the Dreh-Yan process. Top row ~ Compton-Hke 
scattering in QCD; viz., gq — t- 7* (and also gq — ?• 7*, not shown exphcitly). Middle row - QCD 
radiative correction; viz., qq — )• 7*. Bottom diagram - QCD vertex correction. [Figure adapted 
from (Reimer, 2007).] 



in a straightforward manner to yield: 
da Anal 



dxhdxf 



(III.3) 



where qb{t)iixb{t)i) refers to the structure function of the quark of flavor i in the beam (b) 
or target (t) hadron, respectively. It has been shown (Altarelli et ai, 1979) that if the 
structure functions have the same definition as those in deep inelastic scattering; namely, 
that if one employs the DIS factorization scheme, then the NLO part is calculable and 
becomes a multiplicative factor to the expression above, as discussed in the next subsection. 



34 



B. QCD and higher order corrections 



The NLO QCD processes that contribute to Drell-Yan scattering are depicted in 
Fig. III. 13. These processes lead to a modification of the Drell-Yan cross section by in- 
troducing the so-called K factor: 

{NLO) = K^LOj^ (LO) . (III.4) 

With PDFs defined in the DIS factorization scheme, the i^TVLO-factor is given approximately 
by (Altarelli et al, 1979) 

Knlo ^ 1 + ^ (l + l^') (III-5) 
and assumes a value between 1.5 and 2. The consideration of NNLO, as well as NLO 
diagrams, also leads to a simple factorization of the cross-section and an approximate factor 
of two for K. The factorization scheme dependence of the K-factor is described at length 
in (van Neerven and Zijlstra, 1992). We note in addition that the K factor depends on 
kinematics, a fact shown in (Wijcsooriya et ai, 2005) to be important at very high x for 
pionic Drell-Yan studies. 



C. High- a; quark distribution functions 

The Drell-Yan process presents a valuable method for measuring parton distribution func- 
tions in hadrons at very high x. For example, it can be used to probe the quark distribution 
in the beam proton. To see how, consider that if s- and c-quarks are neglected and the 
beam-target kinematics are chosen such that xp '■= Xb — Xt is large, then for proton+proton 
collisions, Eq. (III. 4) can be rewritten (Webb, 2003) 

dfT A-ITCp' r — n 

-r-T ^ -TT- V ^^(^'" + ^^(^^' Q'^YAxt. Q') (ni.6) 

dxbdxt ols L -I 
because qb{xb) ^ <ib{.Xb) and qt{xt) <^ <it{xt) for large-Xi?. Now suppose that the target is a 
deutron with a similar kinematic setup, then Eq. (HI. 4) can be written (Webb, 2003) 

d(7 Attcx'^ K r 1 r — i 

— — — ^ ^ 4:Ub{xb,Q^) + dbixb,Q^) utixt,Q^) + dtixt,Q^) , (III.7) 

CLXbCLXi oiS '- -' '- 

where, as usual herein, q{x) means the distribution of fiavor-g- quarks in the proton; and 
one has assumed isospin symmetry and neglected nuclear binding effects in the weakly- 
bound deuteron. It is thus apparent that this kinematic setup produces Drell-Yan cross 
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FIG. III. 14 (Color online) Ratios of the measured Drell-Yan (squares) and pd (triangles) cross- 
sections to NLO calculations based on the MRST2001 PDF fit (Martin et al, 2002a), plotted as a 
function of beam momentum fraction, here labeled xi, with averaging over the target momentum 
fraction. The shaded area between the solid lines represents the experimental uncertainty ranges 
(Martin et al, 2003) on Au{x) + d{x) in the MRST2001 PDF fit. The data are from the FNAL 
E866 experiment. [Figure adapted from (Webb, 2003; Webb et al, 2003).] 

sections that are primarily sensitive to the valence distributions in the proton beam, and 
the antiquark distributions at small Xt in the proton and deuteron targets. 

Some results from analysis of the FNAL E866 Drell-Yan data are reproduced in Fig. III.14. 
In this figure the E866 results for pp and pd collisions are divided by the appropriate differ- 
ential cross-section computed using the MRST2001 PDFs (Martin et al, 2002a). Analogous 
ratios plotted as a function of Xt (Webb et al, 2003) indicate that the MRST2001 PDFs 
provide a very good description of the pp cross-section's Xj-dependence on the complete 
Xt-domain, and a good description of the pd cross-section for Xt < 0.15. In this light, con- 
sider first the pp data, which hints that the plotted ratio is smaller than one at Xf, ^ 0.6. 
Given that the w-quark is responsible for roughly 80% of the cross-section, this observation 
can be interpreted as an indication that the MRST2001 PDFs overestimate the proton's 
valence w-quark distribution. There is a strong signal from the pd cross-section that the 
plotted ratio is less than unity for ^ 0.5. Given that this cross-section is proportional 
to Au{x) + d{x) and the greater suppression, one can argue that the PDFs overestimate the 
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proton's valence (i-quark distribution. A consideration of the impact of this and other recent 
data on PDF fits is presented in (Owens et al, 2007). 



D. d/u ratio and the Gottfried sum rule 

One of the most celebrated applications of the Drell-Yan process is the measurement of 
the fiavor dependence of the antiquark sea in the proton. This was first suggested more than 
twenty years ago (Bickerstaff et ai, 1984, 1986; Ellis and Stirling, 1991; Garvey, 1987), at 
which time it was usually assumed that the light-quark sea was fiavor symmetric. Such 
experiments led subsequently to a demonstration that this is not true; namely, the light- 
quark sea is fiavor asymmetric. 

One of the first indications that the light-quark sea might be fiavor asymmetric was 
observation of the violation of the Gottfried sum rule. The Gottfried sum rule is defined by 
(Gottfried, 1967) 

Sg:= / -{F!{x)~F^{x)). (111.8) 
Jo X 

Using Eqs. (11.25) - (11.27) and assuming isospin invariance, Sq becomes: 

1 2 



1 2 /•! 
5g = ^ - o / dx 

6 6 Jo 



f dx [{d{x) - u{x)] =:-- -X,-_^ . (III.9) 



If the sea is fiavor symmetric, then the Gottfried sum just evaluates to 1/3. Any other value 
is termed a "violation" of the sum rule. 

The earliest hints of a possible violation of the sum rule can be found in data from SLAG 
(Bodck et ai, 1973) and Fermilab (Ito et ai, 1981). The NMG experiment at GERN gave 
the result (Amaudruz et ai, 1991; Arneodo et ai, 1994, 1997a) 

= 0.235 ± 0.026 , (III.IO) 

significantly different from 1/3. A more recent re-evaluation of the Gottfried sum, using 
a neural network parametrization of all then available data on the nonsinglet structure 
function F^^ = - F^ (Arneodo et al, 1997b; Benvenuti et al, 1989, 1990a), yielded 
Sg = 0.244 ± 0.045 X^_^ = 0.134 t 0.068 (Abbate and Forte, 2005), in agreement with 
the earlier analyses. 

Of course, it was pointed out long ago that Pauli blocking would give some enhancement 
of the ratio d/u (Field and Feynman, 1977). This mechanism might have been sufficient to 
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TABLE III.l Values of which appears in the Gottfried sum rule, Eq. (III.9), compiled in 

Peng (2003). The original sources are: NMC (DIS) (Amaudruz et aL, 1991; Arneodo et al, 1994); 
HERMES (semi-inclusive DIS) (Ackerstaff et al, 1998); and E866 (Drell-Yan) (Hawker et al, 1998; 
Peng et al, 1998; Towell et al, 2001). NB. At NNLO, is almost scale-independent for G 
[1, 100] GeV^, as illustrated in Fig. 3 of (Abbate and Forte, 2005). 



Experiment 


^d-u = lodx [{d{x) - u{x)] 


NMC 

HERMES 

E866 


0.147 ± 0.039 
0.16 ± 0.03 
0.118 ± 0.011 



explain the early SLAG data but Drell-Yan experiments have since given far more infor- 
mation about the magnitude of the effect. Pion cloud models also generate such an effect 
(Thomas, 1983). 

A compilation (Peng, 2003) of values for Xj_jj, which characterizes the second term in 
Eq. (111.9), are given in Table 111.1. The relative agreement is reasonable among these de- 
terminations. The most accurate result quoted is from the Drell-Yan experiment, FNAL 
E866, and the comparison is meaningful because Xj_^ is almost scale-independent for 
g2 e [1,100] GeV^. 

The power of the Drell-Yan technique can be illustrated by making some simple assump- 
tions. Consider the yield per nucleon from a proton beam incident on a target nucleus of 
atomic mass A, charge Z and neutron number N. Assuming that the yield, Yp^, from the 
process factorizes into a simple sum of pp and pn interactions, then 



YpA ~ (Zapp + Napn) /A 



(111.11) 



where cXpp and apn are the proton-proton and proton-neutron Drell-Yan cross-sections, re- 
spectively. If it is further assumed that pp Drell-Yan interactions are dominated by the 
target proton's -quark distribution and those for pn are dominated by the neutron's 
Un = c/p -quark distribution, then 

(111.12) 



pp 



Ur. 



cr. 



pn 



With Eqs. (111.11) and (111.12), it is readily seen that the ratio of yields for proton-nucleus 



38 



to that of proton-deuteron interactions becomes 

Y,A^, , {N-Z){u-d) 



(III.13) 



Ypd A{u + d) 

where we now return to our convention of omitting the subscript "p" when distribution 
functions in the proton are meant. Clearly, if the flavor of the light-quark sea is symmetric; 
i.e., u = d, then this ratio of yields becomes unity. By considering the ratio of pp to pd 
Drell-Yan scattering, then the ratio becomes: 

Ypd u + d 

and the ratio d/u becomes accessible experimentally (Bickerstaff et ai, 1984, 1986; Garvey, 
1987). Of course, the experiments are analyzed in a more sophisticated manner without 
these simplifying assumptions. 

The first Drell-Yan experiment that indicated that the sea was not flavor symmetric 
was experiment NA51 at CERN (Baldit et ai, 1994). This experiment found that u/d = 
0.51 ± 0.04 ± 0.05 at X = 0.18. Experiment E866 made a more comprehensive study by 
extending the x-range of the experiment (Hawker et ai, 1998; Towell et ai, 2001). These 
results are shown in Fig. III. 15. 

As argued, e.g., in (Steffens and Thomas, 1997), Pauli blocking is not sufficient to explain 
the large effect observed in modern experiments. A number of theoretical explanations have 
been advanced to explain this effect in the context of models: pion cloud; chiral quark; 
chiral soliton; and instanton. As real information about this phenomenon is only available 
for X ^ 0.3 and our focus is x ^ 0.4, we do not discuss it further herein but refer the 
interested reader to (Garvey and Peng, 2001; Kumano, 1998; Speth and Thomas, 1997), 
and references therein and thereto. 

From Eq. (III. 13), it should be apparent that the protonic Drell-Yan interaction with 
nuclei is an especially valuable means to measure anti-quark distributions in nuclei. In 
fact, this method has been used (McGaughey et ai, 1992) and proposed (Reimer, 2007) to 
search for an antiquark or, equivalently, a pion excess in nuclei. It is believed (Friman et ai, 
1983; Thomas, 1983) that observation of the pion excess in nuclei would provide a stringent 
confirmation of our understanding of conventional nuclear theory, where the nuclear binding 
is produced by pion exchange. Thus far, no evidence for a pion or antiquark excess in nuclei 
has been discovered (Gaskell et ai, 2001; McGaughey et ai, 1992). 
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FIG. III. 15 (Color online) x-dependence of dju. The solid blue squares represent the measurements 
from FNAL E866 (Hawker et al, 1998; Towcll et ai, 2001), while the red triangle indicates the 
CERN NA51 measurement (Baldit et al, 1994). The yellow shaded region indicates the uncertainty 
from the CTEQ6 evaluation of the world data. 

E. The pion structure function 

The pion plays a key role in nucleon and nuclear structure. It has not only been used 
to explain the long-range nucleon-nucleon interaction, forming a basic part of the Standard 
Model of Nuclear Physics (Pieper and Wiringa, 2001; Wiringa, 2006), but also, e.g., to 
explain the flavor asymmetry observed in the quark sea in the nucleon. However, compared 
to that of other hadrons, the pion mass is anomalously small. This owes to dynamical chiral 
symmetry breaking and any veracious description of the pion must properly account for 
its dual role as a quark- ant iquark bound-state and the Nambu-Goldstone boson associated 
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FIG. III. 16 The u quark distribution function in the vr". The data are from a LO analysis 
of the FNAL E615 experiment (Conway et al, 1989). The curve reproduces a Dyson-Schwinger 
equation calculation (Hecht et al, 2001), described in Sec. VLB. 3, which is consistent with QCD 
expectations. The discrepancy is disturbing and hitherto unresolved. 

with DCSB (Maris et al, 1998). It is this dichotomy and its consequences that makes an 
experimental and theoretical elucidation of pion properties so essential to understanding the 
strong interaction. 

Experimental knowledge of the parton structure of the pion arises primarily from pi- 
onic Drell-Yan scattering from nucleons in heavy nuclei (Badier et al., 1983; Betev et al., 
1985; Conway, 1987; Conway et al., 1989; Falciano et al., 1986; Guanziroli et al., 1988; 
Heinrich et al., 1991). A LO analysis of results from FNALE615 (Conway et al., 1989) 
are shown in Fig. III. 16 but the shape of the empirically extracted pion distribution function 
at high-x is contentious. 

Indeed, theoretical descriptions disagree. The QCD parton model (Ezawa, 1974; 
Farrar and Jackson, 1975), which determines the pion distribution function from the process 
depicted in Fig. III. 17, indicates that at very high x the distribution should behave as (1— x)^. 
Perturbative quantum chromodynamics (pQCD) (Brodsky et al., 1995; Ji et al., 2005) and 
continuum non-perturbative calculations, such as Dyson-Schwinger Equation (DSE) studies 
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FIG. III. 17 (Color online) Feynman diagram of the parton picture of DIS on the pion, where 
one hard gluon exchanges prior to absorption of the highly energetic photon, so that the struck 
parton carries most of the momentum of the pion. Evaluation of this diagram leads to a (1 — x)^ 
dependence of the pion structure function at very high x (Ezawa, 1974; Farrar and Jackson, 1975). 

(Bloch et al, 2000, 1999; Hccht et al, 1999, 2001; Maris and Roberts, 2003), which express 
the momentum-dependence of the dressed-quark mass function that is evident in Fig. 11.11, 
indicate that the high-x behavior should be (1 — xY^"', with an anomalous dimension 7 > 0. 

In contrast, AdS/QCD models using light-front holography (Brodsky and de Tcramond, 
2008) yield (1 — x)" with a = 0, as do Nambu-Jona-Lasinio models when a translation- 
ally invariant regularization is used (Bentz et ai, 1999; Davidson and Ruiz Arriola, 1995; 
Weigel et ai, 1999). On the other hand, NJL models yield a = 1 with a hard cutoff 
(Shigetani et ai, 1993), as do duality arguments (Melnitchouk, 2003). Relativistic con- 
stituent quark models (Frederico and Miller, 1994; Szczepaniak et ai, 1994) give (1 — x)° 
with < a < 2 depending on the form of model wave function; and instanton-based models 
produce (1 — x)"" with 1 < a < 2 (Dorokhov and Tomio, 2000). A full discussion is presented 
in Sec. VI and, in particular. Sec. VLB. 3. 

Given the importance of the shape of the pion distribution function at high a; as a test of 
QCD, a NLO analysis of the FNAL E615 data was performed (Wijesooriya et ai, 2005). The 
results of the this analysis are compared with those of the original LO analysis in Fig. III. 18. 
The solid curve is from the NLO analysis and in comparison with the LO analysis (dashed 
curve) has some strength shifted from the very high x region to the lower x region, as one 
should expect from the gluon radiation involved in the NLO processes depicted in Fig. III. 13. 
Nevertheless, the amount of additional depletion thus uncovered is not yet sufficient to agree 
with the non-perturbative DSE calculation, given by the black dotted curve in the figure, 
or the pQCD prediction, which both give a very-high- a; dependence of the form (1 — x)'^'^^, 
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FIG. III. 18 (Color online) The pion distribution function. Dashed curve - fit from the LO analysis 
of the FNAL E615 data; solid curve - NLO fit to the E615 data (Wijesooriya et al, 2005); and 
dotted curve - calculation of the distribution function using a Dyson-Schwinger equation approach 
(Hecht et al, 2001), which manifestly incorporates the momentum-dependence of the dressed-quark 
mass function that is depicted in Fig. 11.11. 

where 7 > 0. 

This discrepancy remains a crucial mystery for a QCD description of the lightest and 
subtlest hadron, and a number of explanations have been advanced to explain it. These 
range from simple experimental resolution-in-x problems (Wijesooriya et al, 2005); to a 
higher-twist effect that, at fixed Q'^{1 — x) with — 00, accentuates the structure function 
associated with longitudinal photon polarization (Berger and Brodsky, 1979); and theoret- 
ical factorization issues (Kopeliovich et al., 2005) in Drell-Yan at high x. Given this and 
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some questions regarding the Boer- Mulders effect, which are discussed in the next section, 
a new pionic Drell-Yan experiment with better resolution is warranted. 



F. Azimuthal asymmetries 

The decay angular distribution of the lepton pair in Drell-Yan interactions provides in- 
teresting additional insight into the valence structure of the hadron. In the simplest case, 
the decay angular distribution for a purely transversely polarized Drell-Yan photon is given 

by 



1\ ( da\ 
a) Vdn) 



(III.15) 



.167rJ 

where the angles are defined in Fig. 111.19. For the more general case where the Drell-Yan 
photon also has a longitudinal component, the decay angular distribution with angles defined 
in Fig. 111.19 can be written as (Collins and Soper, 1977) 



-) (-] 



3 1 



1 + A cos^ + fi sin 29 cos + sin^ 9 cos 2(f) 



(III.16) 



.47rA + 3j 

This expression is valid in all reference frames. Commonly used reference frames are the: 
u-channel, in which the 2;-axis is chosen antiparallel to the target beam direction; Gottfried- 
Jackson (Gottfried and Jackson, 1964) (t-channel) - z-axis is chosen parallel to the beam 
nucleon; and Collins-Soper (Collins and Soper, 1977) - 2-axis bisects the angle between the 
2;-axes in the other two frames. The quantities A, /i, and u in one frame can readily be 
related to their forms in another (Conway et ai, 1989). 

For pionic Drell-Yan interactions, a very interesting result is that the parameter u in 
the expression above was found (Badicr et ai, 1983; Bctcv et ai, 1985; Conway, 1987; 
Conway et ai, 1989; Falciano et ai, 1986; Guanziroli et ai, 1988; Heinrich et ai, 1991) to be 
large and dependent on the transverse momentum of the lepton pair, as shown in Fig. 111.20. 

In the late 1970's, two processes were investigated that could produce an azimuthal 
asymmetry in Drell-Yan processes. The first was a higher twist effect (Berger and Brodsky, 
1979), while the second was a single gluon radiation process (Collins, 1979). In the former, 
the high-twist diagrams considered gave rise to a cos 0-dependence in the Drell-Yan process. 
This is considered to be important at very high x and relatively low Q^. It is now believed 
(Boer, 1999) that high-twist effects, such as those explored in (Brandenburg et ai, 1994), 
cannot simultaneously describe the observed n and u in the Drell-Yan experiments. 
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FIG. III. 19 (Color online) The decay angular distribution of the Drell-Yan process in the Collins- 
Soper frame. 

In connection with the latter process, diagrams (c) and (d) in Fig. III. 13 were considered. 
The gluon radiation gives rise to the transverse quark momentum. This simple process 
produces a cos 2(f) dependence in the angular distribution of the Drell-Yan lepton pair. For 
example, it was found that (Collins, 1979) 



p2 



+ |p| ' 



^111.17) 



where Pt is the total transverse momentum observed in the process. It appears that this 
simple pQCD process could explain a significant fraction of the large v observed in the pionic 
Drell-Yan experiments, as shown by the dotted curve in Fig. III. 20. 

It is notable that for some time it was believed that if soft gluon resummation is included 
(Chiappetta and Le Bellac, 1986), then this process predicts very small values of v and 
restores in large part the naive Drell-Yan cross section given by Eq. (III. 15). However, it 
was recently shown (Boer and Vogelsang, 2006) that the gluon resummation in that work 
was not applied to both the numerator and denominator in the cos 20 component of the cross 
section. When the gluon resummation is applied correctly, then the effect partially cancels 
out, leaving the simple pQCD process as dominant. This result is verified in (Berger et ai, 
2007). It remains mysterious why the simple pQCD process overestimates u for the proton 
Drell-Yan experiment, depicted by the filled-squares in Fig. (III. 20). 

The mystery increases when the Lam- Tung relation is considered. In the context of the 
parton model, a relationship exists between two of the decay angular distribution parameters 
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FIG. III. 20 The v or cos 2(f) component of the Drell-Yan angular distributions obtained with a 
pion (NAIO and E615) or proton (E866) beam plotted as a function of the transverse momentum 
of the lepton pair. The dotted curve depicts a pQCD prediction (Collins, 1979). [Eq. (III. 17) and 
associated discussion.] The solid curve, which represents the Boer-Mulders effect (Boer, 1999), was 
adjusted to describe the NAlO data. [Eq. (III. 22) and associated discussion.] 

in Eq. III. 16 (Lam and Tung, 1980): 

l-\ = 2u. (III.18) 

The Lam- Tung relation is a consequence of the spin-1/2 nature of quarks. It is an analogue 
of the Callan-Gross relation in DIS, Eq. (11.22), but is less sensitive to QCD corrections. 
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Nevertheless, a violation of the Lam- Tung relation would suggest a rather significant non- 
perturbative process. The pionic Drell-Yan data indicate a large violation of the Lam- Tung 
relation, whereas the proton data do not indicate a violation. A more recent treatment 
(Bcrgcr et al, 2007) demonstrated that a resummed part of the helicity structure functions 
preserves the Lam- Tung relation as a function of Pt to all orders in a^- 

The observed violation of the Lam- Tung relation led to the suggestion of a new non- 
perturbative process (Brandenburg et ai, 1993). Therein, the violation of the Lam- Tung 
relation was parametrized in terms of k 

1-A-2z/ = -4k, (IIL19) 

where an Ansatz was proposed for the Pt dependence of n. After that work, a new structure 
function was advocated (Boer and Mulders, 1998), known now as the Boer- Mulders structure 
function. This quantity, h^^x, kr), represents the T-odd, chiral-odd structure function that 
describes quarks with Bjorken-x and intrinsic transverse momentum in one hadron, while 
h^{x,pT) represents the antiquarks in the second hadron. The cos 20 term arises from a 
double helicity flip process (Brandenburg et al, 1993), and is generally expressed by the 
product of two single chiral-odd, time-reveral odd (T-odd) helicity flip amplitudes, one for 
each of the hadrons involved in the process. Then the expression for u is proportional to 
the product of the two T-odd structure functions: 

z/ oc h^{x, kr) ® hi^x^pr). (111.20) 

[See also (Bodwin et al., 1989; Boer et al., 2003; Collins and Qiu, 2007).] It is notable that a 
similar product gives rise to the cos 20 asymmetry in semi-inclusive DIS, where the antiquark 
distribution is replaced by the Collins fragmentation function. A nonzero Boer-Mulders func- 
tion would signal a correlation between the transverse spin and the transverse momentum 
of quarks inside an unpolarized hadron. 

We note that in contrast to the parton distribution functions with whose properties we 
are primarily concerned, the Boer-Mulders and Collins functions are examples of "dynamic" 
structure functions, which do not have a probabilistic interpretation. As we have indicated 
here, if nonzero, these functions can lead to a wide range of effects that are not usually 
apparent in the parton model (Brodsky, 2009). 

An Ansatz similar to that in (Brandenburg et al., 1993) was used in later work (Boer, 
1999), wherein it was advocated that the parametrization of the Boer-Mulders structure 
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function should assume the same form as that for the Colhns fragmentation function (CoUins, 
1993). This parametrization uses a quark mass scale, Mc in a modified fermion propagator: 

hi{x,pT) = —CH^^^^exp{-aTPT)f{x) , (111.21) 

where Mh is the hadron mass, the parameter Ch is chosen to be unity, a-r = 1 GeV~^, 
and f{x) is a parton distribution function. In assuming that the Boer- Mulders structure 
function and the Collins fragmentation function have exactly the same form, the result for 
u reduces to 

The solid curve in Fig. 111.20 results from the above equation and gives a reasonable descrip- 
tion of the NAlO (pion) data. However, to be consistent with this data, the quark mass 
scale parameter must be unnaturally large [~ 2.3 GeV in (Boer, 1999)] in comparison with 
the chiral symmetry breaking scale, typically 0.3 GeV, as evident in Fig. 11.11. There is no 
sound basis for this and so, clearly, an improved theoretical description of this structure 
function is necessary. In order to accept this simple model as the explanation for the large 
asymmetry in pionic Drell-Yan data, one would also require an understanding of how it may 
correctly be combined with the apparently large pQCD component identified in (Collins, 
1979). 

The Boer-Mulders asymmetry was found (Zhu et ai, 2007) to be extremely small in 
the proton Drell-Yan data of FNALE866 depicted in Fig. 111.20. This small asymmetry is 
not understood, but it might be our first indication that both the Boer-Mulders structure 
function is small for sea quarks and the pQCD part is suppressed. The pionic Drell-Yan 
data involve valence quarks in both the pion and the nucleon target, while these proton 
Drell-Yan data were taken with kinematics that focused attention mainly on the valence 
domain in the beam proton and the sea region in the target. A very interesting test of this 
idea would be a Drell-Yan experiment with an anti-proton beam and a proton target where 
the kinematics were chosen such that the annihilating quarks could both be valence quarks. 
Experiments of this type are planned for the Facility for Antiproton and Ion Research (FAIR) 
at Darmstadt. New studies of pionic Drell-Yan will be initiated in the near future in the 
COMPASS experiment at CERN. 

A new suggestion (Lu et ai, 2006, 2007) is to make use of the cos 20 dependence in 
unpolarized Drell-Yan scattering to measure the flavor dependence of the Boer-Mulders 
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structure function for both quarks and antiquarks. In order to measure the quark (antiquark) 
distribution, a pion (proton) beam on both hydrogen and deuteron targets is proposed to 
perform the flavor separation. 

The Boer-Mulders asymmetry can also be measured in unpolarized semi-inclusive DIS 
(SIDIS). The general expression for this cross section (Ahmed and Gehrmann, 1999), written 
to emphasize the azimuthal dependence, is given by: 

7 2 

^ - (v4 + 5cos(/. + Ccos20), (111.23) 



dxdQ^dzdP^d(j) 2Q^z 

where z is the fraction of the energy transfer imparted to the produced hadron. In this case, 
the Boer- Mulders structure function also gives rise to an azimuthal asymmetry cos 20. In 
fact, the parameter C in Eq. (111.23) has a term that is proportional to the product of the 
Boer-Mulders structure function and the Collins fragmentation function. 

This asymmetry must, however, be disentangled from the Cahn effect (Cahn, 1978, 1989), 
which also gives rise to cos0 and cos 20 dependence. The Cahn effect can arise from the 
transverse momentum of the quark in a manner similar to that for R given by Eq. (11.52). 
At relatively large values of Pt, say > 1 GeV/c, the parameters B and C can also arise 
from gluon radiation effects (Gcorgi and Politzer, 1978). The azimuthal asymmetries arising 
from these effects have been very well established in SlDlS from: the CERN EMC experi- 
ment (Arneodo et ai, 1987); FNAL E665 (Adams et al, 1993); and HERA (Derrick et ai, 
1996). At Jefferson Lab the azimuthal asymmetries were found (Mkrtchyan et ai, 2008) 
to be small, and consistent with the Cahn effect and the Boer-Mulders effect at very low 
values of Pt- Future results from HERMES and COMPASS, and an upgraded JLab facility 
(Avakian et ai, 2006) will probe a more comprehensive kinematic space, so as to pin down 
these effects in the valence region and provide new information about the Boer-Mulders 
asymmetry. 

Given the importance of DCSB in QCD, it is also worth mentioning in connection with 
SlDlS that a leading-twist mechanism within QCD has been identified which can gener- 
ate a transverse spin asymmetry that directly probes partonic structure associated with 
chiral-symmetry breaking (Brodsky et ai, 2002a). This is the so-called "Sivers asymmetry" 
(Sivers, 1990, 1991). In hadron-induced hard processes; e.g., Drell-Yan, this asymmetry 
exists and is reversed in sign, and thereby violates naive universality of parton densities 
(Brodsky et ai, 2002b; Collins, 2002). Indeed, in the Drell-Yan process, even when both the 
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beam and target are unpolarized, the annihilating quark and antiquark have a transverse- 
momentum-dependent transversity. Extensive discussions of distributions associated with 
quark spin asymmetries are presented in (Barone et al, 2002; D'Alesio and Murgia, 2008). 

We close this subsection by reiterating that anti-proton-proton Drell-Yan interactions, of 
the type planned for FAIR, are a particularly powerful method for isolating the Boer- Mulders 
structure function. In this case, the cos 20 term is again given by Eq. (III. 20), where now 
the process involves annihilation of the quarks in the valence region of the proton with the 
valence antiquarks in the antiproton. Similarly, the Collins fragmentation function can be 
determined from e~^e~ — t- hh where h and h refer to the outgoing hadron and anti-hadron. 
Some work (Ogawa et al, 2007) has been performed at Belle in this regard. A recent review 
(D'Alesio and Murgia, 2008) of azimuthal asymmetries and single spin asymmetries in hard 
scattering processes covers these topics in more detail. 

G. The kaon distribution functions 

The valence structure of the kaon is comprised of a light u or d quark/antiquark and a 
strange quark/antiquark. If our understanding of meson structure is correct, then the large 
difference between the current- masses of the s-quark and the u- and d-quarks should give 
rise to some very interesting effects in the kaon structure function. For example, owing to 
its larger mass, the s-quark should carry more of the charged-kaon's momentum than the 
u quark. Then, the Uy quark distribution in the kaon should be weighted to lower values 
in X than that in the pion. A Nambu-Jona-Lasinio model calculation (Shigetani et al, 
1993) exhibits this behavior, as shown in Fig. III.21. A modern DSE prediction, discussed 
in Sec. VLB. 3 and depicted in Fig. VI. 42, confirms this trend and provides an improved 
understanding of its origin. 

Experimentally, a Drell-Yan measurement (Badier et al, 1980) of the ratio of K~ to 
7r~ is consistent with unity over most of the x region, with a suggestion that the ratio is 
dropping at high x. However, the data are not of sufficient quality to test and verify our 
understanding of pion and kaon structure. To explain, we note that one might claim from 
the data in Fig. III. 21 that for x > 0.7 there is a deviation from u^{x) = Uy{x) at the 3.8 a- 
level. However, following (Badier et al, 1980), which relied on simple models for motivation. 
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FIG. III. 21 Ratio of kaonic to pionic Drell-Yan cross sections, which is usually interpreted as the 
ratio of the u quark distribution function in the K~ to that in the tt~ . Data from (Badier et al, 
1980), which is obtained from a sample of dimuon events with invariant mass 4.1 < M < 8.5 GeV; 
and solid-curve from (Shigetani et al, 1993), an NJL-model result evolved from to Qq = 0.5 GeV 
to 4.5 GeV. (See also Fig. VI. 42 and associated discussion.) 



we observe that the entire data set is fitted by (see Fig. VI.42) 



R 



K-/n- 



1.1(1 -X 



,0.22 



(III.24) 



with = 0.45, and this drops to just = 0.032 for data at x > 0.7. Such values 

are unreasonably small. Hence, it is essential to make a high accuracy measurement of the 
structure function of the kaon at x > 0.6 before the behavior in Fig. III. 21 can be viewed as 
anything more than a hint. 
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FIG. IV. 22 Feynman diagram for a possible DIS interaction with a virtual pion in the proton. 
IV. FUTURE EXPERIMENTS AT JLAB, EIC, FNAL, J-PARC AND FAIR 

The 12 GeV upgrade at JLab permits a new opportunity to study structure functions at 
extremely high x. Although the greater energy that will be available at JLab will not be 
sufficient to avoid target-mass corrections for the nucleon at the higher values of x, it does 
present the prospect of obtaining precise data on the proton and neutron structure functions. 
It might also represent a new opportunity to study the pion structure function, although 
assumptions concerning virtual pions would have to be invoked (Wijesooriya et ai, 2001). 
This is because, in determining meson structure functions from deep inelastic scattering, it 
is necessary to consider scattering from virtual mesons, as illustrated in Fig. IV. 22. 

An electron ion collider (EIC) is aimed at polarization studies of the nucleon structure 
function at lower values of x. Nevertheless, it will permit access to the virtual pion and 
kaon structure function through measurements of the forward nucleon structure functions 
(Holt and Reimer, 2001). For example, the scattering process illustrated in Fig. IV.22 can be 
simulated using the RAPGAP Monte Carlo program (Jung, 1995). These processes include 
DIS from an exchanged pion or pomeron (D'Alesio and Pirner, 2000; Holtmann et al, 1996). 
A comparison of results from RAPGAP with HERA data show reasonable agreement for 
fast outgoing neutrons (Adloff et al, 1999). 

The expected accuracy of a possible DIS experiment from a virtual pion was estimated 
using events in which a "spectator" neutron was identified (Holt and Reimer, 2001). In 
general, the neutron is scattered less than 50 mrad from the nominal proton beam axis. 
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FIG. IV. 23 The expected statistical error for a future experiment to measure the vr" structure 
function at an electron-ion collider. [Figure adapted from (Holt and Reimer, 2001).] 

Events were cut on > 1 GeV^. For a 5 GeV electron beam on a 25 GeV proton beam, the 
expected errors are shown in Fig. IV. 23. A modest luminosity of 10^^ cm~^s^^ was assumed 
for a run lasting 10^ s. Such an experiment should be feasible because a similar experiment 
was conducted at HERA (Chckanov et ai, 2007). This experiment gave the very interesting 
result that the shape of the sea distribution for the virtual pion is the same as that for the 
proton. 

The structure function can be measured by considering deep inelastic scattering from 
the kaon cloud surrounding a proton. The basic Feynman diagram would be the same as 
in Fig. IV.22, with the pion replaced by a kaon and the neutron replaced by a A. The 
probability for scattering from the cloud surrounding the proton should be comparable 
to that for the 7r+ because the KNA coupling constant is comparable to that of the nNN 
vertex. In fact, one would only expect about a factor of two reduction in the vertex function 
for the kaon compared to the pion. 

The difficulty with this process is in the detection of the A. The A decays predominantly 
(64%) to a proton and a vr". Thus, a special forward proton spectrometer as well as a 
forward pion spectrometer would be necessary. This should be feasible since the ZEUS and 
HI experiments at HERA have already successfully employed forward proton spectrometers. 

It is also likely that an EIC will permit measurements of the F^/Ff ratio at extremely 
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FIG. IV. 24 (Color online) The absolute Drell-Yan cross section is proportional to 4ii(x) + d{x) in 
the beam proton. The data are from the FNAL E866 experiment. The solid circles are projected 
data from experiment E906 at FNAL. [Figure adapted from (Reimer et al, 2006; Webb, 2003; 
Webb et al, 2003).] 

high X. The technique should be analogous to the method that would be used for the pion 
and kaon structure function, with a deuteron target beam in this case. The scattered electron 
would be detected in coincidence with either a forward-going spectator proton or neutron, 
depending on whether the DIS occurred on a neutron or proton, respectively. In this case, a 
5 GeV electron beam on a 25 GeV deuteron beam at a relatively modest luminosity should 
be sufficient. 

Reliable results at extremely high x are also expected from forthcoming Drell-Yan mea- 
surements. For example, the solid circles in Fig. IV.24 indicate the error limits expected for 
the proposed Drell-Yan experiment at FNAL (E906) (Reimer et ai, 2006). The assumptions 
that went into this simulation were a 120 GeV proton beam from the FNAL main injector 
and a total integrated luminosity of 5 x lO'^^cm"^. 

Finally, the Japanese Proton Accelerator Research Center facility, J-PARC, currently 
under construction at Tokai, Japan, could make polarized Drell-Yan studies possible. Se- 
rious consideration is being given to experiments aimed at measuring parton distribution 
functions in mesons. In addition, the Facility for Antiproton and Ion Research, FAIR, under 
construction at GSI in Darmstadt would make it possible to directly measure the Drell-Yan 
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process from high-x antiquarks in the antiproton annihilating with quarks in the proton. 

V. PRESENT-DAY PARAMETRIZATIONS OF THE PDFS 

Excellent reviews of the parametrization of distribution functions exist. In fact, a good 
starting point for most of the distribution functions can be found at the Durham website 
(http://durpdg.dur.ac.uk/HEPDATA/PDF). At this website, one can access compilations of 
data, codes and grids associated with the following distribution functions: MRST/MSTW 
- Sec.V.A; CTEQ - Sec.V.B; GRV/GJR - Sec.V.C; and ALEKHIN - Sec. V.D, as well 
as an online PDF calculator. Here we consider primarily the CTEQ, MRST and the 
GRV parametrizations, emphasizing the differences between them. Another very useful 
website is that of CTEQ (http://www.phys.psu.edu/~cteq/). Here the CTEQ QCD hand- 
book (Brock et ai, 1995) can readily be accessed. The recent status of the MRST/MSTW 
PDFs can be found in (Thorne et ai, 2009), wherein the MSTW2008 distribution functions 
(Martin et ai, 2009) are recommended. 

The primary source of variations between the parametrizations are the different: 

• data sets used in the fits; 

• selections of data within the data sets; 

• pQCD choices - e.g., evolution order, factorization scheme, renormalization scale, as] 

• parametric forms for the PDF; 

• theoretical assumptions about the z — )■ 1 behavior. 

Although there are also distinct treatments of heavy flavors, assumptions of sea flavor asym- 
metry and X — )■ behavior, these do not have a large impact on the valence region. 

Most parametrizations begin with valence-like input, which means that at some infrared 
scale Qo ~ 1 GeV all distribution functions are represented as 

xfix,Ql)r^x''f{l-xff, (V.l) 

where «/ > and /?/ > are fit parameters, so that even those of the sea and glue distri- 
butions are nonzero but finite at the infrared boundary. Then, typically, the distribution 
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TABLE V.2 For illustration, parameters characterizing the MRST2001 NLO parton distribution 
functions at Q§ = 1 GeV^, defined by Eq. (V.2). 



parton 


Ao 


a 


13 


7 


6 


V 




0.040 


0.27 


3.88 


52.73 


0.5 


30.65 




0.158 


0.25 


3.33 


5.61 


0.5 


55.49 



functions are evolved to higher scales, commensurate with experiment, and the parameters 
fitted to obtain agreement with data. Some detailed starting distributions for the various 
parametrizations are given in the identified sections below. 

A. MRST and MSTW 

The data sets included in Martin-Roberts-Stirling-Thorne (MRST 2002) (Martin et ai, 
2002a, 2003) parametrization are: HI; Zeus, both neutral current and charge current data; 
BCDMS; SLAG; FNAL E665; CCFR; Drell-Yan data (FNAL E605, E772, E866); and FNAL 
CDF and DO data. The MRST parametrization has a starting scale of Qq = 1 GeV^ and 
accepts data down to = 2 GeV^. Thus, the SLAC data have an influence at large values 
of X. 

A parametrization of the parton distribution functions is given by 

xf{x, Ql) = Aox"{l - xf{l + Sx^ + r]x). (V.2) 

The parameters at the starting scale are given in Table V.2 for the MRST 2001 parametriza- 
tion, a NLO fit to the data. 

The chief difference between the MRST 2003 set (Martin et al, 2004c) and MRST2002 
is the exclusion of data below = 10 GeV^ and below x = 0.005; whereas, the MRST 2004 
analysis (Martin et ai, 2004b) included new HERA data at moderate values of x and high 
Q2. The MRST 2004 NNLO set was the first to use the full NNLO splitting functions. 
The newer MSTW parametrization (Martin et ai, 2009; Watt et ai, 2008) represents an 
update of MRST. This update has a number of new theoretical features aimed at the NNLO 
parametrization; e.g., NNLO corrections to the Drell-Yan data. This parametrization also 
includes NuTeV and Chorus data, the CDFII data, HERA inclusive jet data as well as 
direct high-x data on the structure function. The most interesting feature of these 
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FIG. V.25 Five parton distribution evaluations at NLO in the MS factorization scheme at = 
lOGeV^ in the valence region. There are three families of curves: the d{x)/u{x) ratios, the xu{x) 
distributions, and the xd{x) distributions. 

updates is apparent in the valence region of Fig. V.25. The NLO MSTW evaluation gives 
a much smaller xd{x) distribution at very high x than the previous work. Nevertheless, 
it is apparent in the differences between the parametrizations at large x that the d quark 
distribution is poorly constrained above x = 0.6. This further emphasizes the need for new 
data in the very high x region. 

B. CTEQ parameterization 

The most recent parametrizations from the Coordinated Theoretical-Experimental 
Project on QCD (CTEQ) is the CTEQ6 series (Lai et a/., 2007; Nadolsky et a/., 2008; 
Pumplin et al, 2007, 2001, 2002; Stump et al, 2003). The CTEQ collaboration omits data 
for < 4GeV^. In particular, CTEQ6 omits SLAC data as well as some high-Q^ HI 
data. CTEQ uses 10% systematic errors in quadrature with the statistical errors for the 
Drell-Yan data in comparison with 5% systematic errors assumed by MRST2002. CTEQ 
uses a starting scale of Ql = 1.69 GeV^. 
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TABLE V.3 Parameters at Qq — 1-69 GeV^ for the CTEQ6M parton distribution functions, defined 
by Eq. (V.3). 



parton 


Ao 


Ai 


A2 


^3 


^4 


A, 


dv{x, Qq) 


1.4473 


0.6160 


4.9670 


-0.8408 


0.4031 


3.0000 


Uy{X, Qq) 


1.7199 


0.5526 


2.9009 


-2.3502 


1.6123 


1.5917 


9{x,Qo) 


30.4571 


0.5100 


2.3823 


4.3945 


2.3550 


-3.0000 


{u + d){x,Ql) 


0.0616 


-0.2990 


7.7170 


-0.5283 


4.7539 


0.6137 


s(x,Qg) = s{x,Ql) 


0.0123 


-0.2990 


7.7170 


-0.5283 


4.7539 


0.6137 



The functional form of the CTEQ6 parametrization of the distribution functions is given 

by: 

xf{x,Ql) = Aox'^'il-x)^'e^'''{l + e^'x)^\ (V.3) 

The parameters at the starting scale are given in Table V.3 for the CTEQ6M parametriza- 
tion, a NLO fit to the data. 

C. GRV/GJR distribution functions 

The Gluck-Reya-Vogt (GRV) parton distribution functions were developed in a series of 
publications throughout the 1990s (Gliick et al, 1992a, 1993, 1995b, 1998). They are dy- 
namical distributions, which are generated radiatively from valence-like inputs at a low reso- 
lution scale. The latest of this series makes use of the 1994-95 HERA data for > 2 GeV^ 
as well as the SLAG, BGDMS, NMG and E665 data with > iGeV'^ and the simply 
extracted ratios FJ^/Ff from the NMG, BGDMS and E665 experiments. This analysis takes 
into account the Drell-Yan data and the u^/dy ratios extracted from the GERN GDHSW 
and WA21 neutrino data. 

The GRV parton distribution functions are parametrized as 

X/(X, Ql) = AoX"(l - Xf{l + 5y/x + 7]X). (V.4) 

In Table V.4 we report the parameters at the starting scale for the GRV98 LO parametriza- 
tion, a leading order fit to the data; and the parameters at the starting scale for the 
GRV98NLO parametrization, a next-to-leading-order fit to the data in the MS scheme. 
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TABLE V.4 Left panel - Parameters for the GRV98LO parton distribution functions at Qq = 
0.26 GeV^, defined by Eq. (V.4), where A = d — u and xs{x,Qq) = xs{x,Qq). Right panel - 
Parameters for the GRV98NLO parton distribution functions at Qq = 0.40 GeV^. 



parton 


Ao 


a 




<5 


V 


parton 




a 


/3 


5 


7? 


dy{x, Qq) 


0.761 


1.48 


3.62 


-1.8 


9.5 


d^{x, Qq) 


0.394 


1.43 


4.09 




18.2 


Uy(x, Qq) 


1.239 


0.48 


2.72 


-1.8 


9.5 


Uv{x, Qq) 


0.632 


0.43 


3.09 




18.2 


xA{x,Ql) 


0.23 


0.48 


11.3 


-12.0 


50.9 


xA{x,Ql) 


0.20 


0.43 


12.4 


-13.3 


60.0 


x{u + d){x,Ql) 


1.52 


0.15 


9.1 


-3.6 


7.8 


x{u + d){x,Ql) 


1.24 


0.20 


8.5 


-2.3 


5.7 


xg{x,Ql) 


17.47 


1.6 


3.8 






xg{x,Ql) 


20.80 


1.6 


4.1 







It is interesting that the value of (3 for m„ and systematically increases in going from 
LO to NLO. This is expected and consistent with the effect of gluon radiation reducing the 
hardness of the parton distribution. 

These distributions are characterized by a relatively low starting scale for evolution: 
LO, Qq = 0.5 GeV; and NLO, Qq = 0.63 GeV. One might question whether therefrom 
it is valid to employ pQCD evolution equations. In this connection, we observe that it 
is the combination a((5o)/[2vr] which appears in the evolution equations. At leading-order, 
a{QQ = 0.5GeV)/[27r] ^ 0.17 and hence (a(Qo)/[27r])=^ < 0.03. Therefore both the procedure 
and a physical interpretation of the distributions at this scale might be meaningful if higher- 
twist and essentially nonperturbative effects are not too important. At present one cannot 
answer whether the latter is true. Notwithstanding this, the procedure alone can be useful 
fitting and correlating tool. 

A more recent analysis (Gliick et ai, 2008) extends the GRV98 analysis. The distribution 
functions are generated from a radiated valence-like positive input with Qo < I GeV, 
where valence-like refers to Eq. (V.4). With these constraints, the predictions for Fl remain 
positive throughout the accessible kinematic region. Constraining the starting distribution 
to be valence-like leads to smaller quoted uncertainties in the very low-x region than those 
analyses without this constraint. 
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D. Distributions in relief 



A comparison of five modern parton distributions on tfie valence domain at = 10 GeV^ 
is shown in Fig. V.25. All five evaluations are in remarkably good agreement for the xu{x) 
distribution. Despite the large uncertainty in the d/u ratio discussed in sec. II. G, three of the 
distributions (MRST 2002NLO, CTEQ6.5M, GRV98) are in remarkably good agreement for 
the d/u ratio. The newest evaluations, MSTW2008 (Martin et ai, 2009), differ here largely 
because of a more flexible parametrization of the d quark distribution, and the inclusion of 
new data from FNAL on the lepton charge asymmetry from W decays and the Z rapidity 
distribution. 

The ALEKHIN2002 parametrization (Alekhin, 2003) is a somewhat difference case. This 
parametrization is based on the SLAG, BGDMS, NMG, and HERA data existing at the time. 
The cuts on the data included 2.5 GeV^ < < 300 GeV^ and x < 0.75. The model function 
was parametrized for Qq = 9 GeV^, much higher than most other forms. The cuts in 
and X give this parametrization more emphasis on the high-x SLAG data than GTEQ or 
MRST. This could possibly explain the differences in the high-x d-quark distribution. Later 
ALEKHIN parametrizations include Drell-Yan data (Alekhin et al, 2006) as well as even 
lower momentum transfer data (Alekhin et ai, 2008). However, the differences between 
these updated PDFs and the ALEKHIN 2002 distributions do not exceed one standard 
deviation. 

Notwithstanding differences between the Groups' PDFs, it is evident from the figure and 
the tables that all fits are consistent with 

44-^0. (V.5) 

However, this empirical statement has large uncertainties, some of which are discussed in 
(Owens et ai, 2007). It reflects primarily upon the fact that w-quarks dominate proton 
cross-sections and hence it is difficult to tightly constrain the proton's d-quark distribution. 

MRS (Sutton et ai, 1992a) and GRV (Gliick et ai, 1999, 1992b) also produced analy- 
ses of the parton distribution function in the pion. The pion distribution functions were 
extracted from Drell-Yan and prompt photon data. As we discussed in Sec. III.E, a more 
recent NLO re-analysis of the FNAL E615 pionic Drell-Yan data (Wijesooriya et ai, 2005) 
indicated that there was more curvature in the pion distribution at very large values of x 
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than given by either of these fits. As displayed in Fig. VI.41, this result strengthens the 
disagreement between the data and contact-interaction models, but does not produce as 
much curvature as expected from either pQCD (Ezawa, 1974; Farrar and Jackson, 1975) 
or nonperturbative DSE calculations (Hccht et ai, 2001). This discrepancy persists to this 
day and is one of the outstanding problems in understanding the application of QCD to 
real-world data. We therefore discuss it at length in Sec. VI. 

VI. THEORETICAL INTERPRETATION OF THE DISTRIBUTION FUNCTIONS 
A. Evolving insight from models 

We noted at the outset that while parton distribution functions are basic elements of 
factorized perturbative QCD, they are essentially nonperturbative; viz., they can only be 
calculated using a nonperturbative framework. Their calculation is a problem that has been 
with us for more than thirty years; e.g., (Politzer, 1974). The interim has seen the proposal 
and elucidation of numerous models for the phenomena of low-energy QCD, many of which 
have been employed to estimate the valence-quark distribution functions. 

1. MIT bag model 

Soon after it had been proposed, the MIT bag model (Chodos et ai, 1974a,b) was being 
used to calculate structure functions (Jaffe, 1975). It is still being used to make estimates, 
now of generalized parton distribution functions; e.g., (Pasquini et ai, 2005; Scopetta, 2005). 

The basic difficulties with this approach were noted at the outset. The bag is treated as a 
static cavity; namely, a spherical chamber with fixed radius, Rq. Amongst other things, this 
places a lower bound on the domain of x for which distribution functions can be calculated; 
i.e., X ^ l/[2MRo] ~ 0.1, where M is the nucleon's mass. Moreover, the bag boundary 
cannot realistically be static. It should respond to the action of currents on the quarks 
it contains and that would affect the form of the distribution functions. The target is 
treated as being at rest but the framework is not Poincare covariant and hence momentum 
is not conserved through the intermediate Compton scattering state. This entails that the 
calculated structure functions do not possess the correct support; e.g., they do not vanish for 
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X > 1. In addition, in an independent particle picture arising from a spin- flavor symmetric 
Hamiltonian the ratio in Eq. (11.45) is particularly simple; viz., 



because the pointwise behavior of the distributions is the same for both u and d quarks. 
This prediction conflicts with data. Fig. 11.10. Moreover, while the momentum sum rule 
is satisfied, that is achieved with all the target's momentum being carried by the quarks, 
whereas, experimentally, the quarks' momentum fraction should only be ~ 50%; i.e., (x) ~ \. 

From a modern perspective, the last of these outcomes is not seen as a critical impediment. 
The distribution functions are known to be scale dependent in QCD and, in principle, once 
determined at a particular scale, QCD predicts their pointwise form at any and all other 
scales. It should nonetheless be noted that, owing to the role played by gluons in binding 
the compound hadron, there is no scale at which the momentum fraction carried by valence- 
quarks can equal one. 

In practice, one views a model computation as providing the distribution functions at 
some infrared resolving scale, Qq. The evolution equations, discussed in Sec. II. D, are then 
used to obtain their form at some other scale, typically Q = 2 GeV. In applying a model, 
Qo is fixed a posteriori by requiring that the evolved distribution functions yield (x) ~ 0.4, 
which is an empirical feature of the PDFs described in Sec. V, as evident, e.g., in Fig. 2 of 
(Gliick et ai, 1992a). In contrast to modeling, however, in a truly nonperturbative solution 
of QCD the scale Qo would be known a priori. 

On the other hand, the bag model does yield some insight and positive results. For exam- 
ple, as treated, Bjorken scaling is recovered for Wi and {i'/M^)W2 in Eq. (11.13). Moreover, 
in the Bjorken limit, Eq. (II. 9), Wl vanishes and uWl scales and is calculable. Thus one 
preserves the Callan-Gross relation. The structure functions obey quantum number sum 
rules, such as the Adler relation: 



where F!^^'^^ are neutrino scattering analogues of the usual Ff'. In addition, Fi and F2/X 
are peaked at Xq = 1/iV, where is the number of quanta in the target: if three valence 
quarks are all the model possesses, then these functions peak at xq = 1/3. The width of the 
peak is related to the model's confinement radius. 



Vx : 



F^{x) ^ 2 
F^{x) 3 



{SU{&) — symmetric} , 



(VI. 1) 




(VI.2) 
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Primary amongst the problems identified above, in part because it also affects other 
models, is the absence of translational invariance and hence momentum conservation. In 
the usual mean-field approach the quarks move independently within a spherical cavity of 
fixed radius and definite location. The Hamiltonian associated with this problem is not 
translationally invariant. From this arises the so-called center-of-mass problem: particles 
moving independently within the volume cannot constitute an eigenstate of the total mo- 
mentum, in particular, they cannot properly describe a hadron at rest.^ 

In general, the removal of spurious center-of-mass motion from the hadron's wave function 
is not a well-defined problem and has no unique solution (Lipkin, 1958). The illness is typical 
of shell- like models. A number of prescriptions exist, which serve to mitigate the undesirable 
effects that this spurious motion has on physical observables. Amongst these, the Peierls- 
Yoccoz projection (Peierls and Yoccoz, 1957) has been most widely discussed. It projects 
exact eigenstates of the center-of-mass momentum by forming a linear superposition of at- 
rest bag states at different locations. Artifacts remain, however; e.g., some expectation values 
exhibit unexpected momentum dependence. A remedy may be found in the Peierls-Thouless 
projection (Peierls and Thouless, 1962) but that is difficult to implement. NB. The Peierls- 
Yoccoz projection is a nonrelativistic prescription and hence cannot properly account for 
the boost experienced by the intermediate state in the forward Compton scattering process. 
This means that for an initial state at rest, the procedure can only be internally consistent 
for \k\ < Mgg, in the notation of Fig. VI. 26, which corresponds to x ^ 0.6. 

In a little used alternative (Jasiak, 1997; Szymacha, 1984), one may augment the inde- 
pendent particle (bag model) Hamiltonian for the composite state by a fictitious attractive 
potential, which depends on the center-of-mass coordinate; viz., 

^^^^^""^^^^^^ ^^^^ + ^""^ ^ ^""^ + ^^^-^^ 
where i^BM is a bag-model Hamiltonian. With the potential tuned such that the spectra 
of Hqm and i^BM are approximately the same, then one has an equivalence between the 
nonlocal-composite and a translationally-invariant pointlike state. Unlike the Peierls- Yoccoz 
prescription, this method ensures that the valence-quark distribution functions are correctly 
normalized (Jasiak, 1997); namely, Eqs. (11.27) are satisfied independent of model details. 

A perhaps even greater difficulty is that in a interacting relativistic theory, which the model is aimed to 
be, Lorentz boosts change the particle number of a given state. 
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With a prescription for removing the spurious center-of-mass motion one can be confident 
that computed valence-quark distribution functions will possess the correct support; i.e., are 
nonzero only for < x < 1. Typical of such calculations are (Steffens and Thomas, 1994, 
1995), which employ the Peierls-Yoccoz projection following (Signal and Thomas, 1988) and 
achieve a fair degree of phenomenological success. 

In this connection it is notable that having overcome the problem of incorrect support, 
one can ask for more; e.g., models of the nucleon can be classified, and perhaps retained or 
discarded, according to whether or not they comply with the perturbative QCD constraint 
on the valence-quark distribution function (Avakian et ai, 2007; Brodsky et al, 1995) 



where n is the minimal number of fermion spectator lines in the quark-level scattering 
process, and 6X = \X'^ — is the difference in helicities between the struck parton and 
the hadron; e.g., SX is zero for a struck-quark with helicity parallel to that of the J = ^ 
hadron and one for a struck-quark with helicity antiparallel. These counting rules are first- 
principles predictions of QCD, in which the power increases logarithmically under evolution 
to an higher scale, and they satisfy "Gribov-Lipatov reciprocity" (Gribov and Lipatov, 1971, 
1972); namely, that at leading-order, spacelike and timelike parton cascades are identical so 
that structure and fragmentation functions are algebraically related. 
For the proton, Eq. (VI. 4) means 



for both u- and (i-quarks, unless nonperturbative effects somehow preclude either the valence 
M-quarks or rf-quark from sharing in the proton's helicity on a domain that includes the 
neighborhood of x = 1. NB. Subsequently, in connection with Eq. (VI. 98), we will describe 
a quantitative bound on the x-domain within which Eq. (VI. 5) should be observable. 

A naive consideration of the isospin and helicity structure of the proton's light-front 
quark wave function at x ~ 1 leads one to expect that (i-quarks are five-times less likely 
than ti-quarks to possess the same helicity as the proton they comprise. This leads obviously 
to the prediction (Farrar and Jackson, 1975) 




(VI.4) 




(VI.5) 




(VI.6) 
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q 




FIG. VI. 26 "Handbag diagram" contribution to the virtual photon-nucleon forward Compton 
scattering amplitude: ^*{q)N{p) — )• ^*{q)N{p). In this figure the intermediate state (or nucleon 
remnant in the DIS process) is assumed to be a scalar or axial-vector diquark, with masses Mg and 
M„, respectively. 

However, as we saw in Sec. V, whilst fitted PDFs produce /3„ ~ 3, they yield 4 < /3d < 5, 
the feature which leads to Eq. (V.5). This may be understood as an indication that the 
probability for a d-quark at x = 1 to possess the same helicity as the proton is actually 
much less than the naive expectation. Given the large uncertainties in the fits, however, it 
does not necessarily require that Pd > Pu- A reliable nonperturbative approach to nucleon 
structure is required in order to determine the relative strength of the /S^ = 3 and = 5 
components in the proton. 

As explained in (Brodsky et ai, 1995) and we subsequently elucidate,'^ Eq. (VI. 4) is ap- 
plicable at the infrared resolving scale, Qq. The exponent is a lower bound and increases 
under QCD evolution to a larger momentum scale. This is described in Sec. II. D. Equa- 
tion (VI. 5) can therefore serve as a stringent discriminator. It is satisfied with bag model 
valence-quark wave functions that fall as for large momentum \p\ (Margolis et ai, 

1993), whereas the oft used independent particle wave functions with Gaussian decay fail 
this test [see Eq. (VI. 11) and associated discussion]. 

The computation described in (Steffens and Thomas, 1995) is formulated with four pa- 
rameters: a bag radius R; the masses of the two valence- diquark spectator states in an 

^ A quantitative definition of Qq is presented in Sec. VLB. 3; and specifically in connection witii Eq. (VI. 83). 
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FIG. VI. 27 A bag model result for the x- weighted sum of valence-quark distributions in the proton. 
The model parameters are: bag radius R = 0.8 fm; scalar and axial-vector diquark masses, Mg = 
0.55 GeV, = 0.75 GeV; and infrared resolving scale Qo = ^ = 0.26 GeV. [Figure adapted from 
(Steffens and Thomas, 1995).] 

impulse approximation to the forward Compton scattering amplitude - scalar, M^, and 
axial- vector. My (see Fig. VI. 26); and the infrared resolving scale at which the computation 
is supposed to be valid, denoted by Qq above and fi therein.^ Naturally, the calculated 
valence-quark distribution functions depend on these parameters, which may therefore be 
fixed by applying leading- or next-to-leading-order evolution to the calculated distributions 
and requiring a good least-squares fit to the MRS parametrization (Martin et ai, 1993) of 
experimental data at = 10 GeV^. 

At leading-order. Fig. II.4, the values R = 0.8 fm, = 0.55 GeV, My = 0.75 GeV, 
Qo = fJ' = 0.26 GeV yield the results depicted in Fig. VI. 27. The diquark masses are 
roughly 20% lower than estimates derived from QCD's Dyson- Schwinger equations, Eq. (B.l) 



^ A comparison between bag model results and those of other, kindred non-topological solition models is 
presented in (Bate and Signal, 1992). Marked quantitative differences can be found in spin-dependent 
structure functions. 
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(Burden et al, 1997; Maris, 2002). Of more material concern is the small value of Qo, in 
which case the evolution equations involve the coupling a{Qo)/[2'K] = 0.4; cf. the values 
considered reasonable in Sec. V.C. 

Since this value of a{Qo) might be considered too large for the fitting- by-evolution pro- 
cedure to be valid, (Steffcns and Thomas, 1995) repeated the analysis using the next-to- 
leading-order evolution equations, Eqs. (11.30) - (11.35), and in this way obtained R = 0.8 fm, 
Ms = 0.7 GeV, = 0.9 GeV, Qq = fx = 0.34 GeV. While the value of Qo is still small, 
these diquark masses agree with contemporary estimates and now a{Qo) /[27r] = 0.1. In 
comparison with Fig. VI. 27, the NLO prescription cannot significantly change the evolved 
distribution because that is the quantity fitted. On the other hand, there is a signifi- 
cant difference in the valence-quark distribution at the new infrared scale Qq = 0.34 GeV: 
x[uv{x] Qo) + dy{x]Qo)] peaks at x = 0.4 (reduced by 30%), has a maximum value of 1.4 
(reduced by 20%) and material support only for x < 0.8 (domain contracted by 10%). These 
differences owe primarily to the accelerated production of sea-quarks through valence-quark 
depletion in NLO evolution as compared with the rate of this effect in LO evolution. The 
change in diquark masses serves to fine-tune the initial condition: more massive intermediate 
states produce a Qq distribution that peaks at smaller x, a feature required to fit the data 
parametrization when evolving from a larger value of Qo. 

As one might have expected, the computation reviewed here does not ensure correct 
normalization of the valence-quark distribution functions: roughly 20% of the strength is 
missing. The authors argue plausibly that this defect owes to modeling and parametrization 
of the intermediate state as opposed to its internally consistent calculation. They assert 
that the missing strength is probably localized in the sea-quark region; viz., x < 0.3, and 
one might therefore be justified in supposing that this weakness affects only modestly the 
pointwise behavior of the valence-quark distribution on the domain of interest. 

Upon the domain within which perturbative QCD is unarguably valid the parametriza- 
tions of valence-quark distribution functions are infinite ai x = while xg„(x; Qp^^j.^) = 0. 
The divergence is associated with parton splitting, see Fig. II. 4, and may be understood 
intuitively by observing that infinitely many x ~ partons are required to carry any fixed 
fraction of the proton's momentum. On the other hand, by assumption, model calculations, 
such as that just described, generate g„(x; Qq) < oo because only a fixed number of valence 
quarks (three for the proton) interact with the photon probe, and the valence-quark's struc- 
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ture is not resolved; i.e., there is no mechanism for them to fragment. While the property 
qv{x] Ql) < oo is not incompatible with the evolution equations, it does disagree with data 
parametrizations; e.g., see Sec. V. 

2. Constituent quark models 

In the minds of many, potential models framed in terms of constituent-quark degrees- 
of- freedom are exemplified by (Isgur and Karl, 1978). The application of such mod- 
els to the spectra, and strong and electromagnetic couplings of baryons is reviewed in 
(Capstick and Roberts, 2000). It is noteworthy that one can exactly eliminate the spu- 
rious center-of-mass motion in nonrelativistic constituent-quark models (CQMs) and hence 
they can be used effectively as a tool to study the spectrum. However, one should bear 
in mind that the absence of Poincare covariance is a problem when studying scattering, in 
general, and deep inelastic scattering, in particular. 

These models are defined by an Hamiltonian, which is typically diagonalized in an har- 
monic oscillator basis. The Hamiltonian contains a hyperfine interaction, whose structure 
is modeled according to a practitioner's taste. In many cases it is based on notions of one- 
gluon exchange [e.g., (Isgur and Karl, 1978)], in others, on pseudoscalar meson exchange 
[e.g., (Glozman et ai, 1998)]. Herein, it is only important that the Hamiltonian's eigenvec- 
tors do not exhibit SU{Q) spin- flavor symmetry and hence results other than Eq. (VI. 1) are 
possible. 

Relevant here is a relationship that appears empirically to hold between the anomalous 
magnetic moments of the neutron and proton, k„ p , and the valence-quark distribution 
functions [see Sec. 4.2.3 of (Goeke et ai, 2001)]; viz.. 



is the fraction of the proton's momentum carried by a valence-quark of type q (= u, d, etc.), 
at the scale Qq. Direct computation shows that even though the numerator and denominator 
on the right-hand-side (rhs) in Eq. (VI.7) are separately scale-dependent, the ratio is not. 
NB. We reiterate and emphasize that Eq. (VI.7) holds in a phenomenological sense: it is 




(VI.7) 



where, as usual. 
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satisfied to within 1% by modern parametrizations of parton distribution functions but has 
not been proven theoretically. 

The model described in Sec. VI.A.l yields = 2(x)rf^ and hence, unsurprisingly, 

from Eq. (VI.7). It is therefore evident that the magnitude of any deviation from 
this result can serve as one gauge of the amount by which SU (6) symmetry is broken within 
the nucleon. Experimentally, Hp/nn ~ —0.94, so the breaking is small by this measure. 

It has been argued (Traini et ai, 1997) that in CQMs, if one assumes 

1 



Mk)\' = m,n,{\k\) , (VI.9) 



where (pqik) is the momentum-space wave function for a constituent- quark of type q within 
a nucleon of mass M, and ^^(Qo) and ng{\k\) are, respectively, the constituent-quark's mass 
and momentum density distribution, then the virtual photon-nucleon Compton scattering 
amplitude, illustrated in Fig. VI. 26, can yield = ko + k^, see App. A) 

xq.{x;Qo) = '^^^^Jd'knM)Si^-^)- (VLIO) 

This expression can satisfy the momentum sum rule if sea-quarks and gluons are absent 
at Qo and ^^(Qo) = M/3. However, the distribution does not have correct support: it is 
nonzero for x > 1, and nor does it satisfy the normalization conditions, Eqs. (11.27). 
Both flaws are overcome if one adopts the prescription (Traini et ai, 1997) 

1 f k~^ ,:7*:s ^, X fc"*", 

1^ 1 



xq.ix; Qo) = ^^3^ J d'k - n,{\k\) - -) , (VI.ll) 

which exhibits two modifications. The simpler is the replacement mf ^ k~^, which provides 
the struck parton with an amount k~^ of the target's momentum, whereas it was at rest in 
Eq. (VI.IO). This ensures that Eqs. (11.27) are satisfied.^ More complicated is the modifica- 
tion of the Dirac (5-function's argument. It was inspired by an observation that only inelastic 
processes can contribute to the distribution functions. In this case the intermediate fermion 
line in Fig. VI. 26 should carry momentum k + g', where q' = (z/', q') with u ' = u — and 
z/ei = Q'^/[2M]. With this prescription the valence-quark distribution function is defined 
such that the energy transfer does not include the part absorbed in elastic scattering. An 
additional factor of 1/(1 — x)"^ is necessary to preserve the normalization. 



^ In the Bjorken limit, ~ xM in the target rest frame. 
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Equation (VI.ll) is a reasonable definition, which can be apphed to a large class of CQMs. 
That of (Isgur and Karl, 1978) is included. In this case the wave functions are Gaussian, 
from which it follows immediately that the model yields valence-quark distribution functions 
that violate Eq. (VI. 5); namely, ~ exp(— a;^/(l — x)^)/(l — x)^ for a; ~ 1. 

It is plainly the high-momentum components of constituent-quark relative-momentum 
wave functions that determine this behavior and Eq. (VI. 5) can be obtained from Eq. (VI. 11) 
if, and only if, 

nm ^'-^ (VI.12) 
\k\ 

(NB. The behavior l/|/c|^ gives (1 — x)^, while l/|/c|^ gives (1 — xY ) Hence, while 
the constituent-quark wave functions in the CQMs of (Bijker et ai, 1994; Ferraris et al, 
1995; Giannini et ai, 2001) provide greater support at high-momentum than those of 
(Isgur and Karl, 1978), in all cases it is still insufficient to give agreement with Eq. (VI. 5) 
via Eq. (VI.ll). We know of no CQM that can. 

Notwithstanding this, an illustration of the impact of SU{6) symmetry breaking via 
CQMs might be valuable. Such illustrations are provided in (Pasquini et al, 2002; 
Traini et ai, 1992). The more recent study compares distribution functions calculated in the 
models of (Glozman et ai, 1998; Isgur and Karl, 1978). Results are depicted in Fig. VI. 28. 
It is apparent that evolution does not dramatically affect the computed ratio; e.g., com- 
pare with the large effect evident in Fig. VI. 27. Hence a comparison between the models' 
predictions is straightforward. The Hamiltonian defining each model breaks SU (6) symme- 
try, albeit in slightly different ways, and therefore both produce curves which deviate from 
Eq. (VI. 1). Notably, the models exhibit quite different behavior at large-x. This is primarily 
because, in contrast to (Isgur and Karl, 1978), a relativistic kinetic energy operator is used 
in (Glozman et ai, 1998) and that increases the wave function's support at high-momentum. 

Along with the two CQM calculations. Fig. VI. 28 depicts the behavior of -F^/Ff that 
is inferred from the CTEQ5 parametrizations of experimental data (Lai et ai, 2000). In 
this connection it is instructive to recall Fig. 11.10 and the discussion in Sec. II. G, which 
emphasizes that F2/F2 is well known for x < 0.7. The parametrization ratio can therefore 
be viewed as a good representation of data on this domain. The absence of a free neutron 
target only introduces material uncertainty for a; ^ 0.7, a domain not reproduced in the 
figure. 
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FIG. VI.28 F^{x)/F^{x). Left panel - computed in the model of (Isgur and Karl, 1978): dot- 
dashed curve, at hadronic scale Qq ~ 0.4 GeV; dotted curve - leading-order evolution to Qq = 
lOGeV^, solid curve - next-to-leading-order evolution. Right panel - computed in the model 
developed from (Glozman et al, 1998). In both panels - Triangles: CTEQ5 fit to experimental 
data (Lai et al., 2000); and straight line: Eq. (VI. 1). NB. Using Eq. (11.45) and assuming St,(x) = 
on the valence-quark domain, the large-x behavior in these panels corresponds to dy{x) /uy{x) > 1.0. 
[Figure adapted from (Pasquini et al., 2002).] 

It is therefore significant that neither CQM yields behavior for F^/Fj' that is consis- 
tent with the parametrization ratio, not even for x ^ 0.4; i.e., the valence-quark domain, 
whereupon they are supposed by proponents to provide a veracious description of the nu- 
cleon via its salient degrees-of-freedom. For 0.4 < x < 0.7 the ratio in both panels lies well 
above the data parametrization. This is in contrast to models wherein the struck quark 
is partnered by scalar- and axial-vector-diquark intermediate states (Carlitz, 1975; Close, 
1973; Close and Thomas, 1988; Meyer and Mulders, 1991), to which we shall return. In 
(Pasquini et al., 2002) it is argued that the large-x behaviour of the the distribution func- 
tions is greatly influenced by proper implementation of the Pauli principle. That is neglected 
in (Carlitz, 1975; Close, 1973; Close and Thomas, 1988), which treat the diquarks as an ele- 
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mentary degree of freedom whose quark constituents are not active in the scattering process. 
Naturally, this is a weakness of those models, whose impact should be explored. (NB. That 
can be done within a Faddeev equation treatment of the nucleon; see, e.g.. Sec. VI. A. 7 and 



With the information described hitherto one can examine Eq. (VI. 7) within CQMs. To 
this end three different S'f/(6)-breaking Hamiltonians were considered in (Giannini et ai, 
2003); namely, those associated with the models of (Ferraris et ai, 1995; Giannini et al, 
2001; Isgur and Karl, 1978): Eq. (VI. 7) is satisfied to within 1% in each case. However, 
despite their differences, this collection of models belongs to a class'' whose characteristic is 
precisely. It seems therefore that Nature employs a different mechanism to break 
SU (6) symmetry than that expressed in this class of models. 

3. Pion cloud 

Chiral symmetry and the pattern by which it is broken are fundamental to the nature of 
hadron structure. Neither the MIT bag model nor constituent-quark models express QCD's 
chiral symmetry correctly. In working toward a theoretical understanding of distribution 
functions it is therefore important to explore the influence of correcting this defect. 

Pseudoscalar mesons come naturally to mind when considering chiral symmetry in QCD 
and mesons mean antiquarks. The addition of qq pairs or correlations to the intermediate 
state that is depicted in Fig. VI. 26 has long been considered a plausible mechanism by which 
to solve the normalization problem encountered in computing the valence-quark distribution 
functions (Schreiber et a/., 1991). 

This idea is canvassed in (Kulagin et al., 1996) within the framework of a dispersion 
representation of the valence-quark distribution function (Kulagin et al., 1994) 



Here p is the probability density associated with an intermediate spectator state of invariant 



^ It appears that this class includes the model of (Glozman et al., 1998) if a point-form spectator model 
(PFSM) is employed for the electromagnetic current operator (Meldc et al, 2007). Of the currents con- 
sidered in this connection, only the PFSM can produce results in fair agreement with experiment. 



App.B.) 




(VI.13) 
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ABC 

FIG. VI. 29 Depiction of three invariant-mass sectors conjectured in (Kulagin et al, 1996) to con- 
tribute to deep inelastic scattering. Left image - A: \ow-w - diquark spectators. Central image 
- B: intermediate- It; - compound spectator system involving a pion (dashed line), constituent- 
quarks and diquarks. Right image - C: large--u; - Regge trajectory model for qq scattering from a 
constituent-quark. [Figure adapted from (Kulagin et a/., 1996).] 

mass w = {P — ky and 

A;LxKx)=x(m2-^), (VI.14) 

where M is the nucleoli mass, is the maximum kinematically-allowed value of the struck- 
quark's squared-four-momentum, which is equivalent to its virtuality if the quark is massless. 
In connection with Eq. (VI. 13), convergence and internal consistency suggest that the prob- 
ability density can only have material support for /c^ ^ — M^. The condition 

kL^i^^x) = -M^ =^ = M^{l-x^)/x. (VI.15) 

Hence, in this case it follows that only intermediate states with w ^ are important 
for X > 0.6 . This explains the spectator diquark assumption of Fig. VI. 26. However, as x 
decreases below x = 0.6, intermediate states with s > will play an increasingly important 
role and must be included in order to obtain a pointwise-accurate result for qy{x; Qq). 

A concrete model is used in (Kulagin et al, 1996) to illustrate these points. It is notable 
that, by combining the components indicated below, this study aims to produce a valence- 
quark distribution function that is valid at a resolving scale Qq = 1 GeV; i.e., at a scale which 
may lie within the perturbative domain. We judge that perturbative QCD evolution from 
such a scale is sounder than from the manifestly infrared scales employed in Sees. VI. A. 1 
and VI.A.2. 
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The study divides the space of intermediate states into three sectors, which are illustrated 
in Fig. VI.29 and described here: 

A : w < M^. This is the dia gram of Fig. VI. 26, involving diquark spectators, which the 
studies described in Sees. VI. A. 1 and VI. A. 2 computed as the sole process contributing 
to the valence-quark distribution function. It is represented by a probabihty density 

pAiw,p^,x;Qo) = ^(«^ - "^z)) /Q/7v(P^a:; Qo) , (VI.16) 

D=0+,1+ 

with fq/j^jip'^, X] Qo) identified as the probability distribution for finding a constituent- 
quark in the nucleon with invariant mass p'^, proton light-cone momentum fraction 
X and partnered by a diquark with = D. Faddeev equation models of nucleon 
structure provide an appropriate framework for its evaluation. (See, e.g., App. B.) 

B : ^ w ^ Wq, Wq ^ 2 M^. One model for the addition of qq pairs or correlations to 
the intermediate state is to dress constituent-quarks with pseudoscalar mesons.^ In 
(Kulagin et ai, 1996) the dressing sum is truncated at just one pion rung, depicted 
in the central image of Fig. VI.29. The remainder is shifted to sector C. 

C : w > Wq. The large invariant mass spectator component is modeled via a Regge trajec- 
tory, with intercept ~ 1/2, to describe qq scattering from a constituent-quark. This 
piece provides a contribution to the valence-quark distribution function that behaves 
as x~°'^ for small-x. As we reported in Sec. VI. A. 1, this is a desirable feature. 

The model is elaborate and has numerous parameters. Their values are fixed through 
comparison with data and an eye to avoiding large discrepancies with legitimate theo- 
retical constraints. Important amongst the parameters are: the constituent- quark mass, 
rriQ = 0.45 GeV; the diquark masses, with the large values^ mo+ = 1.0 GeV and mi+ = 

® This can be sensible even when one unfolds and understands the structure of constituent-quarks 
(Blaschkc et ai, 1996; Cloct and Roberts, 2008; Ocrtcl et at, 2001). Indeed, it is an integral step in 
a systematic truncation of QCD's Dyson-Schwinger equations. 

^ These values are roughly 20% larger than extant computations of the mass-scales associated with diquark 
correlations, Eq. (B.l) (Burden et ai, 1997; Maris, 2002). 
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1.2 GeV appearing to be favored by data in this framework; momentum space widths used 
in parametrizations of constituent-quark-diquark Faddeev wave functions, Ao+ = 1.0 GeV 
and Ai+ = 1.2 GeV, and exponents for these vertex functions, no+ = 2.0 and ni+ = 3.5, 
chosen to guarantee Eq. (VI. 5); a hard- cut off- scale for the constituent-quark-pion vertex, 
Aqtt ~ 1.0 GeV; and a power-law form factor to model the p^- and fc^-dependence of 
the Regge exchange quark-quark scattering amplitude, which is defined by a mass-scale, 
Af( = 0.25 GeV, and a power, ri/j = 4. 

Combining these ingredients, the valence-quark distribution functions are computed at 
Qo = 1 GeV from the sum 

qv{x; Qo) = Zg{Qo){qA{x) + qsix) - mix) + qc{x) - qc{x)] , q = u,d, (VI.17) 

with the normalization constants fixed by Eqs. (11.27) subject to the additional constraints 
Jdxu_A{x) = 2, J dx d^lx) = 1. The relative normalizations of the contributions in 
Fig. VI.29 are determined once the model's parameters are fixed; i.e., they are not addi- 
tional independent quantities. Hence, in this model, defined at Qo = 1 GeV, the total 
valence-quark normalization is constituted from the diagrams in Fig. VI.29 as follows: 

dx {u,{x; Qo) + d,{x; Qo)} = A{56%) + B{16%) + C(28%) . (VI.18) 

Jo 

Moreover, 

dx X {u,ix;Qo) + d,{x;Qo)} = A{A1%) + B{6%) + C{Om%) = 0.47, (VI.19) 

Jo 

which leaves 53% of the proton's momentum to be carried by non- valence-quark degrees-of- 
freedom. 

Figure VI. 30 depicts features of the model's distribution functions. Although the sector- 
A contribution is always dominant on the valence-quark domain, it is evident from the left 
panel that the renormalization effected by adding the sector-i3 and -C contributions plays 
an important part within the model: it is a surrogate for evolution, shifting support to lower 
X at the model's mass-scale of Qo = 1 GeV. The right panel exhibits that behavior for the 
valence-quark flavor-ratio which is typically identified with diquark models: dy{x)/uv{x) de- 
creases monotonically with increasing x. The ratio plotted is formed only from the sector-^ 
contributions. On the valence-quark domain the full result may be estimated by multiplying 
the value plotted in Fig. VI. 30 by the ratio of normalization constants; namely, Z^jZ^ ^ 2/3. 
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FIG. VI. 30 Results illustrative of the pion-cloud model described in connection with Fig. VI. 29. 
Left panel: Total valence-quark distribution function. Complete calculation, combining all sectors 
[Solid curve]; A + B [dashed curve]; A alone [dash-dotted curve]. Right panel: Contribution from 
sector-^ to the ratio of valence-quark distribution functions. mQ+ = 1.0 GeV, mi+ = 1.2 GeV 
[solid curve]; and mQ+ = 0.6 GeV, mi+ = 0.8 GeV [dashed curve]. Both panels: Computed results 
are reported at Qq = lOGeV^, after LO evolution from Qo = 1 GeV, whereat a/[27r] ~ 0.03. 
Parametrizations of data (Lai et al, 1995; Martin et al., 1994) are provided for comparison [shaded 
band]. [Figure adapted from (Kulagin et al., 1996).] 

It is now clear that by implementing a scheme which resolves some of the structure of a 
constituent- quark it is possible to shift the model-scale Qo to a point wherefrom perturbative 
evolution may plausibly be justified a priori. In this way one can reproduce experimental 
data and provide an interpretation. However, the price appears high: with many ingredients 
needing to be carefully constrained and combined, the predictive power of an elaborate model 
is diminished. 

4. Topological soliton models 

Attempts to describe baryons as topological solitons are usually motivated by: empha- 
sizing the importance of chiral symmetry and the pattern by which it is broken in QCD; 
the notion that when a carefully-defined large-iVc limit is considered, QCD is a theory of 
weakly interacting mesons (and glueballs) and strongly interacting baryons; and the obser- 
vation that a simple form of classical nonlinear cr-model supports strongly-interacting static 
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topological soliton solutions, which can be quantized as fermions. Given this starting point, 
which possesses no explicit quark degrees-of-freedom, a connection with the parton model 
and DIS looks problematic. 

As in (Jaffc, 1975), a straightforward calculation of distribution functions for localized 
field configurations is possible when the current operator is at most quadratic in the fun- 
damental fields and the associated propagators are free-particle-like. The first of these 
conditions is difficult to satisfy in topological soliton models expressed in terms of el- 
ementary meson fields because the nucleon is a nonperturbative object involving all or- 
ders of the pion field and hence the current operator cannot be quadratic. This difficulty 
can be circumvented (Diakonov et al, 1996; Weigel et al, 1997) by realizing the nonlinear 
cr-model through a truncated bosonization of a four-fermion theory (Cahill and Roberts, 
1985; Ebert and Reinhardt, 1986; Roberts et al, 1988, 1989). Point coupling models of the 
Nambu-Jona-Lasinio type are often used. A regularization prescription and an associated 
ultraviolet mass-scale, Auv, are an essential part of the definition of this model. In general 
one cannot guarantee a priori that Auv plays no role in the Bjorken limit; namely, that 
Q^/Auv scaling violations are absent. 

It has been argued that within this framework the singlet unpolarized quark distribution 
function in a nucleon at rest may be written (Diakonov et al, 1996) 

-(^) + di^) = ^ / ^ E ^mil^ + 7V)]$n(p)| ^ , (VI.20) 

J y^'i ) n occup. " 

wherein: $„ is an eigenfunction, with energy En, of a Dirac Hamiltonian constructed using a 
so-called hedgehog soliton Ansatz for the pion field; the sum runs over the Dirac sea and the 
single discrete (valence) level bound in the soliton background, with —mq < < mq, con- 
taining Nc constituent- like quarks; and x G [—1, 1]. In defining Eq. (VI.20) one must follow 
the usual soliton model steps of projection and quantization to arrive at states with the cor- 
rect spin and isospin, and three-momentum. Moreover, it is necessary to employ a Poincare 
invariant Pauli-Villars regularization of the underlying Nambu-Jona-Lasinio model, or an 
equivalent, in order to preserve, amongst other things, completeness of the eigenfunctions 
and positivity of the rhs in Eq. (VI.20). NB. The singlet case is simplest in this approach 
because effects associated with cranking the topological soliton can be neglected. The dis- 
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FIG. VI. 31 Results obtained in a topological soliton model whose infrared resolving scale is ar- 
gued to be Qq ~ 0.6 GeV. Left panel: {x/2){uv{x;Qq) + dy{x;QQ)}, proportional to the valence- 
quark distribution: contribution from the single discrete (valence) energy level bound in the 
soliton background [Solid curve] and from the Dirac continuum [dotted curve]. Right panel: 
{x/2){u[x]Qq) +d(x;(5o)}) proportional to the antiquark distribution. Both panels: squares - 
NLO data parametrization drawn from (Gliick et ai, 1995a), associated with Qq = 0.56 GeV. 
[Figure adapted from (Diakonov et al, 1996): Auv = 0.56 GeV, the calculated nucleon mass is 
1170 MeV, and Nc = 3 with l/Nc corrections neglected.] 



tributions in Eq. (VI. 20) satisfy 



q{x) 



q{x) 



X > 



(VI.21) 



— g(— x), X < 
It can be shown algebraically within the model that 

J ^ dx [u{x; Qo) + d{x; Qo)] = dx [uy{x; Qo) + d^{x; Qo)] = Nc , (VL22) 

which may be compared with Eqs. (11.27). It is notable that the rhs of Eq. (VI. 22) is saturated 
by the discrete level alone and hence that the Dirac continuum, distorted by the presence of 
the large pion field, contributes nothing to the soliton's baryon number. In addition, with 
the Pauli-Villars regularization, 



dx x[u{x; Qo) + d{x; Qo)] 



dxxJ2{x] Qo) 



dxx[u{x; Qo) + u{x; Qo) + d{x; Qo) + d{x; Qo)] 



(VI.23) 
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so that at the resolving scale for which the model is assumed to be defined; viz., Qo ~ Auv = 
0.6 GeV, the proton's complete momentum is carried by the model's constituent-quarks and 
-antiquarks. 

Distributions typical of a topological soliton model are exemplified in Fig. VI. 31: left 
panel, valence-quark distribution; and right panel, the nonzero antiquark distribution. In 
this model the discrete level contributes to both the DIS valence-quark and antiquark dis- 
tributions. Furthermore, the Dirac continuum, distorted by the soliton, produces unequal 
distributions of quarks and antiquarks, and hence should not be solely identified with the 
sea. It will be observed that while the model's antiquark distribution is incorrect in detail, it 
does have the important property of possessing little support on the valence-quark domain. 
Moreover, owing to the nonzero antiquark distribution, which only vanishes as Nc — )■ oo, 
the valence-quark distribution peaks at x ^ 

A comparison between the left panels of Figs. VI. 30 and VI. 31 is instructive. Allow- 
ing for the difference in normalization, there is a semi-quantitative similarity between the 
soliton model result and that produced by the sum of sectors A and B in the model of 
(Kulagin et al, 1996); i.e., constituent-quark plus pion (see Fig. VI. 29). This highlights the 
omission of sea-quark effects in the soliton model at its resolving scale and also the role that 
sea-quarks play in general; e.g., with their inclusion, xq{x) has far greater support at low-x 
and commensurately less support on the valence-quark domain. 

It is interesting to compare the approach discussed here with that in Sec. VI. A. 1. In 
contrast to the static cavity treatment of bag models, wherein the three constituent- quarks 
contribute only to quark distributions, here they also generate an antiquark distribution 
owing to the strong pion mean-field, which distorts the negative-energy continuum of the 
Dirac Hamiltonian. In the context of Sec. VI. A. 1, this might be mimicked by incorporating 
effects of feedback between the constituent-quarks and the bag surface. 

In the approach outlined here the baryon number and momentum sum rules are automati- 
cally preserved. However, again owing to a lack of translational invariance in the formulation, 
the model's distribution functions do not vanish at x > 1 but behave as exp(— const. Ncx). 

That problem was considered in (Gambcrg et al, 1998), which argues that distribution 
functions in the infinite momentum frame (IMF) may be obtained from those calculated in 
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the nucleoli rest frame via the substitution:^'^ 

giMF(x) = ^(1 - a;) 'ZrfI- ln[l - x]) . (VI.24) 

The same formula was derived in formulating a translationally invariant bag model in 1 + 1 
dimensions (Jaffe, 1981), with the operation x — )■ — ln(l — x) identified as expressing Lorentz 
contraction of the bag-like object. Akin to Eq. (VI.ll), Eq. (VI.24) maps a function with 
support on [0, oo) to one with support restricted to [0, 1] and preserves normalization; i.e., 

/ dxqiMF{x)= dzq^piz). (VI. 25) 

Jo Jo 

For a typical distribution, it increases the peak height and suppresses the tail. The effects 
are large, perhaps too large for favorable comparison with experiment, and soliton rotation 
was not considered in the derivation (Wakamatsu and Kubota, 1998). While the alternative 
of merely neglecting the support problem is not meritorious, the prescription of Eq. (VI.24) 
is not widely employed. 



5. Pion structure function in the topological soliton model 

In topological soliton models chiral symmetry is emphasized and a strong pion mean- 
field is indispensable in order to realize a baryon. This said, it is natural to ask what 
the model predicts for the pion's valence-quark distribution functions. Suppose, as above, 
one specifies the theory via a bosonization of a point-coupling Nambu-Jona-Lasinio model, 
and in addition completes its definition by using a two-subtraction Pauli-Villars regular- 
ization scheme for the real part of the action. Under these conditions it has been argued 
(Gamberg and Weigel, 2000) that the pion's valence-quark distribution function is 



with 



V{mlx-Qo) = ^g2n(g^Qo) 



(VI.27) 



^'^ Recall that in the infinite momentuni frame the target's momentum is equivalent to its light-front mo- 
mentum and hence Bjorken-x coincides with the fraction of the target's light-front momentum carried by 
a parton. It is only in the IMF that g„ (x) can be interpreted as a single-parton probability density. 



80 



where {cj,Aiuv} are the Pauh-Villars parameters. It follows that in the chiral limit; viz., 
= 0, an analytic point in this calculation, ql{x]Qo) is a constant whose integral over 
X G [0, 1] is unity. Hence 

gr°''*°"(a^;Qo) = l, (VI.29) 

which, strikingly, means that at the model's resolving scale it is equally likely for the w-quark 
to carry any fraction of the pion's momentum: all, none, or some fraction in between, with 
the same probability. 

To place Eq. (VI.29) in context, the parton model prediction for the pion's valence-quark 
distribution function is (Berger and Brodsky, 1979; Brodsky and Lepage, 1979; Ezawa, 1974; 
Farrar and Jackson, 1975):^^ 

g:(a;;Qo) ~'(l-a;)l (VI.30) 

Plainly, Eq. (VI. 4) is also applicable to J = hadrons since n = 1 and 2 = 1 in this case 
(Brodsky, 2005). The extra multiplicative factor of (1 — x) is present because x = 1 corre- 
sponds to elastic scattering and hence must vanish for a spin- zero target. That imposes 
an helicity selection rule at the parton level, which causes (Tt to vanish a single power faster 
as X — 7- 1 than one would anticipate from dimensional counting (Llewellyn Smith, 1980). 
Owing to this, Eq. (VI.30) is consistent with the Drell-Yan-West analyses (Drell and Yan, 
1970b; Lepage and Brodsky, 1980; West, 1970). The soliton model prediction in Eq. (VI.29) 
conflicts conspicuously with the parton model result. 

6. Pion in the Nambu-Jona-Lasinio model 

Models of the Nambu-Jona-Lasinio (NJL) type have long been used in connection with 
the strong interaction. As we described in Sec. VI. A. 4, the NJL model provides a framework 
within which quark distribution functions may be defined in topological soliton models. 
However, it is more common for NJL studies to emphasize dressed-quark degrees-of-freedom 
and realize hadrons as poles in n-point Green functions. Following this path the model can 
be viewed from a Dyson-Schwinger equation (DSE) perspective (Gutierrez-Guerrero et ai, 

Equation (VI.30) generalizes to {x; Qo) ~ (1 — .t)^" for a pseudoscalar meson with n + 1 valence-quarks. 
Again, the exponent, 2n, is a lower bound, saturated at the infrared resolving scale. The value is increased 
by QCD evolution. 
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FIG. VI.32 Depiction of the general "handbag diagram" contribution to deep inelastic scattering 
from the pion. As we saw in Fig. VI. 26, with an explicit expression for the single-flavor removal 
amplitude, T-L{k,i), one defines a model framework within which to complete the calculation. 

2010; Roberts and Williams, 1994). 

The continuing appeal of the NJL model can be explained by two features; namely, it is 
defined by a simple interaction and exhibits dynamical chiral symmetry breaking (DCSB). 
This latter aspect, which the models of Sees. VI.A.1-VI.A.3 do not possess, suggests that it 
might sensibly be used to explore the structure function of the pion. As we noted above, 
the pion is special in QCD because it is both a dressed-quark-antiquark bound-state and 
the Goldstone mode arising from DCSB (Maris et al, 1998). 

By analyzing the Bjorken limit of the parton model expression for DIS from a pion 
target, one can derive a formula for the quark distribution function, which describes the 
quark-parton number density in the IMF: 



where k"^ = or, equivalently, k = {k'j_ + m\)/k^ . It follows from Eq. (VI.31) that 



independent of Qo? because the last matrix element is simply the contribution of a single 
quark flavor to the pion's electromagnetic form factor. 

Equation (VI.31) can be expressed in momentum space as follows: 




(VI.31) 




(VI.32) 




(VI.33) 
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Here the matrix trace applies to color and spinor indices, and 'HTr{k,i; Qq) is the amplitude 
defined by illustration in Fig. VI. 32: it is partially-amputated; i.e., the external pion legs are 
removed. In this form, 

J^' dx q:{x- Qo) = ^J -0^ trcD 7+^.(^, ^; Qo) (VI.35) 

and one recovers Eq. (VI. 32) because the integral on the rhs is merely an expression of the 
matrix element describing a single flavor's contribution to the form factor. As we may 
readily establish, this identity is preserved in the commonly-used rainbow-ladder truncation 
so long as in constructing T-Lnik, i; Qo) the quark-photon vertex satisfles the Ward-Takahashi 
identity and the pion Bethe-Salpeter amplitude is canonically normalized (Roberts, 1996). 

Within the NJL model, Eqs. (VI. 33) and (VI. 34) are not yet meaningful. A regularization 
scheme must be specifled and the result for gj(x; Qo) depends critically on the prescription 
adopted. To illustrate this we write explicitly the formal NJL model expression for the rhs 
of Eq. (VI.34): 

qvix;Qo) = ^ / ^^(^2n)^ ^'^CDi759nggS{i)-f^ S{i)i-f5g^qqS{k - ^)|^+^^^+ , (VI.36) 

where QT^qq is the Bethe-Salpeter amplitude for the pion^^ and 7"'" is the quark-photon vertex. 
Both are independent of the quark-antiquark relative momentum in a symmetry-preserving 
regularization of the NJL model owing to the momentum-independent interaction. 

One feature of a symmetry-preserving regularization is translational invariance, following 
from which one can use the Ward identity to rewrite 



Qvix; Qo) = \j ~(^r ^"^CD il5gngqS{e)-f+S{i)i-f5g^qqS{k - I) ^^^^^^ (VI.37) 



-Ii{k,x) 



8k^— ' ^^'-^'^ 

Il{k, X) = k+ J -^2^ n59nqqS{i)il59nqgS{k - ^)|^+=^fc+ • (VI.39) 



It is generally true that parton distribution functions are determined by a hadron's Bethe-Salpeter am- 
plitude. In the light-front quantization of QCD the PDFs can be computed from the absolute squares 
of light-front wave functions, integrated over the transverse momentum up to the resolution scale Qo 
(Brodsky et al, 1998). 



83 



Although written in a shghtly different form, Eqs. (VI. 38) and (VI. 39) are simply a restate- 
ment of Eq. (VI. 26). This is natural, given that the NJL model provides the framework 
within which parton distribution functions are defined in topological soliton models. 

At this point it is worth emphasizing that light-front concepts have hitherto appeared 
only as a result of a change in integration variables and it is generally assumed that the 
change of variables is nugatory. Hence one is not actually working within a light-front 
formulation of QCD. This makes no difference in perturbative analyses. However, it does 
involve significant, hidden assumptions in connection with important nonperturbative phe- 
nomena. For example, in order to evaluate Eq. (VI. 39) one requires the light-front form of the 
nonperturbatively-dressed quark propagator. That quantity is, however, unknown. Thus, in 
proceeding from Eq. (VI. 39), (Bentz et ai, 1999), for example, assumes that a dressed-quark 
mass, two orders-of-magnitude larger than the w-quark current-mass, is generated nonper- 
turbatively and that Qj^qg is related to the dressed-quark mass via a Goldberger-Treiman 
relation. Neither of these results has yet been proved in the light-front formulation. This 
problem was highlighted and considered in (Bentz et ai, 1999), see, e.g., Eq. (3.31) therein, 
but it was not solved, merely defined away. Thus, an internally consistent treatment of 
Eq. (VI. 39) is wanting, so that a covariant treatment following (Landshoff et ai, 1971) is 
preferable. 

Nonetheless, if one overlooks these caveats and proceeds as is customary, then it follows 
from Eq. (VI. 39) that in the chiral limit with a translationally invariant regularization: 

g:'''^-(x; go) = gr°'*°"(a:; Qo) = i, (vi.40) 

for all values of Qo, and consequently (Dorokhov and Tomio, 2000): 

^^n^.NJL,: f dx ql''''^ {x' Q,) = (VI.41) 

These simple results are themselves problematic and can directly be traced to the absence 
of a length-scale characterizing the pion's transverse size when a translationally-invariant 
regularization is employed for a chiral-limit NJL model. In this case both the interaction 
and the regularization procedure are momentum-independent, and hence neither can act 
to bound the relative momentum of the pion's constituents. This explains why the model 
disagrees so markedly with the QCD parton model prediction, Eq. (VI. 30): in QCD there 
cannot exist a resolving scale at which the relative momentum of the pion's constituents is 
unbounded. 
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FIG. VI. 33 Pion's valence-quark distribution function calculated in a Nambu-Jona-Lasinio model 
that is regularized by using a hard cutoff (Shigetani et al, 1993). Under one-loop evolution, a 
resolving scale of Qo = 0.5 GeV optimizes agreement with data parametrizations. 

One need not naively employ a translationally invariant regularization scheme, however. 
An ultraviolet cutoff, Auv, can serve to mimic asymptotic freedom in QCD, in the sense that 
it defines an NJL-model quark-antiquark interaction which vanishes for momenta larger than 
Auv- This procedure, adopted in (Shigetani et ai, 1993) for analysis of the virtual photon- 
pion forward Compton scattering amplitude, yields 

2 2 

g^^JLHc g^) ^ _g2_^ r^K e{Al^ - Ke) ^-^Q-^^^. e{mlx{l - x) - xM' - {1 - x)n) , 

JO [K fTT-Q ) 

(VI.42) 

where ke = — k,+ xm^ — xitlq /(I — x) and the constant of proportionality is fixed by requiring 
Jq dxq^{x;Qo) = 1. 

The valence-quark distribution function obtained from Eq. (VI.42), with Auv = 0.9 GeV 
and rriQ = 0.35 GeV, is depicted in Fig. VI. 33. It differs markedly from Eq. (VI.40). As 
will become clear, this could not have been otherwise because the behavior of q^{x) for 
a; ^ 1 depends on the precise form of the cutoff function. If this were not the case, then the 
valence-quark distribution function would be ineffectual as a discriminator between models 
for the strong interaction and, furthermore, no behavior could be uniquely identified with a 
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prediction of QCD. With the hard cutoff 



(VI. 43) 



From the distribution function displayed in Fig. VI.33 one obtains 

{^Vq^^'^'''' = 0-39 , = 0.22 , (VI.44) 

values which are significantly smaller than those obtained in the absence of a bound-state 
length-scale, Eq. (VI.41). It is apparent that in a model whose interaction acts at the 
resolving scale Qq to bound the relative momentum of the pion's constituents, the valence- 
quarks do not carry all the bound-state's momentum. The remainder may be associated 
with dressed-gluons, which effect binding within the pion bound state and are invisible to 
the electromagnetic probe (Hecht et ai, 2001). (See also Sec. V.C) 

The kaon's valence-quark distribution functions were also calculated in (Shigetani et ai, 
1993). With a current-quark-mass ratio of rris/ma ~ 25, they obtained ms/m£, ^ 1.5; 
and given this value for the ratio of constituent-quark-like masses, it is unsurprising that 
within the kaon the valence s-quark distribution is harder than that of the valence u-quark: 
the heavier valence s-quark must typically carry more of the bound-state's momentum. It 
follows from momentum conservation that u^{x)/u1{x) < 1 on the valence-quark domain. 
This outcome is consistent with experiment (Badier et ai, 1980), as illustrated in Figs. III. 21 
and VI.42. 

7. Nucleon in the Nambu-Jona-Lasinio model 

When emphasizing dressed-quark degrees of freedom, it is natural to describe baryons 
using a Poincare covariant Faddeev equation, aspects of which are described in App. B. 
Early attempts to compute the valence-quark distribution functions within this framework 
are described in (Kusaka et ai, 1997; Mineo et ai, 1999). These studies retained only a 
scalar diquark correlation and the latter employed a particularly severe truncation of the 
quark-exchange kernel in Fig. B. 43; viz., the so-called static-approximation, in which all 
momentum dependence is ignored in the propagator of the exchanged quark. This impacts 
significantly on the results. The studies also employ a simple current. In the former, that 
current is built only from Diagrams 1 and 2 in Fig. B.44. These diagrams are also the basis 
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for the current in (Mineo et ai, 1999) but therein, departing from the purely valence-quark 
setup, it is augmented by two diagrams argued to provide an estimate of the contribution 
from a pion cloud. 

Unsurprisingly, there are many similarities between the nucleon study in (Minco et ai, 
1999) and that of the pion in (Bcntz et ai, 1999). In particular, a translationally invariant 
regularization prescription is employed so that, in the absence of a pion cloud, the valence- 
quarks carry all the nucleon's momentum at the model's resolving scale, in analogy with 
Eq. (VI. 41). Some of this momentum is transferred to the pions following their addition to 
the current. The amount is related to the pion content of a dressed-quark, which is typically 
~ 10% (Cloet and Roberts, 2008). 

As found in (Kulagin et ai, 1996) and illustrated in Fig. VI. 29, (Mineo et ai, 1999) report 
that a pion cloud works to soften the valence-quark distribution and increase its support 
at low X. Unlike (Kulagin et ai, 1996), however, (Mineo et ai, 1999) still need a very 
low resolving scale; namely, Qo = mq = 0.4 GeV, if agreement with the parametrizations 
of (Martin et ai, 1994) is to be obtained under NLO evolution. The NJL model, with 
its momentum-independent interaction, must naturally produce hard distributions but the 
effect is compounded by the static-truncation used in (Minco et ai, 1999). As a result 
of this truncation, the nucleon's Faddeev amplitude is independent of the quark-diquark 
relative momentum, just as the pion's Bethe-Salpeter amplitude is independent of the quark- 
antiquark relative momentum in the NJL model formulated with a translationally invariant 
cutoff. 

As evident in Fig. VI. 34, the results are also sensitive to details of the regularization 
scheme. The favored scheme is that adopted in (Bentz et ai, 1999), wherewith a cutoff is 
placed on i±_ and qv{x;Qo) 7^ on < x < 1. However, in this case qv{x;Qo) exhibits a 
disconcerting, sharp and physically unacceptable increase at x ~ 0.9, which must artificially 
be suppressed before the evolution equations can be used. In practice, the distribution 
was forced to behave as (1 — x)^°. Owing to the behavior produced by both regularization 
schemes, this model cannot say anything meaningful about the distributions at Qo for x ^ 
0.8. This is a serious limitation because evolution is an area-preserving operation that shifts 
strength from large- to small-x. Hence, after evolution the behavior in this model of the 
distributions on the valence-quark domain is primarily determined by artifacts induced by 
the regularization. This is neither a sound basis for prediction nor for comparison with 
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FIG. VI. 34 Valence w-quark (left-panel) and d-quark (right-panel) distributions computed in a 
scalar-diquark picture of the nucleon based on a NJL model (Mineo et ai, 1999). The results 
exhibit a marked sensitivity to the manner in which the model is regularized. In this study the 
"TR" scheme is judged best. Both schemes lead to unphysical behavior on x > 0.8. NLO evolution 
of the model results is performed from Qq = 0.4 GeV. The parametrization is that of (Martin et al, 
1994). [Figure adapted from (Mineo et al, 1999).] 

parametrizations of data. 

Some general features of the distributions are nonetheless physically reasonable. For 
example, in this representation of the proton, the (i-quark appears only as a constituent 
of a [u, d] scalar diquark with mass uiq < mo+ < 2mQ. Hence, the probability of finding 
a valence d-quark in the proton is obtained by convoluting two probabilities; viz., that of 
striking a O"*" diquark in the nucleon with that of striking a (i-quark in the scalar diquark. 
Naturally, therefore, the valence ci-quark distribution peaks at smaller x than the valence 
M-quark distribution and is softer for 0.6 ^ x ^ 0.8. 

A realistic picture of the nucleon must include axial-vector diquark correlations, which 
generate significant attraction (Hecht et a/., 2002), and in (Mineo et al, 2002) they are 
added to the model just described. The Faddeev equation is still solved in the static- 
truncation but the regularization procedure denoted by "LB" in Fig. VI. 34 is adopted. As 
apparent, it forces the distribution functions to vanish at a; ~ 0.8, which, while still unphys- 
ical, is less difficult to overlook than the sharp rise produced by the "TR" scheme. The LB 
regularization procedure defines a model in which the distributions are extremely soft at Qo- 
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FIG. VI.35 F^{x)/F^{x) computed in (Mineo et al, 2002). r a measures the strength of axial- 
vector diquark correlations in the proton. It is a free parameter. Used in conjunction with mg = 
0.4 GeV and mQ+ = 0.68 GeV, the model-preferred value is = 0.25, at which the axial-vector 
correlation contributes just 7% to the proton's normalization, = reproduces the scalar- 
diquark only model, Fig. VI. 34. With = 2/3, the axial- vector diquark contributes 39% to the 
proton's normalization. NLO evolution of the model results is performed from Qq = 0.4 GeV. An 
explanation of the experimental points is provided in (Melnitchouk and Thomas, 1996). [Figure 
adapted from (Minco et al, 2002).] 

Needless to say, such behavior is inconsistent with Eq. (VI. 5), which specifies the behavior 
of distribution functions on a domain that includes the neighborhood of x = 1. 

The model has other internal difficulties. Notably, the axial-vector diquark is unbound 
and, owing to the absence of confinement, the A resonance cannot be formed without re- 
quiring the proton to possess what, within the model, is an unrealistically large axial- vector 
diquark component. Notwithstanding these problems, one might hope for useful informa- 
tion about the ratio -F2'^(x)/-Ff (x) on the valence-quark domain. Figure VI.35 displays this 
ratio's sensitivity to the strength of the proton's axial-vector diquark component. While 
the pointwise behavior in x is unlikely to be accurate, we judge that the general trend is a 
reliable indication of one effect of axial- vector diquark correlations within the nucleon; viz., 
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FIG. VI. 36 (Color online) Valence-quark distribution functions computed in (Cloet et al, 2005), 
which differs from (Mineo et al, 1999) primarily in the choice of regularization prescription. The 
marked sensitivity to regularization scheme in NJL-type models is evident through comparison 
with Fig. VI. 34. NLO evolution of the model results is performed from Qo = 0.4 GeV. "Empirical" 
denotes the parametrization of (Martin et al, 2002b). [Figure adapted from (Cloet et al, 2005).] 

a (i-quark in an axial- vector diquark can possess the same helicity as the proton target. 

The analysis of (Mineo et ai, 1999) is updated in (Cloet et al, 2005), wherein the model 
is redefined through a proper-time regularization procedure. In addition to the ultraviolet 
regulator, an infrared cutoff is also used in the reformulation. That expedient serves to 
mock-up confinement (Ebert et al, 1996) and thereby enables an improved description of 
the A resonance. 

In effecting a soft cutoff at Qq, the proper-time prescription yields much improved be- 
havior of the distribution functions on the valence-quark domain. This is illustrated in 
Fig. VI. 36, which was computed with: in GeV, Air = 0.28, Auv = 0.66, uiq = 0.4, 
mo+ = 0.65 mi+ = 1.2; and diquark coupling strengths = 0.5, = 0.08. NB. The 
splitting between the scalar and axial-vector diquarks is more than twice the usual size [cf. 
Eq. (B.l)] and the values of r^, suggest that the axial-vector diquark contributes < 2% to 
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the proton's normalization. 

Both the u- and d-quark distributions in Fig. VI. 36 behave as (Cloet, 2009) 

Qo) ~' (1 - x)^ e ~ 1.5 , (VI.45) 

a result which, however, still conflicts with the parton model constraint, Eq. (VI. 5): (1 —x)^. 
On the other hand, for the flavor ratio one flnds 

which corresponds to F^/Ff — )■ 0.43 as a; — )■ 1. [See Eq. (11.45) and cf. Fig. VI. 35.] How- 
ever, the agreement between Eq. (VI. 46) and the leading-order vector-gluon prediction of 
(Farrar and Jackson, 1975) is incidental. Instead, the value reflects the role and strength of 
axial-vector diquarks in the proton. 



8. Pion in the Instanton liquid 

The instanton liquid model derives from the notion that topological soliton solutions of 
the classical Euclidean-space Yang-Mills equations might play a role in describing QCD's 
ground state (Schafer and Shuryak, 1998). This possibility is being actively discussed 
(Ilgenfritz et ai, 2008; Lian and Thacker, 2007) and the model meanwhile continues to be 
used. 

In the present context, the instanton liquid model may be viewed as corresponding 
to a type of NJL model in which the interaction is momentum-dependent. With this 
feature it generates, e.g., a dressed-quark mass and pion Bethe-Salpeter amplitude that 
are momentum-dependent, and thereby overcomes a material weakness of the contact- 
interaction NJL model described above. On the other hand, the momentum-dependence 
differs signiflcantly from that of QCD: it is exponential, not power-law, and hence possesses 
too little support away from zero relative momentum, which is the domain most relevant to 
the valence-quark distribution. 

In this form the instanton liquid model has been employed to calculate the pion's valence- 
quark distribution function (Dorokhov and Tomio, 2000). As should be anticipated, the 
momentum dependence of the interaction leads to an x-dependent distribution at Qo = 
0.55 GeV, which vanishes at x = 0, 1, in contrast to that obtained with some regularizations 
of the contact-interaction NJL model; e.g., Eqs. (VI. 29) and (VI. 40). 
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Following leading-order evolution from Qo = 0.55 GeV to Q = 2 GeV, the model's re- 
sult for the distribution function on 0.2 ^ x ^ 0.9 is uniformly larger than experiment 
(see Fig. VI. 41), and the parametrization in (Sutton et ai, 1992b) which is based upon it.^^ 
The authors attribute this to their omission of gluon, sea and higher Fock state contribu- 
tions at Qo, which is plausible given the content of Fig. VI. 30. If the curve is multiplied by 
0.65, a factor representing a crude estimate of the renormalization attendant upon inclusion 
of the omitted effects, then the magnitude is in better agreement with the parametriza- 
tion. Following this rescaling, one finds that the large-x behavior is softer than that of the 
parametrization. Indeed, in this model 

q^^^ix; Qo) "~ (1 - xY, e ~ 1.53 - 1.56 . (VI.47) 

While this behavior also conflicts with the QCD parton model prediction, Eq. (VI. 30), it is 
a less striking disagreement than any other result described hitherto. 

B. Modern QCD perspective 

Today there are a few approaches to the calculation of parton distribution functions 
that are firmly founded in QCD. There is the numerical simulation of lattice- regularized 
QCD, of course. There are also frameworks which go beyond modeling by incorporating and 
expressing effects that are rigorously established features of QCD. In the best cases these 
approaches are directly connected with QCD via a symmetry-preserving truncation so that 
quantitative calculations can be performed with readily quantifiable errors, which can be 
estimated a priori. 

1. Generalized parton distributions 

These quantities (GPDs) were considered in (Dittes et al, 1988; Miiller et al, 1998) but 
have risen to prominence in hadron physics owing to their connection with deeply virtual 
Compton scattering (Ji, 1997; Radyushkin, 1996b) and hard meson production (Radyushkin, 
1996a). Unlike the usual parton distribution functions, GPDs are non-diagonal matrix 

Recall the critical reanalysis of the conclusions drawn from the E615 experiment, which is reviewed in 
Sec. III.E. 
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FIG. VI. 37 "Handbag diagram" contribution to the virtual Compton scattering process: 
7*(g)T(p) —7- 'y{q')T'{p'). Left panel - Kinematic variables in a symmetric reference frame. Right 
panel - Light-front representation of H'^ in Eq. (VI. 48). [Figure adapted from (Boffi and Pasquini, 
2007).] 

elements and hence do not have a simple probability interpretation. Instead, GPDs describe 
interference between amplitudes that represent different partonic configurations within a 
hadron and provide a means by which to chart correlations between different momentum 
states of a given parton. They are of interest here because the usual parton distribution 
functions can be expressed in terms of the GPDs.^"^ It follows that either: (1) valence-quark 
distribution functions can be used to constrain the parametrization of GPDs; or (2) complete 
knowledge of GPDs fixes the valence-quark distributions. The former link is used in present- 
day practice and that will not soon change. Herein we nevertheless elucidate connections 
between GPDs and PDFs. 

In Fig. VI. 37 we depict a handbag diagram contribution to a deeply virtual Compton 



scattering process, in which an high-energy photon with large spacelike virtuality, 



impacts upon an hadronic target, and produces a real photon and an hadronic final state with 
invariant mass <^ Q^. The amplitude depends on the Bjorken variable Xb = Q'^/{'2p ■ q), 
defined in Eq. (II. 5), and the momentum transfer t = [p' — pY = A^. It is assumed that 

For a spin-i hadron, h, with initial and final helicities A and A', respectively, one can 
define eight GPDs. For our purposes it is enough to consider only two of them explicitly; 
namely, the unpolarized distributions if^ and El: 

dz~ 



Q2- 



-p+ 



{h{p\\'m-^)^^qmh{pA)) 



2+=0,2x=0 



NB. For our purposes, we need not consider transversity distributions. 
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^ :u,{p', A') mix, e, th"- + Elix, e, t)7^a+"A A), (VL48) 



2P+ ' / L .^w , nv ... 

where x = k~^/P~^ is the average quark longitudinal-momentum-fraction and 
^ = [p'^ — p''^) I {j>^ + p'"*") is the so-called skewness. The GPDs also depend on the resolving 
scale Qo- 

These two GPDs are relevant herein because in the forward limit; i.e., p' = p, one has 

Hl{x = XB,^ = 0,t = 0;Qo) = { '^'^^ ^ , (VI.49) 

[ -qhi-XB] Qo), xb <0 

from which if follows, e.g., that 

rfxx^" Hl{x, 0, 0; Qo) = c/xx'" q'^{x; Qo) = {x^%, ■ (VI.50) 



Furthermore, 

dxHl{x,^,t-Qo) = Fl{t) 



1 



(VI.51) 



j\xEl{x,i,t- Qo) = FUt), 

where, e.g., F^f^{t) is the contribution to the hadron's Dirac form factor owing to quarks of 
flavor q. The standard normalization is such that 

Fi;(t = 0) = 2, (VI.52) 

F^it = 0) = 1 , (VI.53) 

F2pit = 0) = = 2Kp + K„ = 1.673 , (VI.54) 

F^j^{t = 0) = = + 2k„ = -2.033 . (VI.55) 

It is noteworthy that the right-hand-sides in Eqs. (VI.51) are independent of ^ because 
integrating over x eliminates all information about the particular light-front direction with 
respect to which ^ is defined; viz., owing to Lorentz invariance. 

Poincare invariance entails that there is only one generalized quark-parton distribution 
for the pion (Diehl, 2003): 



H^^{x,U;Qo) = I ^ et--^" (^-^(P')l^(-f)7+g(f)knp)) 



(VI.56) 

z+=0,zx=0 



Naturally, coupled with the fact that the pion's polarized GPD is identically zero; i.e., 
= 0,^^ it follows that Ag^(x) = 0; i.e., the pion's spin-dependent parton distribution 



15 



With the replacement 7+ — > 7^75 in Eq. (VI. 48), one arrives at the definition of Hh 
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function vanishes. When isospin is a good symmetry, one has 

Hl^^" := - Ht+ = -H^Z'' , (VI.58) 
= . (VI.59) 

The analogues of Eqs. (VI. 49) and (VI. 50) are correct, and, moreover, 

dx Hl{x, t- Qo) = F^{t) , F:{t) - F^{t) = 2F^{t) . (VI.60) 

2. Lattice QCD 

An operator expression for the pion's quark distribution function is given in Eq. (VI.31). 
We have seen in Sec. VLB. 1 that it is generally true; namely, in a spinless hadron "h" with 
total momentum /c, the distribution function is given by 

q\x-Qo) = — / dre-2^''^^'{hik)\,pMl^Mrmk)). (VI.ei) 

This expression is usually understood as representing the distribution in light-cone gauge; 
i.e., A'^ = 0, where is the gluon field. Equation (VI. 61) can readily be generalized to 
bound-states with spin. 

It is straightforward to modify Eq. (VI. 61) so that the expectation value is gauge invariant 
and yet unchanged in light-cone gauge; viz., one introduces a path-ordered exponential 
(Wilson line) 

Sir) :=Pexp|i^?^°d2;-yA+(0,z-,0^)| (VI.62) 

as follows 

q\x-S;Qo) := — / dCe'^^'^^^' {h{k)\ q{0)l^S{r)q{rmk)) . (VI.63) 

AtT J-oo 

With Eq. (VI.63) one has a gauge-invariant expectation value of a bilocal operator evaluated 
along a light-like line. It is mathematically precise and provides the pointwise behavior of 
the distribution function. However, it cannot be evaluated using the numerical approach of 
lattice-regularized QCD. 

Lattice methods can, however, be used to evaluate expectation values of local operators. 
Consider therefore 

:= jyxx''[q\x;£;Qo) - {-irq\x;£;Qo)] . (VIM) 
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It is notable that 

{x%, = £ dx q'iix; S; Qo) = < ; (VL65) 
namely, the number of valence-quarks of type q in hadron h. Furthermore, in general 

= £ dxx'-[q^{x;£;Qo)-q^ix;£;Qo)]= £ dxx'-q':{x;£;Qo) (VL66) 

(^'"^')qo = /' [^'(^; Qo) + ^'(^; ^; Qo)] ■ (vi.e?) 

It is apparent from these two results that the even moments express properties of a nonsinglet 
distribution whereas the odd moments probe the singlet distribution. Therefore they are 
independent and do not mix under evolution. Importantly, unlike the singlet distribution, 
a nonsinglet distribution does not receive contributions from disconnected diagrams. 

Introducing the Wilson line does not alter the distribution function's domain of support 
and, moreover, one still has q^{x) = —q^{—x). Hence, 

/I roo 
dxx''q\x;S;Qo)= dx x"" q\x;S;Qo) . (VI.68) 
-1 J —oo 

Now, upon insertion of Eq. (VI. 63), one obtains 

(^")qo = (A:+)-'^-^(Mfc)l V',(0)7+[^^+(0)]>,(0)|MA;)) , (VI.69) 

where 

D,{z) = ^-tg^A^,{z) (VI.70) 

is QCD's covariant derivative. (NB. A connection with Eq. (VI. 33) is readily established.) 
The moment on the left of Eq. (VI.69) is Poincare invariant but it only has a probabilistic 
interpretation in the infinite momentum frame. Owing to Poincare invariance, one may 
equally write 

{k,,k,,...k,J{x%^^ = {h{k)\^P,mr"^D^'{0)...^D^-{0)}^p,mh{k)), (VI.71) 

where the parenthesis indicate a symmetrization of indices and the subtraction of traces; 
e.g., {fc/io^Mi} ~ ^Mo^/ii ~ ^/ioMi^^- With Eq. (VI.71) one has an expression for the moments 
of quark distribution functions that can be evaluated using modern lattice-QCD methods 
(Soper, 1997). 

Computations are readily performed for low-order moments; i.e., n = 0, 1, 2, 3. Of course, 
owing to Eq. (VI. 65), there is no information in {x^): it merely provides a check on the 
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numerical procedure. One must also bear in mind that discretized spacetime does not 
possess the full rotational symmetries of the Euclidean continuum. Hence, for n > 4 it is 
impossible to define lattice-operators that are precisely equivalent to those of Eq. (VI.71): 
the lattice analogues cannot transform irreducibly and computations therefore require a 
numerical determination of the coefficients that describe mixing with operators of lower 
dimension. 

Suppose therefore that one focuses on the n < 3 moments; namely, those most directly 
accessible within the lattice framework. A typical lattice-QCD result for any one of these 
moments will contain, in addition to statistical errors, systematic errors from the following 
sources (Guagnelli et al, 2005): the so-called "quenched approximation" - which is not 
an approximation, but a truncation, that is still common today but will become less so 
in future; nonperturbative renormalization - the method used to obtain a renormalization 
group invariant moment from which continuum results can be inferred; discretization - 
extrapolation to zero lattice spacing; finite volume - extrapolation to infinite lattice size; and 
chiral extrapolation - inferring the value of a moment at the real-world current-quark masses 
from that computed at large, unphysical values, which is discussed, e.g., in (Detmold et ai, 
2003). The extent to which these errors have been accounted for and the nature of the 
procedure employed for each type of error must be considered when comparing the value 
of a moment obtained from a lattice simulation with its empirical or phenomenologically 
determined value. 

Lattice results are available for the first three moments of the pion distribution function; 
e.g., (Best et ai, 1997; Brommel et ai, 2007), with the results reported in Table VI. 5. The 
more recent study (Brommel et ai, 2007) employed 0(a)-improved Clover- Wilson fermions 
with two flavors of dynamical quarks and pion masses 430, 600, 800 MeV. A linear extrap- 
olation in was used to infer a result at the physical pion mass. The authors estimated 
rudimentarily that finite-size effects are of the order of 10% at the smallest pion mass used. 
Since disconnected diagrams were neglected in calculating all the moments and the simula- 
tion pion masses are large, it is reasonable to identify the results with moments of the pion's 
valence-quark distribution. This interpretation turns lattice artifacts to advantage in order 
to circumvent practically a problem of principle; viz., it is plain from Eqs. (VI. 66), (VI. 67) 
that odd moments of the valence-quark distribution are not readily accessible via the vector 
current operator. The valence-quark distribution is directly accessible though moments of 
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Qo (GeV) 


n = 1 


n = 2 


77, = 3 


Par.92 


2.0 


0.24 ±0.01 


0.10 ±0.01 


0.058 ± 0.004 


Emp.05 


5.2 


0.217 ±0.011 


0.087 ±0.005 


0.045 ± 0.003 


Lat. 07 


2.0 


0.271 ±0.010 


0.128 ±0.008 


0.074 ±0.010 


Lat. 


2.4 


0.273 ±0.012 


0.107 ±0.035 


0.048 ± 0.020 


Lat. 03Q 


2.4 


0.24 ±0.02 


0.09 ±0.03 


0.043 ±0.015 


DSEOl 


2.0 


0.24 


0.098 


0.049 




Qo (GeV) 


n = 1 


n = 2 


77, = 3 


Lat. 97Q 


2.4 


0.334 ± 0.021 


0.174 ±0.047 


0.066 ± 0.039 



TABLE VI. 5 Lowest nontrivial moments of the pion's valence u = d-quark distribution function. 
The rows report: "Par.92" - Moments at Qo = 2 GeV determined from the fits in (Sutton et ai, 
1992a); and "Emp.05" - Empirical result determined in (Wijesooriya et al, 2005) at Qo = 5.2 GeV. 
[NB. After evolution, these results are consistent with those in "Par.92".] "Lat. 07" - Lattice-QCD 
computation reported in (Brommcl et al, 2007) with nonperturbative evolution to Qq = 2.0 GeV 
(statistical and renormalisation error added in quadrature); "Lat. 97*^" - Quenched computation 
reported in (Best et al, 1997), with renormalization scale Qo = 2.4 GeV fixed via one-loop, chiral- 
limit lattice perturbation theory (statistical error only); and "Lat. 03" - Reanalysis of the results in 
(Best et al, 1997) using a different extrapolation to physical current-quark masses (Detmold et al, 
2003). The row labeled "DSEOl" is the result obtained using the Dyson- Schwinger equation 
framework in (Hecht et al, 2001). 

the more complicated V — A current operator. 

The earher study (Best et al, 1997) used Wilson fermions, in a quenched simulation, with 
current-quark masses that correspond to ~ 480, 650, 780 MeV. Values for the moments 
at a physical light-quark current-mass were obtained through linear extrapolation in m^. 
This could plausibly therein lead to overestimation of the first three moments by 10-20% 
(Detmold et al., 2003) (see Table VI. 5), and therefore also in (Brommel et al, 2007). 

It is curious and disturbing that there is little material difference between the results 
in Rows 3-5: all agree within quoted errors, even though they were obtained over a ten- 
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year time-span in which lattice methods improved greatly. In our view, for this reason and 
because the pion is so important to strong interaction physics, the evaluation of the pion's 
moments is a problem worth renewed attention. 

Low moments of the p-meson's valence-quark distribution functions are also inferred in 
(Best et al, 1997). There are three different structure functions for the spin-1 p: unpo- 
larized, with moments denoted by a„; polarized, with moments denoted by r„; and spin- 
projection, with moments denoted by dn- The spin-projection structure function is absent 
for a spin-i hadron. In the infinite momentum frame it measures the difference between 
quark distributions in = 1 and = targets. Herein, however, we are concerned with 
unpolarized distributions and so only ai.2,3 are listed in Table VI. 5. With the parameters of 
this lattice simulation, the vr and p moments differ noticeably and the first moments do not 
agree within statistical errors. It appears to us, therefore, that it is obtuse to claim that 
the results in (Best et al, 1997) support an assumption that Fi{x) ~ as was done 

therein. 

More attention has focused on the nucleon, with the first computations in 
(Martinelli and Sachrajda, 1989) and a brief contemporary overview in (Zanotti, 2008), but 
it remains true that disconnected diagrams are neglected in almost all extant studies. It fol- 
lows that only results for isovector quantities are directly comparable with those determined 
empirically; i.e.. 



For the reason noted above, calculations are restricted to at-most the first three nontrivial 
moments. We list some results in Table VI. 6. 

We present two sets of determinations from (Dolgov et al, 2002). Both were obtained 
with Wilson fermions. However, one set was computed in quenched-QCD; whilst the other 
was calculated using two degenerate flavors of dynamical sea quarks, with four current- 
quark mass values that provide for m^r/mp G [0.69,0.83]. Perturbative renormalization was 
employed to quote results at Qo = 2 GeV. The studies possess significant but unquantified 
finite-size and -volume errors, and a linear extrapolation in was used to infer results at 
physical current-quark masses. 

The fifth row of results in Table VI. 6 is compiled from (Gockeler et al, 2005), a quenched- 
QCD study that employed an 0(a)-improved Wilson action and fully nonperturbative renor- 
malization to reach Qq = 2.0 GeV. The authors argued that finite-volume errors were unim- 




(VI.72) 
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\^ /Qo 


n = 1 


n = 2 


n = 3 


Emp. 90 


0.16 


0.054 


0.023 


Emp. 02 


0.157 ±0.009 






Lat. 02Q 


0.251 ±0.018 


0.098 ± 0.068 


0.028 ± 0.049 


Lat. 02 


0.269 ± 0.023 


0.145 ± 0.069 


0.078 ± 0.041 


Lat. 04Q 


0.245 ± 0.009 


0.083 ±0.017 


0.059 ±0.018 


Lat. 08 


0.157 ±0.010 







TABLE VI. 6 Lowest nontrivial moments of the isovector nucleoli distribution function, all evolved 
to Qq = 2 GeV, albeit using different methods. Whilst they are all moments of nonsinglet distribu- 
tions, in principle only n = 2 is connected with a difference of valence-quark distributions. The rows 
report: "Emp. 90" - Empirical result determined from the parametrizations of rE(n^(x) ±(i^(x)) and 
xd^{x) in (Kwiecinski et al, 1990), which were fitted to then-extant data via next-to-leading-order 
evolution from a starting scale of Qo = 2.0 GeV [these values are little changed in more modern 
parametrizations, e.g., using MRST2002, available through the U.Durham data base (Durham, 
2009), {x)2Gcy = 0.157 ± 0.009, as indicated in the row labeled "Emp.02"]; "Lat.02Q" and 
"Lat. 02"- Lattice-QCD computations reported in (Dolgov et al, 2002), the former in quenched- 
QCD, both with perturbative evolution to Qq and linear chiral extrapolation (statistical error 
only); "Lat. 04*^" - Quenched-QCD study described in (Gockeler et al., 2005), with renormaliza- 
tion scale, Qq, fixed via nonperturbative renormalization and linear chiral extrapolation (error 
is statistical combined with estimates of systematic errors arising from renormalization, operator 
choice, and continuum and chiral extrapolations); "Lat. 08" - Analysis reported in (Hagler et al, 
2008), with renormalization scale, Qq, fixed via nonperturbatively-improved one-loop renormaliza- 
tion and extrapolation via covariant baryon chiral perturbation theory (Dorati et al, 2008) [before 
extrapolation: (x)"^'^ ~ 0.205 ± 0.015, statistical error only]. 

portant, and were attentive in estimating systematic errors arising from from nonperturba- 
tive renormalization, operator choice, and continuum and linear chiral extrapolations. 

We also list the estimate made in (Hagler et al, 2008) for the n = 1 moment, which 
was the only moment considered therein. This analysis employed a mixed-action approach; 
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FIG. VI. 38 (Color online) Current-quark-mass-dependence of , Qo = 2 GeV, computed di- 

rectly through numerical simulations of lattice-regularized QCD: Circles - (Ohta and Yamazaki, 
2008); Diamonds - (Orginos et al, 2006); Boxes - (Hagler et a/., 2008); Right-triangles - 
(Brommel et al., 2008); and Left-triangles - (Gockeler et al., 2005). The empirical value is marked 
by an asterisk. The Clover and Domain Wall Fermion results are in excellent agreement, except 
at the lightest pion mass, a discrepancy which can likely be attributed to finite size effects. The 
mixed-action results appear to differ from the others by a normalisation factor, which is perhaps 
understandable given that (Hagler et al., 2008) does not use the fully nonperturbative renormal- 
ization procedure employed in the other studies. [Figure adapted from (Zanotti, 2008).] 

namely, domain wall valence- quarks but 2 + 1 flavors of staggered sea-quarks, with light- 
quark masses corresponding to = 350 MeV. A nonperturbatively-improved perturbative 
renormalization scheme was used to quote a result at Qo = 2 GeV. It was argued that finite- 
volume errors are negligible with respect to statistical errors, whilst errors associated with 
nonzero lattice spacing (a = 0.124 fm) were not discussed. In order to quote a result at the 
physical pion mass, covariant baryon chiral perturbation theory was employed (Dorati et al., 
2008). As anticipated in (Detmold et al., 2001), an extrapolation nonlinear in can have 
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an enormous impact: in this instance, a 25% reduction in the value of the moment. 

A consideration of Fig. VI. 38 and the results listed in Table VI. 6 demonstrates that agree- 
ment between the empirical value of and that obtained using lattice-QCD is only pos- 
sible if there is significant nonlinearity in the m^-dependence of {x)^~'^ on < 0.1 GeV^. 
On the other hand, given that a chiral expansion of physical quantities is invalid for 
> 0.2 GeV^ (Chang et al, 2007) and reliable lattice results are restricted to this do- 
main, it is unsurprising that extant results from lattice-QCD do not exhibit nonlinearity. 

It is evident from the discourse in this subsection that contemporary lattice-QCD provides 
little aid in understanding the vast amount of DIS data, in general, and the valence-quark 
distributions, in particular. Practitioners express hope that this will change. However, it 
should be noted that low moments of parton distribution functions are only very weakly 
sensitive to the pointwise behavior of q'^ix^Qo) on x ^ 0.4. For example, consider two 
valence-quark distributions with large-x behavior (1 — x)^'^^ (Sutton et al, 1992a) and (1 — 
(^fjgcj-^^ fl/^ 2001). The exponents differ by 60% but the n = 3 moments, which 
are the highest accessible in modern lattice-QCD calculations, are indistinguishable within 
experimental error. One must compute the n = 5 moments before the difference exceeds 
20%. Plainly, pointwise calculations of the distribution functions are vital in order to probe 
the valence region. 

3. Dyson-Schwinger equations 

As elucidated in (Roberts and Williams, 1994), the Dyson-Schwinger equations (DSEs) 
are a system of coupled integral equations (DSEs) that have long been used in nuclear and 
particle physics. The DSEs provide a nonperturbative approach to QCD in the continuum 
and are particularly well suited to the study of QCD because of asymptotic freedom, which 
entails that mo del- dependence can always, in principle, be restricted to infrared momenta; 
viz., p ^ 1 GeV. Today, some elements of the approach are well-constrained even on that 
domain owing to positive feedback between DSE and lattice-QCD studies. Applications 
of DSEs to hadron physics are reviewed in (Alkofer and von Smekal, 2001; Fischer, 2006; 
Maris and Roberts, 2003; Roberts et al, 2007; Roberts and Schmidt, 2000). 

The physics of hadrons is ruled by two emergent phenomena; namely, confinement and 
dynamical chiral symmetry breaking. Confinement is the empirical fact that no quarks have 
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been detected in isolation. Dynamical chiral symmetry breaking (DCSB), which is respon- 
sible, amongst many other things, for the large mass splitting between parity partners in 
the spectrum of light-quark hadrons, which is present even though the relevant Lagrangian 
current-quark masses are small, explains the origin of constituent-quark masses and under- 
lies the success of chiral effective field theory. Neither confinement nor DCSB is apparent in 
QCD's Lagrangian but yet they play a dominant role in determining the observable charac- 
teristics of real-world QCD. 

With respect to confinement, it is important to appreciate that the static potential mea- 
sured in quenched lattice-QCD is not related in any simple way to the question of light-quark 
confinement. It is a basic feature of QCD that light-quark creation and annihilation effects 
are nonperturbative and thus it is impossible in principle to compute a potential between 
two light quarks (Bali et ai, 2005). 

On the other hand, confinement can be related to the analytic properties of QCD's 
Schwinger functions (Krein et ai, 1992; Roberts, 2008; Roberts and Williams, 1994), which 
are the basic elements of the DSE approach. Hence the question of light-quark confinement 
can be translated into the challenge of charting the infrared behavior of QCD's univer- 
sal /3-function.^^ Solving this well-posed problem is an elemental goal of modern hadron 
physics. It can be addressed in any framework enabling the nonperturbative evaluation of 
renormalization constants. 

Through the DSEs, the pointwise behaviour of the /9-function determines the pattern 
of chiral symmetry breaking. Moreover, the DSEs connect the /3-function to experimen- 
tal observables, so that a comparison between computations and observations of, e.g., the 
valence-quark distribution functions, can be used to learn about the evolution of the /3- 
function into the nonperturbative domain. 

A notable and relevant recent success of the DSE approach to hadron physics is its provi- 
sion of an understanding of DCSB via QCD's gap equation. The dressed-quark propagator 

This function may depend on the scheme chosen to renormalizc the theory but it is unique within a given 
scheme (Cclmastcr and Gonsalvcs, 1979). 
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can be written in a number of equivalent forms; e.g./^ 

1 



t-fpA{p\C) + B{p\C) 
^"'•«') (VI.73) 



i'-f ■ p + M{p^) 

It is important that the mass function, M{p^) = B{p^, C)/A{p^, C), illustrated in Fig. 11.11, 
is independent of the renormalisation point, Q. Furthermore, this propagator and, indeed, 
all colored Schwinger functions are gauge covariant. Therefore, when contracted in the 
expression for a color-singlet and gauge-invariant observable, all information about the evo- 
lution of the quark mass with momentum, displayed in Fig. 11.11, is retained. Hence, this 
manifestation of DCSB has an impact on cross-sections that is truly observable. 

That impact is greatly amplified in observables involving the pion, which, as QCD's 
Goldstone mode, is innately connected with DCSB. Indeed, the DSEs explain the pion as, 
simultaneously, both a Goldstone mode and an intricate bound-state of a dressed-quark and 
-antiquark, described by a Bethe-Salpeter amplitude: 

r,(A;;P) = 75 iE^{k- P) + ^ ■ PF^{k- P) + ^ ■ kk ■ P G^{k- P) + a^^k^P, H-,{k- P) . 

(VI.74) 

Through the axial-vector Ward-Takahashi identity, which is the statement of chiral sym- 
metry and the pattern by which it is broken, the DSEs yield (Maris et ai, 1998) a set of 
Goldberger-Treiman relations for the pion, which are exact in the chiral limit; viz., 

f^E^{k;0) = (VI.75) 

FR{k- 0) + 2 UF^{k- 0) = , (VI.76) 

Gn{k; 0) + 2 UG^ik; 0) = 2A'(P) , (VI.77) 

Hnik; 0) + 2 UH^{k; 0) = , (VI.78) 

where: the functions on the right-hand-sides are defined in Eq. (VI.73); here represents the 
pion's leptonic decay constant; and Fr, Gr, are associated with terms in the axial-vector 
vertex that are regular in the neighborhood of the pion pole. 



In order to maintain easy contact with existing literature, in this subsection we employ a Euclidean 
metric. In concrete terms that means: for Dirac matrices, {7/j,7iy} = 2(S^,y, 7^ = 'y^, = 74717273, 
tr757;^7,y7p7CT = —'^s^ivpa'i and a ■ b = X]i=i ^.A- A timelike vector, p^, has < 0. More information can 
be found in Sec. 2.1 of (Roberts and Schmidt, 2000). 
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The first of tfiese relations, Eq. (VI. 75), states that in the chiral limit the dominant, 
pseudoscalar part of the pion's Bethe-Salpeter amplitude is completely determined by the 
scalar piece of the dressed-quark self-energy. This exact result in QCD expresses an extraor- 
dinary fact; viz., owing to DCSB, the solution of the two-body problem in the pseudoscalar 
channel is known, almost completely, once the quark one-body problem is solved.^® In com- 
bination with the feature that the dressed-quark propagator has a spectral representation 
when considered as a function of current-quark mass (Langfeld et ai, 2003), Eq. (VI. 75) 
entails that the chiral-limit dressed-quark scalar self-energy provides an accurate pointwise 
approximation to E.,r{k; P) at physical light-quark current-masses. 

It will now be evident that, through the mesonic analogue of Fig. VI. 26, illustrated in 
Fig. VI. 39, the leading x-dependence of the physical pion's valence-quark distribution func- 
tion is determined by the chiral-limit dressed-quark scalar self-energy. As with the electro- 
magnetic pion form factor, the leading-order contribution is accurate apart from small cor- 
rections to the anomalous dimension (Maris and Roberts, 1998; Maris and Tandy, 2000a). 
This has far reaching consequences, which can be revealed by considering the gap equation. 

Let us imagine a theory in which the kernel of the gap equation behaves as 

where k > and M|, is some intrinsic, characteristic scale. Then, using the method of 
(Higashijima, 1984; Roberts and McKellar, 1990), one finds 



The case of QCD is similar: the kernel is 



B{e) (x-(-) . (VI.80) 



g2 g2 In[g2/A2^^p] ' 

where A = 12/[33 — 2nj], with n/ the number of active fermion flavours, and 



(VI.81) 



The following three relations establish that the pion necessarily has pseudovector components because the 
dressed-quark wave function renormalization is not identically one. This fact has an extremely important 
impact on the pion's electromagnetic form factor at large momentum transfer (Gutierrez-Guerrero et al., 
2010; Maris and Roberts, 1998). 
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The general pattern should now be clear; namely, the momentum-dependence of the chiral- 
limit dressed-quark self-energy replicates, up to an anomalous dimension, that of the quark- 
quark interaction and this behavior is manifest for k"^ ^ M|,. 

With this perspective, the analysis of (Ezawa, 1974) can now be understood to predict 
that, in a theory with an interaction of the type in Eq. (VI. 79), the pion's valence-quark 
distribution function must evolve according to 

q:ix; Qo) "S:' (1 - xf^^^''^ \/Qo > Mn ■ (VI.83) 

This indicates and interprets the scale at which "counting-rule" behavior should be evident; 
viz., it is the mass-scale at which the asymptotic form of the evolution of the chiral- limit 
dressed-quark scalar self-energy is manifest. 

One can readily translate the preceding argument into a QCD prediction; namely, 

g:(x;Qo)"5^'(l-x)'+^ Qo =^ 1 GeV, (VI.84) 

where < 7 ^ 1 is an anomalous dimension, because, as evident in Fig. 11.11, 1 GeV^ ~ 
20 Aqqj-) is the mass-scale at which the chiral-limit dressed-quark mass-function has plainly 
assumed its asymptotic form, Eq. (VI. 82). Given Eq. (VI. 83), it is unsurprising that the ex- 
ponent "2" in Eq. (VI.84) is a lower-bound. As we have indicated, the additional logarithmic 
correction to the momentum-dependence of the mass-function, a manifestation of asymp- 
totic freedom, leads truly to an exponent greater-than two at the counting-rule resolving 
scale. Moreover, this exponent increases under QCD evolution to Q > Qq. 

It is now instructive to reconsider Eq. (VI.83) in connection with models of the NJL-type, 
Sec. VI. A. 6. Such a model, in the chiral limit and regularized in a translationally invariant 
manner, may be realized through an analytic continuation of Eq. (VI. 79) to n = —1. One 
then reads the mass-function from Eq. (VI. 80): V/c^, B{k'^) = mq] i.e., a constant. The pion's 
valence-quark distribution function in this model follows from Eq. (VI.83): q'^^'^^'^^{x) = 1. 
In the absence of any regularization scale, this is the chiral-limit prediction at all resolving 
scales. Of course, in concrete applications of such models, an ultraviolet regularization scale 
is finally introduced: typically Ayy ~ 1 GeV. The model's practitioners describe this scale as 
the boundary between the nonperturbative and perturbative domains. As the only scale in 
the model's description of a chiral limit pion, it should properly define the model's intrinsic 
resolving scale. Hence, one arrives at the prediction 

g:^^^-U^;Qo=Auv) = l. (VI.85) 
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FIG. VI. 39 "Handbag" contributions to the virtual photon-pion forward Compton scattering am- 
plitude, vr, dashed-line; 7, wavy- line; S, internal solid-line, dressed-quark propagator, Eq. (VI. 73). 
The filled circles represent the pion's Bethe-Salpeter amplitude, r,r in Eq. (VI. 74), and the dressed- 
quark-photon vertex, discussed in App. B, depending on which external line they begin/end. 
[Figure adapted from (Hecht et al., 2001).] 

Our modernization of the analysis in (Ezawa, 1974) serves to eliminate much of the 
ambiguity in understanding what is meant by the "QCD prediction" for q'^{x). We have 
provided a clean definition of the resolving-scale at which this prediction should be valid. 
Moreover, with Eqs. (VI.83) - (VI. 85) and the associated discussion, we have demonstrated 
concretely that the pointwise behavior of the pion's quark distribution function on the 
valence-quark domain is a sensitive probe of the nature of the quark-quark interaction. 
To be concrete, if the conclusion of (Conway et al, 1989) is confirmed, then the theory 
underlying hadron physics is not QCD. The arguments can be generalized to other hadrons. 

Additional information can be explicated through a recapitulation of a calculation of 
the pion's valence-quark distribution function within the DSE framework. The diagrams in 
Fig. VI. 39, an explicit realization of Fig. VI. 32, provide the starting point for the computation 
in (Hecht et al, 2001). The upper diagram represents the renormalized matrix element 



T+(g,P) = tr 



/p^|r-r.(.,-P, 
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X S{kt) leQT^ih, k) S{k) ieQT^{k, h) S{kt) T+T^{kr; P) S{h) , (VI.86) 



where: r± = |(ri ±2x2) is a combination of Pauli matrices; S{i) = dia.g[S u{i), Sii{i)], 
with Su = Sd = S, assuming isospin symmetry; and kr = k — q — P/2, kt = k — q, 
ks = k — q — P. The fully-dressed quark-photon vertex, Tf^{ii,i2), appeared in Eq. (VI.34) 
and here Q = diag(2/3, —1/3) is the quark-charge matrix. The matrix element represented 
by the lower diagram is the crossing partner of Eq. (VI.86) and is obvious by analogy. The 
hadronic tensor relevant to inclusive deep inelastic lepton-pion scattering can be obtained 
from the forward Compton process via the optical theorem: 



W,.{q; P) = W'^M; P) + P) 



27r 



-Im 



T+(g;P) + r-(g; P) 



In the Bjorken limit one finds explicitly (Hecht et ai, 2001) 



t 



W+{q-P)=Fnx)t,. + F+{x: 



2x 



Qi^Qu/q ,qj, = qii + 2a;P^, and 



F+{x) = 2xF+{x) , ^1+2(2;) ^ as X ^ 1 



(VI.87) 



(VI.88) 



(VI.89) 



Combining these results with their analogues for W~^, one recovers Bjorken scaling of the 
deep inelastic cross section; namely, the cross section depends only on x, and not separately 
on P ■ q and q^. One may therefore write 



4 1 - 

F2'^{x) = F^{x) + F2 (x) = -[xu{x) + xu{x)] + -[xd{x) + xd{x)] + 

9 9 



(VI.90) 



where the ellipsis denotes contributions from heavier quarks, which are small. We emphasize 
that this result is only valid at lowest order in the strict Bjorken limit. It is not preserved 
under QCD evolution, see Sec. II. D. 

In an explicit calculation of the handbag diagrams. Fig. VI. 39, the resolving scale is 
exposed through the act of taking the imaginary part when employing the optical theorem 
in Eq. (VI.87). At this point a snapshot is taken of the dressed-quark propagators underlined 
in Eq. (VI.86). Their dressing functions are sampled at the particular renormalization point 
and this defines the resolving scale Qq. One can recover the analysis of (Ezawa, 1974), 
augmented by perturbative QCD evolution, by pursuing this route at a resolving scale deep 
into the perturbative domain. 
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On the other hand, (Hecht et ai, 2001) employ algebraic Ansdtze for the elements in 
Eq. (VI. 86) whose form is known even on the nonperturbative domain. They are determined 
by studies of meson properties (Burden et ai, 1996; Maris and Roberts, 1998; Roberts, 1996) 
and exhibit behavior that is broadly consistent with that of QCD's n-point functions; viz., 
the ultraviolet power-laws are explicitly expressed but the logarithms are suppressed in order 
to achieve a level of simplicity. Hence, following the path described above, (Hecht et ai, 
2001) delivers a representation of the pion's distribution functions at an infrared resolving 
scale through the definitions 



How should Eqs. (VI. 91) be understood? In (Hecht et ai, 2001), the computation is in- 
terpreted within the setting of a rainbow-ladder DSE truncation. Hence, sea-quark contribu- 
tions are absent because they cannot appear without nonperturbative dressing of the quark- 
gluon vertex (Chang et ai, 2009; Chang and Roberts, 2009; Cloet and Roberts, 2008). Thus 
Eqs. (VI. 91) describe valence-quark distribution functions and one should have 



viz., the TT"^ contains one, and only one, n- valence-quark and one d- valence-quark. Note that 
one has d'^'^ {x; Qq) = u'^,^ {x; Qq) = d1 (x; Qq) in this calculation. 

It is a deficiency of this and kindred calculations that the model's resolving scale is not 
determined a priori. That can be overcome by calculating the moments of the pion's distribu- 
tion via Eq. (VI. 33), in which case the resolving scale is identical to the renormalization point 
used in computing the dressed Schwinger functions. The cost of that approach, however, 
is a loss of direct knowledge about the distribution's pointwise evolution. In (Hecht et ai, 
2001), Qo was chosen so that the computed distribution, when evolved to Qq = 2 GeV 
using leading-order formulae, produced first and second moments in agreement with those 
reported in (Sutton et ai, 1992a). This yields^^ Qo = 0.54 GeV so that Eqs. (VI. 92) are 
satisfied with a valence-quark mass M = 0.30 GeV~ M{Qo). 

It is apparent that this procedure views the distributions defined in Eqs. (VI. 91) as in- 
frared boundary value input for the valence-quark evolution equations. This is strictly valid 

The quantitative similarity between this and the mass-scale for LO-evolution in Sec. V.C is noteworthy. 



Qo) := 2 Qo) , d^x; Qo) := 9 Ff (x; Qo) • 



(VI.91) 




(VI.92) 
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FIG. VI. 40 (Color online) Solid curve - Pion's valence-quark momentum distribution computed 
in (Hecht et al, 2001). Dashed curve - A(l — x)^^ , with the parameters drawn from Eq. (VI. 97). 

only if all nonperturbative corrections are negligible for Q > Qq, a. constraint whose faithful 
implementation would require Qo > 1 GeV. On the other hand, it might be an efficacious 
approximation so long as a(Qo)^/(4'7r^) ^ 1. In this connection, it is notable that 

^^^^ = 0.017 , Qo = 0.54 GeV; (VI.93) 

whereas a{0.5 Qq)"^ / {An'^) = 0.60. Plainly, there is a narrow domain upon which this per- 
spective may be reasonable. NB. This view mirrors that described in Sec. V.C. Furthermore, 
and significantly, these considerations do not affect the pointwise behavior of the distribution 
function at Qq. 

The valence-quark distribution function computed by (Hecht et al, 2001) is depicted in 
Fig. VI. 40. From this form one finds the momentum-fraction carried by the valence- quarks; 
viz., 

(^.)Qo=o.54Gev = //^ ^ K(^; Qo) + <(x; Qo)] = 0.72 . (VI.94) 

The remainder is carried by the gluons; namely, {xg)'Q^ = 0.28, which are invisible to the 
electromagnetic probe, since the sea distribution is zero at Qq. The pion is a bound-state 
of finite extent formed through the exchange of dressed-gluons between a dressed-quark and 
-antiquark, and hence gluons are necessarily always present. It is therefore plain that at 
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no resolving scale could all the pion's momentum be carried by the valence-quarks. For 
comparison, the parametrized valence-like pion parton distributions in (Gliick et ai, 1999) 
yield a gluon momentum-fraction of^° (^9)^0 ~ ^-^^ ^0 ~ 0.51 GeV. In (Hecht et ai, 
2001), the second and third moments are 



{X. 



u+d)Qo=0.54GeV ~ O-^S , {^u+d)Qo=0.54GeV ~ ^-^l • (VI. 95) 



The low moments are primarily determined by the distribution function's behavior at 
small-x (see the close of Sec. VLB. 2). Sensitivity to the nature of QCD's interaction is 
found at large x, which corresponds to large relative momentum between the dressed-quark 
and -antiquark. Two questions come immediately to mind, namely: what should one expect 
for the x-dependence of the distribution in Fig. VI. 40 at large-x; and what truly constitutes 
the large-x domain? 

The first question is readily answered. In this case Ql ^ 5Aq(^j-), which corresponds to 
a scale whereat the chiral limit mass function is dropping rapidly but does not yet exhibit 
the behavior associated with its truly asymptotic momentum-dependence: Qo = 0.54 GeV 
does not lie beyond the inflexion point of the chiral-limit mass-function, see Fig. 11.11. One 
would therefore anticipate that, for x ~ 1, Ui,{x; Qo = 0.54 GeV) ^ (1 — x)'', with rj ^ 2. It 
was found in (Hecht et ai, 2001) that the distribution in Fig. VI. 40 is pointwise accurately 
interpolated by the function 

x<(x; Qo) = A x"' (1 - €y/^ + 7x) (1 - x)"^^ (VI.96) 

with the fit-parameters taking the values 

11.24 1.43 2.44 2.54 1.90 

These parameters depend on Qq and the value of 772 is fully consistent with expectation. 

The second question posed above can now also be answered quantitatively using 
Eqs. (VI.96) and (VI. 97). The dashed-curve in Fig. VI. 40 is that component of Eq. (VI.96) 
which dominates the x-dependence of m^(x; Qo) at "large x." On the domain 

= {x\x > 0.86} , (VI.98) 



A novel perspective on the magnitude of an hadron's gluon momentum fraction is discussed in (Chen et ai, 
2009a,b; Ji, 2009). 
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FIG. VI. 41 (Color online) Compilation of results for the valence-n-quark distribution in the pion: 
Solid curve - DSE result (Hecht et al, 2001); dot-dashed curve: NJL model (see Sec. VI. A. 6); 
short-dashed curve: instanton model (Dorokhov and Tomio, 2000); dash- dot- dotted curve: light- 
front constituent-quark model (Frederico and Miller, 1994); circles - Drell-Yan data presented 
in (Conway et al, 1989); and long-dashed curve - reanalysis of that Drell-Yan data described 
in (Wijesooriya et al, 2005), which is also depicted in Fig. III. 18. All calculations evolved at 
leading-order to Qo = 4.0 GeV using a four-flavour value of Aqcd = 0.204 GeV, except that from 
(Wijesooriya et al., 2005), which is reported at Qo = 5.2 GeV. 

the dominant component agrees at the level of 20% or better with the full curve. NB. The 
extent of this domain depends weakly on the mass-scale M^: it is a little larger in a model 
with a smaller value of Md; and the disagreement increases to 37% at x = 0.76 

In (Hecht et ai, 2001), using leading-order evolution, the distribution in Fig. VI. 40 is 
evolved to Qo = 2 GeV and Qo = 4 GeV. The w-quark moments at the former scale are 
presented in Table VI. 5 and the curve at the latter scale is presented in Fig. III. 16. 

In Fig. VI. 41 we display a compilation of results for the valence-w-quark distribution in 
the pion. The elucidation herein explains why only the DSE prediction exhibits behavior at 
large-x that is consistent with the QCD parton model, Eq. (VI. 84), first derived in (Ezawa, 
1974; Farrar and Jackson, 1975). The discussion in (Wijesooriya et ai, 2005) and Sec. III. E 
shows that one cannot draw firm conclusions about the large-x behavior of the pion's valence- 
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FIG. VI. 42 (Color online) Solid-curve - Preliminary Dyson- Schwinger equation prediction for the 
ratio of kaon-to-pion valence-quark distribution functions at Qo = 5 GeV, evolved at leading- 
order from the dashed curve, which is the DSE prediction at Qq = 0.57 GeV (Nguyen, 2010; 
Nguyen et al., 2010). Dotted and dash-dotted curves - our re-evaluation of [x) /u^{x) using the 
model in (Shigetani et al., 1993); filled circles - the ratio of kaon-to-pion Drell-Yan cross-sections 
obtained from a sample of dimuon events with invariant mass 4.1 < M < 8.5 GeV (Badier et al., 
1980); thin dashed line - the curve u^{x) = u1{x); and thin dash-dot curve - fit to the data given 
in Eq. (III.24). 

quark distribution function from the single extant vrA^ Drell-Yan experiment (Conway et al., 
1989). The status of QCD as the strong interaction piece of the Standard Model will seriously 
be challenged if an improved experiment, such as that canvassed in (Wijesooriya et al., 2001), 
is also incompatible with Eq. (VI. 84). 

Computations of the valence-quark distribution functions in other mesons are underway 
(Nguyen, 2010; Nguyen et al, 2010) using the DSE approach that successfully predicted the 
pion's electromagnetic form factor (Maris and Tandy, 2000a,b). The impact of the dressed- 
quark mass function on the ratio {x) / u1{x) is depicted in Fig. VI. 42. In comparison with 
the nonpointlike-pion- regularized NJL result in (Shigetani et al., 1993), one finds that the 
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momentum-dependent mass function markedly affects the separate behaviors of u'^(x) and 
u^{x), especially on the valence-quark domain. However, the preliminary indication is that 
it does not materially affect the ratio; e.g., 



We anticipated a nonzero value for the ratio because, in the neighborhood of x = 1, the 
M-quark distribution should have the same pointwise behavior in all pseudoscalar mesons. 
These features and predictions add additional emphasis to the need for a much improved 
measurement of the kaon structure functions. Before claiming understanding of QCD, it 
is crucial to verify the predicted large-x behavior, especially those properties which are 
environment-dependent and those which are not. 

The calculation of baryon valence-quark distribution functions is also possible; e.g., using 
the Poincare covariant Faddeev equation described in App. B. This approach capitalizes on 
the importance of scalar and axial-vector diquark correlations within the nucleon. In QCD 
these correlations are essentially nonpointlike (Alexandrou et ai, 2006; Maris, 2004). While 
no DSE computations of the pointwise behavior of the nucleon structure functions are yet 
available, based on the diquark-correlation probabilities presented in Table 2 of the DSE 
study of nucleon electromagnetic form factors described in (Cloet et ai, 2009) and the fact 
that the ratio is a fixed-point under evolution, one can estimate 



cf. Eqs. (V.5), (VI.6) and (VI.46), and Figs. 11.10, VI.28 and VI.35. In addition, the nonzero 
value highlights the important role of axial-vector diquark correlations: they enable a truly 
valence d-quark to carry the proton's helicity, which is impossible if scalar diquarks are the 
only correlations present in the proton's Faddeev amplitude. 

VII. PERSPECTIVE AND PROSPECTS 

Understanding the physics of hadrons on the valence-quark domain is a definitive task 
for hadron physics. Indeed, a given hadron is defined by its flavor content and that is a 
valence-quark property. Whilst significant continuing and new investments in experiment 
and theory are required in order to acquire this understanding, the potential rewards are 
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great. For example, this information is essential to both: an elucidation of the Standard 
Model; and the provision of reliable parton distribution functions for use in the analysis of 
large collider experiments - our best tool for uncovering phenomena outside this paradigm. 

We have highlighted numerous outstanding and unresolved challenges in the valence 
region, which experiments must confront. Important amongst them are the x > 0.8 region 
of the nucleon and the pion, and the kaon distribution function. The upgraded Jefferson 
Lab facility is well-suited to provide new data for the valence region in the nucleon, from 
both inclusive and semi-inclusive deep inelastic scattering. Moreover, a future electron ion 
collider could provide crucial new data for the mesons. There is also great potential in using 
Drell-Yan interactions, at FNAL, J-PARC and GSI, to push into the large-x domain for 
both mesons and nucleons. 

The challenges for theory are exciting and equally great. The valence-quark domain is 
the purview of nonperturbative methods in QCD. Hence, true understanding, in contrast to 
modeling and parametrization as a simple means of identifying and highlighting key features, 
requires a widespread acceptance of the need to adapt theory: to the lessons learnt already 
from the methods of nonperturbative quantum field theory; and a fuller exploitation of those 
methods. Thus, before an elucidation of hadron structure can be achieved, theory must 
accept the most conspicuous of QCD's challenges. The Standard Model will not otherwise 
be solved. 
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Appendix A: Light-front conventions 



Unless otherwise specified, we use the following light-front notation: 

^+ ^ ^- ^ e) , (A.i) 

where is any quantity that transforms as a contravariant four-vector. Naturally, ^± = 
and f = £+i- -i±-i±. It follows that 

j d^^^\j de"de+dH± . (A.2) 

The spinor describing a spin-i hadron is orthonormalized according to 

A')w(p, A) = 2p+(27r)35,,,5(p' - p)6'{p'^ - p^), (A.3) 

where A, A' are helicities. 

NB. If one considers a particle's momentum, then A;"'" > in the light-front frame. Hence, 
the only way to make a zero momentum Fock state is for every particle it contains to have 
= 0. Although this defines a set of measure zero in the light-front phase space, such 
states cannot casually be discarded because operators exist which are singular at = 0. 
For example, to throw away these states is to preclude dynamical chiral symmetry breaking, 
a keystone of nonperturbative QCD and hadron physics. 

Appendix B: Dressed-quarks, the Faddeev equation and the current 

When emphasizing dressed-quark degrees of freedom it is now natural to describe baryons 
using a Poincare covariant Faddeev equation. This approach sits squarely within the ambit 
of the application of Dyson-Schwinger equations (DSEs) in QCD (Maris and Roberts, 2003; 
Roberts et ai, 2007; Roberts and Schmidt, 2000). 

One arrives at the Faddeev equation by noticing that in quantum field theory a nu- 
cleon appears as a pole in a six-point quark Green function. The pole's residue is pro- 
portional to the nucleon's Faddeev amplitude, which is obtained from a Poincare covariant 
Faddeev equation that sums all possible exchanges and interactions that can take place 
between three dressed-quarks. A tractable Faddeev equation for baryons was formulated 
in Ref. (Cahill et ai, 1989). Depicted in Fig. B. 43, it is founded on the observation that 
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FIG. B.43 Poincare covariant Faddeev equation. ^ is the Faddeev amplitude for a nucleon of total 
momentum P = pq + pa- It expresses the relative momentum correlation between the dressed- 
quark and -diquarks within the nucleon. The shaded region demarcates the kernel of the Faddeev 
equation, in which: the single line denotes the dressed-quark propagator; T is the diquark Bethe- 
Salpeter-like amplitude; and the double line is the diquark propagator. [Full explanation provided 
in App. A of (Cloet et al., 2009), from which this figure is adapted.] 

an interaction which describes colour- singlet mesons also generates quark-quark (diquark) 
correlations in the colour-3 (antitriplet) channel (Cahill et al, 1987). 

The dominant correlations for ground state octet and decuplet baryons are scalar (0"'") 
and axial- vector (1+) diquarks because, for example, the associated mass-scales are smaller 
than the baryons' masses (Burden et al, 1997; Maris, 2002), namely (in GeV) 

m[ud]g+ = 0.7 - 0.8 , rrnuu)^+ = m(_ud)^+ = m{dd)^+ = 0.9 - 1.0 . (B.l) 

While diquarks do not appear in the strong interaction spectrum; e.g., Refs. (Bender et ai, 
2002, 1996; Bhagwat et al, 2004), the attraction between quarks in this channel justifies a 
picture of baryons in which two quarks are always correlated as a colour-3 diquark pseu- 
doparticle, and binding is effected by the iterated exchange of roles between the bystander 
and diquark-participant quarks. 

The kernel of the Faddeev equation is completed by specifying that the quarks are dressed, 
with two of the three dressed-quarks correlated always as a colour-3 diquark. As illustrated 
in Fig. B.43, binding is then effected by the iterated exchange of roles between the bystander 
and diquark-participant quarks. The Faddeev equation yields the nucleon's mass and am- 
plitude. 

In order to calculate nucleon distribution functions one must know the manner in which 
the nucleon described by the Faddeev equation couples to a photon. That is derived in 
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FIG. B.44 Vertex which ensures a conserved current for on-shell nucleons described by the Fad- 
deev amphtudes, ^ij, obtained from the Faddeev equation depicted in Fig. B.43. The single hne 
represents the dressed-quark propagator; the double line, the diquark propagator; and F is the 
diquark Bethe-Salpeter amplitude. The remaining vertices are described briefly in the text: the 
top-left image is Diagram 1; the top-right, Diagram 2; and so on, with the bottom-right image, 
Diagram 6. In connection with DIS, the photon line is equated with a zero-momentum insertion. 
[A full explanation of the diagrams depicted here is provided in App. C of (Cloet et al, 2009), from 
which this figure is adapted.] 

(Oettel et al, 2000) and illustrated in Fig. B.44. Naturally, the current depends on the elec- 
tromagnetic properties of the diquark correlations. A detailed explanation of the diagrams 
in Fig. B.44 is presented in App. C of (Cloet et al, 2009). Here we only provide a brief 
explanation. 

Diagram 1 represents the photon coupling directly to the bystander quark. It is a nec- 
essary condition for current conservation that the quark-photon vertex satisfy the Ward- 
Takahashi identity. Since the quark is dressed, the vertex is not bare. It can be obtained by 
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solving an inhomogeneous Bethe-Salpeter equation. 

Diagram 2 depicts the photon couphng directly to a diquark correlation. Naturally, the 
diquark propagators match the line to which they are attached. Moreover, the interaction 
vertices satisfy Ward-Takahashi identities and obey QCD constraints. 

Diagram 3 shows a photon coupling to the quark that is exchanged as one diquark breaks 
up and another is formed. It is a two-loop diagram. It is noteworthy that the process 
of quark exchange provides the attraction necessary in the Faddeev equation to bind the 
nucleon. It also guarantees that the Faddeev amplitude has the correct antisymmetry under 
the exchange of any two dressed-quarks. This key feature is absent in models with elementary 
(noncomposite) diquarks. 

Diagram 4 differs from Diagram 2 in expressing the contribution to the nucleons' form 
factors owing to an electromagnetically induced transition between scalar and axial-vector 
diquarks. This transition vertex is a rank-2 pseudotensor, kindred, e.g., to the matrix 
element describing the p7*7r° transition (Maris and Tandy, 2002). 

Diagrams 5 & 6 are the so-called "seagull" terms, which appear as partners to Diagram 3 
and arise because binding in the Faddeev equations is effected by the exchange of a dressed- 
quark between nonpointlike diquark correlations (Octtcl et ai, 2000). The new elements in 
these diagrams are the couplings of a photon to two dressed-quarks as they either separate 
from (Diagram 5) or combine to form (Diagram 6) a diquark correlation. As such they 
are components of the five point Schwinger function which describes the coupling of a pho- 
ton to the quark-quark scattering kernel. These terms vanish if the diquark correlation is 
represented by a momentum-independent Bethe-Salpeter-like amplitude; i.e., the diquark is 
pointlike. 
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